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ABSTRACT 


This thesis examines the pool boiling heat transfer 
characteristics of HCFC-124 and HCFC-124/oil mixtures with up 
to 10% (by weight) miscible alkylbenzene oil. One smooth and 
4 enhanced tubes were tested: a 19 and 26 low integral-fin 
tube (GEWA-K); a modified finned tube (TURBO-B); and a porous 
coated tube (HIGH FLUX). The tests were carried out using the 
procedure used for CFC-114 at the same saturation temperature 
Ome 2:2 CC. This allowed for direct comparison of the pool 
boiling heat transfer characteristics between the two 
refrigerants. 

The smooth and GEWA-K 19 fin per inch tube performance in 
pure HCFC-124 Ame HCFC-124/oil mixtures ranged between 10 to 
50% better than in pure CFC-114 and CFC-114/oil mixtures for 
all heat fluxes. The HIGH FLUX and TURBO-B tubes were similar 
in performance. With pure HCFC-124, the finned tubes typically 
provided enhancements in the heat transfer coefficient between 
2and 3 times that of a smooth tube. The HIGH FLUX and TURBO-B 
Surfaces typically provided additional enhancements 2 times 
that of the finned tubes. With the addition of oil, the heat 
transfer increased from the smooth and finned tubes but 
decreased from the HIGH FLUX and TURBO-B tubes. The HIGH FLUX 
and TURBO-B tubes therefore exhibited enhancements less than 
the finned tubes at high oil concentrations and high heat 


fluxes. 
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NOMENCLATURE 


area 

tube outside surface area of active boiling section 
cross-sectional area of tube 

specific heat 

diameter 

tube inside diameter 

outside diameter of unenhanced ends 

diameter at the position of the thermocouple 
diameter of test section to the base of fins 
gravitational acceleration 

heat-transfer coefficient 

latent heat of vaporization 

Gule rent 

current reading by AC Current Sensor 

thermal conductivity of liquid retvigerwane 
thermal conductivity Gf copper: 

active boiling tube length 

non-boiling length of the test tube 

Nusselt number 

tube outside wall perimeter 

Prandtl Number 

heat-transfer rate from boiling surface 
heat-transfer rate through one non-boiling end 
heat-transfer rate from cartridge heater 

heat flux 

Rayleigh number 

temperature 

average wall temperature at the thermocouple location 


Critical temperature 


Atiel 


sat 


= 
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Seon ote See eed 


qs 


film temperature 

temperature of the respective (n) thermocouple 
location 

Saturation temperature 

outer wall temperature of the boiling test tube 
voltage across the cartridge heater 

voltage reading by AC-DC true RMS converter 
thermal diffusivity 

volumetric thermal expansion coefficient 
uncertainty in measurement and calibration 
sSuperheat (Two-Tsat) 

dynamic viscosity 

kinematic viscosity 

density 


surface tension 


5 all 





I. INTRODUCTION 


ini eSonOnCOlig@mmerrOnre. to | eliminate the use of 
chlorofluorocarbons (CFCs), which have been linked to the 
depletion of the earth’s protective ozone layer, the United 
States Navy is looking into using recently developed 
alternative refrigerants. There are 2 principal alternative 
refrigerants: NYawotinerocarbons earCs). and 
MyemociloroLlucorocarbons (HCFCs). HFCs contain no chlorine and 
have zero ozone depletion potential (ODP). HCFCs do contain 
chlorine, but with the attached hydrogen atom to the CFC 
structure, the molecule is less stable and the ozone damaging 
chlorine is released in the lower atmosphere prior to reaching 
the ozone layer. Ties OCOP= of ~PAGPC-124 1s about 23, 
Significantly less than the 70% ODP for CFC-114. Figure 1.1, 
provided by Baehr [Ref. 1], displays the advantages of HCFC- 
124 over CFC-114 with respect to both ozone depletion 
potential (ODP) and global warming potential (GWP). 

The US Navy is currently using CFC-114 (R-114) in all of 
Meeelarger capacity centrifugal air conditioning units (>100 
tons). The search for a ‘drop-in’ replacement has identified 
HCFC-124 (R-124) as a possible alternative to R-114. ‘Drop-in’ 
here means operating at the same evaporating temperature as 
existing systems. Therefore the operating parameters of R-124 


under similar conditions as present R-114 chiller systems 
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must be taken into account. For example, the shell must be 
capable of withstanding higher pressures (11 psig for R-124 at 
fee -C compared to Owpsig for R-114 at 2.2 °C). The shutdown 
system pressure is also a significant factor when the system 
is subjected to higher temperatures from a surrounding engine 
room or perhaps a casualty condition. A temperature of 100 °F 
would result in a system vapor pressure of 31 psig for R-114 
compared to 68 psig for R-124. If the heat exchanger shell 
can handle the increased system pressure, hardware 
modifications should be limited to alterations such as 
compressor impeller size and changes to sealing materials to 
Maintain leak-tight integrity. 

If hardware modifications can be made, then the other 
major consideration is the heat transfer characteristics of R- 
124 compared to those of R-114. This is not such a simple 
comparison as the two fluids may exhibit varying boiling 
characteristics for each type of enhanced surface. [In 
addition, the oil used for each fluid is different since each 
must be miscible in its particular refrigerant and this may 
Significantly affect heat transfer performance. 

To evaluate and compare the heat transfer performance of 
R-124 with that of R-114, the following thesis objectives were 
established: 

1. Install a new fiberglass-strengthened evaporator and 


condenser section into existing apparatus to withstand 
the higher saturation pressure associated with R-124. 


Collect pool boiling heat transfer data of pure R-114 and 
R-114/oil mixtures over a range of increasing and 
decreasing heat fluxes for the newly manufactured tubes 
and check their repeatability with existing R-114 data 
using the newly modified apparatus. 


Modify existing data reduction program to account foie 
thermophysical properties of R-124. 


Collect pool boiling heat transfer data of pure R-124 and 
R-124/oil mixtures over a range of increasing and 
decreasing heat fluxes usSing a smooth tube and four 
enhanced surface tubes. 


Study effect of oil on boiling hysteresis for a smooth 
aNnG DOnROUS COareoa Eube. 


Compare single tube heat transfer characteristics of R- 
124 with previously obtained R-114 data taken by Sugiyama 
(Ref. 2]. 


II. POOL BOILING OF REFRIGERANTS 


A. POOL BOILING CURVE 

The pool boiling curve, as provided by Bar-Cohen [Ref.3], 
is shown in Figure 2.1. The applicable regions within the 
curve where single tube pool boiling heat transfer tests of R- 


114 and R-124 refrigerants were conducted are reviewed. 


Heat flux, q” 





(ATs ), (ATs buen, FRA 
Wall Superbea:. T,-T, 


Figure 2.1 Pool Boiling Curve 


1. Natural Convection 
This region (a-b) occurs during the increéasinguees 
flux phase where there are no vapor bubbles being produced at 
the tube surface. In this region, heat transfer is dueveeueuss 
to single-phase free convection effects. Several correlations 
are available to predict heat transfer performance in this 
region, including the Churchill and Chu (Ref. 4] corredeae 


shown below in equation (1). 


h=k/D,[.6+.387(g°B-O-D,7/u-a) 1/8/ [14+ (0.559/Pr) 2/15] 9/27) 


2. tncapiene, Bolling 

The transition from the natural convection reqiuGaiie. 
the nucleate boiling region is the subject of numerous 
studies. The superheat (T,.31-Tea-) at which bubble ine@epaaias 
(point b) occurs is referred to as the onset of nucleate 
boiling (ONB) or incipient boiling. Considerable reduce renee 
tube wall temperature occurs at this point, resulting ina 
considerable reduction in superheat noted as ATL, in page 
2 elie This ‘jump’ in performance only occurs damage 
increasing heat flux and is referred to as the thermal 
overshoot. A hysteresis effect occurs when decreasing heat 
flux data (from point c to point a) are compared )eomiaaes 
initial ascending path (a-b-c). Several factors believed to 


influence the thermal overshoot at the ONB for highly-wetting 


liquids were studied by Bar-Cohen [Ref. 3]. Probable 
Pifluential factors included: vibration,~ surface tension, 
enhanced commercial surfaces, pressure, time between test 
runs, oil mixture, temperature attained at end of decreasing 
heat flux phase of the cycle, and variation with cool-down 
procedure. Because of these numerous influences on the boiling 
incipience superheat, past studies, including Bar-Cohen [Ref. 
3], reveal highly variable and widely distributed superheats 
at ONB. Furthermore, the process of first nucleation is 
inherently stochastic, leading to variations in behavior from 
one day to the next. 
3. Nucleate Boiling Region 

Once nucleate boiling commences, during increasing 
heat flux (point c), sustained improvement in heat transfer is 
realized throughout the entire fully-developed nucleate 
boiling region (c-d). As vapor bubbles form at nucleation 
sites and separate from the surface, dramatic increases in 
heat transfer are attained. There are several factors 
affecting where this region is located on the pool boiling 
curve, the most important of which is the number of active 
Sites on the surface, as may occur on different boiling 
surfaces. Thome [Ref. 5] provides detailed insight into 
enhanced boiling heat transfer from various commercially 
available surfaces. For plain tubes, several correlations have 


been proposed to predict boiling performance. Two. such 


correlations were provided by Stephan and Abdelsalam [Ref. 6] 
and Chongrungreong and Sauer [Ref. 7]. 
B. SINGLE TUBE BOILING EXPERIMENTS AT 

NAVAL POSTGRADUATE SCHOOL 

Heat transfer measurements were collected by McManus 
(Ref.8] for single tube pool boiling of R-114 at a satiated 
temperature of 13.8°C with up to 6% oil concentration. The 
boiling tubes in these tests were heated internally by warm 
water. HIGH FLUX, TURBO-B and smooth tubes were tested in this 
manner. McManus [Ref. 8] reported enhancements (compared to 
the smooth tube) in pure R-114 as high as 14.6 and 6.4 for the 
HIGH FLUX and TURBO-B tubes respectively, while performance 
decreased to enhancements of 7.0 and 4.9 respectively with a 
6s Oll- CONnGenerar lon: 

Pool boiling heat transfer coefficients in pure R-114 and 
R-114/oil mixtures up to 10% (by weight) were measured by 
Wanniarachchi et al [Ref. 9] for electrically-heated smooth 
and enhanced tubes. These tests were run at a saturation 
temperature of 2.2°C. At 30 kW/m* Wanniarachchi et al [Rememen 
reported enhancement results for the HIGH FLUX, Thermoexcel-E, 
Thermoexcel-HE, and GEWA-T tubes of 9.1, 8.2, 6.8, and 4.4 
respectively. Also reported were the enhancements for these 
tubes with 3%, and 10% oil concentration. 

Sugiyama [Ref. 2] conducted extensive pool boiling tests 
with R-114 and R-114/oil mixtures using a large variety of 


enhanced surface tubes. For the re-entrant cavity tubes with 


pure R-114, Sugiyama [Ref. 2] reported enhancements of 6.4 and 
6.2 for the HIGH FLUX and TURBO-B tubes respectively at 35 
kW/m. With 10% oil these enhancements were reduced to 3.8 and 
5.6 respectively. This degradation in performance was 
partially attributed to the clogging of surface pores with 
Oil. For GEWA-K 26 and 40 fpi tubes, enhancements of 2.4 and 
3.1 were reported with pure R-114 at 35 kW/m. With 10% oil 
these same tubeS were surprisingly reported to have 
enhancements of 3.2 and 4.6 respectively. This indicated an 
improved heat transfer performance with addition of oil in the 
case of finned tubes, believed to be caused by foaming action 
created within the pool that was most prevalent at higher oil 
concentrations (3-10%) and heat fluxes (>30 kW/m-) . 

The effect of oil on the incipience of nucleate boiling 
for pure R-114 and R-114/oil mixtures was studied by Memory 
and Marto [Ref. 10] with smooth, finned and re-entrant cavity 
tubes. They observed a significantly lower heat flux at 
incipience for the HIGH FLUX and TURBO-B tubes compared with 
the smooth and finned tubes in pure R-114. However, consistent 
delay in the boiling incipience of the HIGH FLUX and TURBO-B 
tubes was exhibited in the presence of 3 and 10% oil 
concentrations. At these oil concentrations, peiilang 
incipience for the smooth and finned tubes showed no 
dependence with oil. Also emphasized was the inherent 
randomness and variation connected with temperature overshoot 


associated with boiling incipience. 


C. COMPARISON OF R-114 TO R-124 


1. Property Comparison of R-114/R-124 


Table 1 is provided to show direct comparison of the 


properties of R-114 and R-124 ap T=2e2 8 


Table 1. R-114/R-124 PROPERTY COMPARISON 





Property 







Psat / (kPa) 










p (liq) / (kg/m?) 1526 1429 


(g/m:s) 0.4449 0.3375 
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2. Need for Oil Mixtures 


The oil mixtures used are to simulate the aceeuaw 


environment in refrigeration systems which can collect lube 


Oil within a unit's flooded evaporator to concentrations as 


6 


feegm@eects 910s. The O1lS used for R-124 and R-114 are both 
miscible within their respective refrigerants. However, with 
R-114, a York-C mineral oil was used, while with R-124, an 
alkylbenzene oil was used. 
3. Need for R-124 Data 

There 1s presently no data on single tube pool boiling 
heat transfer with R-124 or R-124/oil mixtures. Tests must 
therefore be conducted with R-124 to parallel R-114 results 


obtained to date. 
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III. DESCRIPTION OF EXPERIMENTAL APPARATUS 


A. SYSTEM DESCRIPTION 

Details of the relocation and early modifications to the 
Original apparatus are given by Sugiyama [{Ref. 2] and Lake 
[Ref. 11]. To accommodate the increased operating pressures 
associated with R-124, stronger evaporator and condenser units 
were installed. The data reduction program DRP71 was modified 
(DRPGB) to include the thermophysical properties of R-124. 

A schematic diagram representing the experimental 
apparatus is shown in Figure 3.1 and comprises the following 
components: 1) an evaporator tee used to boil the refrigerant; 
(2) a condenser tee used to condense the refrigerant; (3) a 
reservoir for liquid refrigerant; (4) a refrigerant/oil 
subsystem; (5) a refrigeration cooling support system; (6) a 
vacuum pump; (7) a data acquisition and instrumentation 
System; (8) an aluminum/plexiglas framework housing items (1) 
Co 14 )2 

The general operation of the system was closed loop with 
vapor generated in the evaporator by the heated tube located 
within the refrigerant pool cf ~ne evapecarar The vapor 
passed through an aluminum ‘L-shaped’ tube to the condenser 


where the vapor was condensed by a chilled water/ethylene 
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Figure 3.1 Schematic of Apparatus 
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glycol mixture circulated through the condenser coil by ieee 
gpm turbine type positive displacement pump. The R-114 or R- 
124 condensate was then gravity fed back to the evaporator 
through a return line. R-502 and R-12 refrigeration systems 
maintained the water/ethylene glycol mixture in the sump 
between -12 and -18 degrees Celsius. Sump cooling was 
accomplished using a counter flow heat exchanger by 
recirculating the sump coolant (using an identical p@seumeae 
displacement pump) over the R-502 side of the heat exchanger. 
To further maintain this coolant temperature, the R-12 system 
was located directly within the water/ethylene glycol sump. 

The addition of oil to the refrigerant withiiieae 
evaporator was accomplished by using a graduated cylinder to 
measure the precise quantity. This simply drained into th 
evaporator by gravity. 

Accurate horizontal positioning of the boiling tubewwas 
critical and accomplished by supporting each end of the tube 
by Teflon inserts. The tube was sealed uSing O-rings. For 
safety, a relief valve on the aluminum pipe between the 


evaporator and condenser was set at a gage pressure of 49 psi. 


B. BOILING TEST SECTION 
1. Evaporator 
The evaporator was a Pyrex glass 'T’ section coated 
with a continuous winding of fiberglass filaments impregnated 


with a modified epoxy resin. This upgraded evaporator shown in 
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Figure 3.2 was designed for a working pressure of 50 psi gage 
and initially hydrostatically safety tested to 75 psi gage. 
The original glass was designed for a working pressure of 30 
psi gage. The inside surface was the same as before and 
maintained the advantages of minimizing nucleate boiling at 
the surface as well as being corrosion resistant. The two ends 
of the evaporator were fitted with a cast-iron flange and 
gasket assembly as shown in Figure 3.2. They mated with 
aluminum flanges at each end of the evaporator, which in turn 
contained the thermocouple housings, oil entry connection and 
Teflon inserts containing the evaporator tube. 
2. Evaporator tubes 

Figure 3.3 shows a schematic of a boiling tube. These were 
accurately held in place using Teflon inserts and sealed by O- 
rings. A 4-stud Sema eee secured the Teflon inserts inside 
the aluminum flanges. 

Five different tubes were tested within the apparatus. 
Dimensions of each tube are provided in Table 5 of Appendix C 
while an expanded view of the enhanced surfaces are shown in 


Figure 3.4. The tubes used were: 


im omooth tube 

2. GEWA-K (19 fpi) finned tube 

3. GEWA-K (26 fpi) finned tube 

4. HIGH FLUX (porous coated) enhanced tube 


5. TURBO-B (deformed surface) enhanced tube 
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a) Corning Pyrex Glass Evaporator (D x Ri = 402x510 


L = 178 mm, ly= 127 mn) 





INSERT 


b) Cast Iron Flange and Gasket (d,= 190mm, dg= 210mm, 


L,= lumm, A = Piles) 


Figure 3.2 Sketch of Pyrex Glass Vessels and Flanges 
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Schematic of boiling tube 


Figure 3.3 
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TURBO-B SURFACE 





HIGH FLUX SURFACE 


Figure 3.4 Surface Profiles of Enhanced Tubes 
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The heater used within each boiling tube was nominally a 
1 kW, 240-Volt stainless-steel cartridge heater, measuring 
6.35 mm in outside diameter and 203.2 mm in length. The actual 
heated length of the heater was 190 mm. The HIGH FLUX and 
GEWA-K 26 fpi tubes were previously manufactured while the 
TURBO-B, GEWA-K 19 fpi and smooth tubes required assembly 
@eror tO testing with R-124 or R-114. The manufacturing 
process of these tubes is described below. 

The heater was placed inside a copper sleeve (Figure 3.5 
provides dimensions of smooth tube sleeve) with a tight 
mechanical fit of .001 inch clearance. The copper sleeve was 
meaenhined with 6 grooves (1.27mm x 1.27mm) on its outer 
diameter of various lengths. Thermocouples were laid in each 
of these grooves and the edges of the grooves were peened over 
with a punch. Any resulting imperfections were smoothed out 
with 250 grit emery paper. The thermocouple wires and heater 
leads all exited the copper sleeve at the same end. The copper 
sleeve containing the heater and thermocouples was then placed 
on 2 aluminum V-blocks and the outer surface of the sleeve was 
covered with a flux and solder paste (50/50 tin/lead); great 
care was taken to ensure complete coverage of the thermocouple 
wires lying in the grooves. The thermocouple and heater leads 
were then passed through the enhanced boiling tube (the outer 
diameter of the copper sleeve was machined to fit tightly into 
each respective tube). To prevent shorting between the heater 


and outer tube shell, a high temperature RTV silicone rubber 
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Figure 3.5 Schematic of copper sleeve 
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mas applied to the heater and thermocouple leads at the point 
where they leave the sleeve. The copper sleeve was then 
tapped gently into the tube and positioned accurately so as to 
line up precisely with the 8 inch active length of the 
enhanced surface. This enhanced surface length was different 
from the actual heated length which was only 7.5 inches. 

The entire assembly was then placed horizontally in an 
oven and heated to between 400°F and SOO°F. This was the 
temperature ‘window’ at which the flux melted, but the 
thermocouple leads did not. The temperature of the tube was 
monitored by attaching the thermocouples to ae simple 
temperature indicator. After stabilizing the oven temperature 
between these temperatures and ensuring that the solder had 
melted, the tube was taken out and quenched with a damp cloth 
at one end (as close to the end as possible) where the 
thermocouple and heater leads protrude. The oven was then 
placed vertically and the tube reinserted from below keeping 
the lower end quenched while holding both tube and cloth with 
a pair of channel locks. When the temperature increased to 
>400°F, the solder again began to melt. Additional 40/60 tin- 
lead solder was then applied by hand from the end which 
protrudes from the top of the oven. Remelting and addition of 
extra solder in this way ensured a uniform circumferential 
layer of solder and helped to prevent formation of air 
pockets. The tube was then removed from the oven and quenched 


over its whole length with a damp cloth while in a vertical 


waa 


position. It was then left to cool in air for 24 Noumeeeiaeee 
procedure differed from earlier manufactured tubes in tie 
previous copper sleeves had 4 thermocouples leaving each end 
of the sleeve. Trying to repeat this previous pre@eemem. 
accurately led to large discrepancies in the themmi@e@aaiea. 
readings (220°C) as heater power was increased. It was felt 
this was due to air pockets created in the fabrication 
procedure (the tube could never be placed vertically ise 
oven with leads protruding from both ends). The above 
procedure was much simpler and led to excellent uniformity in 


the thermocouple values at all heater powers. 


C. CONDENSER SECTION 

The condenser was also a fiberglass coated Corning Pyrex 
glass Tee identical in size and shape to the evaporator shown 
in Figure 3.2. The refrigerant was condensed on a 3/8 inch 
outside diameter copper tube that was coiled within the 
condenser. The condenser was mounted vertically. An O-rin 
sealed aluminum flange mounted at the top of the condenser 
contained not only the water ethylene/glycol coolant 
inlet/outlet, but also a 3/8 inch outside diameter copper aa 
vent, which was connected either to a vacuum pump via valve 
V1ll or to a charging cylinder via valve V10. The bottonoreer 
condenser was also capped with an O-ring sealed aluminum 
flange and contained the gravity feed return line to the 


evaporator via valve VS. The refrigerant vapor ftroneeme 


Ze 


evaporator entered the condenser through an aluminum L-shaped 
tube, connected to the condenser by O-ring and gasketed 
flange. A pressure gage ranging from 30 inches mercury vacuum 
to 150 psig was mounted on the L-shaped tube; this was used 


only as a rough visual check of pressure within the facility. 


D. OIL ADDITION SECTION 

Heme provide Gil igaddition into the «mefrigerant, a 
cylindrical aluminum reservoir, 6 inches in diameter and 6 
inches high, and a graduated cylinder 355 mm in length and 
25.4 mm in diameter is used as shown in Figure 3.1. The oil 
used with R-124 was a miscible alkylbenzene oil while the oil 
Meeqd with R-114 was a York-C mineral oil, also miscible. The 
graduated cylinder was filled from the reservoir via valve V2, 
while the reservoir was replenished through valve V3 at the 
top. Oil entered the evaporator from the graduated cylinder by 
Gravity feed via valve V1. Oil concentrations (by weight) of 


0,1,2,3,6, and 10% were used for each of the tubes. 


E. COOLING SYSTEM 
1. Water/Ethylene-Glycol Coolant Tank 
At the base of the apparatus, a 5.35 cubic foot tank 
stored Se gallons of ethylene-glycol/water mixture. 
Recirculation and discharge ports were connected to the tank 


together with a thermocouple to monitor sump temperature. 


ne 


2. Refrigeration Plants 

To provide the cooling of the water/ethylene glycol 
mixture, a 1/2 ton R-502 and 1/4 ton R-12 refrigeration system 
were installed near the sump as seen in Figure 3.1. Both 
plants contained a hermetically sealed compressor assembly, 
an-air cooled condenser, receiver, filter-dryer, pressure 
regulator, temperature control switch, and thermoesteaua 
expansion valve. The R-502 plant (Figure 3.6) used a counter 
current heat exchanger while the R-12 plant used a coiled 
evaporator placed within the sump. The sump temperature was 
maintained at approximately -17 °C when both these systems 
were used. Each plant was controlled by a temperature control 
Switch and thermostatic expansion valve. 

3. Pump and Control Valve 

An 8 gpm, positive displacement pump was installed 
with a 1 inch diameter PVC pipe connected to the 
water/ethylene glycol sump. The discharge of the pump was 
piped to the condenser via valve VC as shown in Figure 3.1. A 
by-pass valve V9 was placed upstream of the control valve VC 
to avoid overloading the pump when control valve VC was closed 
and provides for mixing of the return coolant from the 
condenser. At high heat fluxes, the by-pass valve needed to be 
closed down to increase flow of coolant through the condenser. 
An identical pump discharged the ethylene glycol/water coolant 


through the R-502 counter-current heatVexetamgeanue 
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R-502 Refrigerant plant 


Bugure 3.6 


F. REFRIGERANT RESERVOIR 

The reservoir used to store the R-114 or R-124 was made 
from an aluminum cylinder, 9 inches in diameter and 10 inches 
in height. A sight glass attached on the reservoir allowed for 
the liquid level of refrigerant to be monitored. The reservoir 
was connected to the ‘'L’ shaped tube through valve V7 and to 
the liquid line into the evaporator through valve V6 as shown 


Lie Serres. 4, 


G. FRAME 

To house the entire assembly, an aluminum frame was 
constructed measuring 1070.1 mm x 509.8 mm x 610.1 mm. Witha 
view to keeping the apparatus in an isothermal atmosphere, as 
well as for safety considerations, the sides were covered with 
12.7 mm thick plexiglas with hinges mounted on the 2 sidesmas 
provide easy access to the Reser ie The top of the frame was 
covered by plywood while the bottom was covered by an aluminum 
plate. The entire frame rested on the coolant sump tank. The 
front side of plexiglas had small holes to accommodate the 
valve stems of valves V1 through V8, so they can be operated 


while maintaining the integrity of the enclosed frame. 


ear INSTRUMENTATION 
1. Power Measurement 
The boiling tube heater was powered by a 240 volt AC 
source, adjusted to 0-220 volts and 0-5 amps by a variac 


control. Power input to the heater was measured with an AC 
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Data acquisition/control unit 
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Current inductive sensor and a voltage sensor. The output was 
a proportional signal in volts having been converted thvemed 
an AC-DC true R.M.S. converter. Both the AC current sensa@ ages 
the AC-DC true R.M.S. converter were connected to the data 
acquisition/control unit as shown in Piguize =e 
2. Temperature Measurement 

The primary means of obtaining temperatures was the 
use of 30 gage copper-constantan thermocouples. Specie 
temperatures recorded were from 6 thermocouples in the bowie 
tube (8 in the earlier manufactured tubes), 3 thermocouples in 
the evaporator pool (housed within separate sealed 
thermocouple wells within the evaporator pool), and one 
thermocouple to monitor the water/ethylene glycol sump tank. 
A Hewlett-Packard 3497A data acquisition system controlled by 
a Hewlett-Packard 9826 computer read these thermocouple 


measurements. 


I. DATA ACQUISITION AND REDUCTION PROGRAM 

Data collection and reduction program DRP71, used by 
Sugiyama [Ref.2] for use with R-114 was modified for use with 
R-124. The thermophysical properties of saturated R-124 were 
obtained from RefProp software [Ref. 12] and were used to 
generate equations for use in the calculations of the 
reduction program. Appendix A contains plots of Gach prepemay 
versus temperature and the respective equation. The modified 


program was renamed DRPGB. 
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The Hewlett-Packard 9826 computer reduced all the 
thermocouple, current, and R.M.S. voltage readings using the 
DRPGB program. The program was run through the keyboard 
allowing regulation of the data acquisition unit. A printout 
of the reduced data was provided by a Hewlett-Packard Inkjet 
printer. Table 2 lists the channel assignments for the sensor 


inputs to the data acquisition system. 


J. PROCEDURE FROM KEYBOARD 
The following procedure outlines the procedure used to 
collect and process data from the Hewlett Packard data 


acquisition system: 


a Selleacred = ene taking data’ option after initial 
loading of program DRPGB. 

ae Selected the refrigerant being tested, tube heating 
mode, (1e electrically or water heated) and 


thermocouple type from DRPGB. 


a. Assigned name to file of data to be processed and 
stored. 
4. Identified and entered any defective thermocouples by 


channel # 1-8, noting that the new tubes have only 6 
thermocouples and channels 7 & 8 were always 


‘defective’. 
ae Selected test tube type. 
6. SemGceteccmsatiration temperature (2.22°C) “for 


present tests. 


a. Set desired heat fiux desired using variac adjustment. 

a. Attained desired saturation temperature by adjusting 
Plow @mecoolant through condenser coils with control 
valve VC. 
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Once saturation temperature was achieved, waited 5 
minutes for steady state conditions prior to taking 
data. 


Prompted data acquisition unit to scan all the 
channels listed in Table 2. All channel readings were 
made in volts and stored in user specified fields. 

Calculated temperature and power from these voltages. 


Calculated heat transfer rate from cartridge heater. 


Calculated the average wall temperature of the test 
tube and calculated wall superheat (Twall- Tsat). 


Calculated physical properties of R-124 using property 
correlations from RefProp [Ref.12] at film temperature 
(Twall+Tsat) /2. 


Calculated the natural convection heat-transfer 
coefficient of R-124 from the unheated ends Gime 
test tube. 


Computed heat losses from unheated ends of the test 
tube. 


Calculated heat flux from heated length of the test 
tube. 


Calculated the heat-transfer coefficient from the 
heated length of the test tube. 


Stored heat flux and wall superheat for each data set 
in user specified fields. 


Plotted data on available software. 


Sample calculations of the above procedure are given in 


Appendix C. 
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TABLE 2. HP 3497A CHANNEL ASSIGNMENTS 
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IV. EXPERIMENTAL PROCEDURE 


A. ASSEMBLY PREPARATION 
1. Vacuum test of the apparatus 

After placing the desired boiler tube in position 
within the evaporator, the system was evacuated down to a 
pressure of 29 inch Hg. vacuum using a portable vacuum pump 
Via valves V1l and V8 in Figure 3.1. The system was then 
isolated for approximately 30 minutes to check for leak-tight 
integrity (indicated by any significant vacuum loss). If a 
leak was detected, a positive air pressure (20 pSig) was 
placed on the system via valves V10 and V8. A soap/water 
solution was applied on the sealing areas to identify any 


leaks (if there is a leak, bubbles are seen). If identified, 


the leak source waS corrected and the system integrity 
retested. 
2. Charging apparatus with refrigerant 


After system integrity was ensured, the evaporator was 
filled with refrigerant to a level of 20 mm above the test 
tube. The refrigerant was drained by gravity from the 
reservoir into the evaporator through valves V6 and V5. Valve 
V7 was also open to equalize pressure in the system. The 
reservoir was then isolated from the system with the closure 


of valves V5 through V7. 
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3. Degassing and Data Acquisition Channel Check 

Prior to connecting all thermocouple leads, the 
resistance between the heater leads and outer tube wall was 
checked with a simple multimeter. Readings were usually 
greater than 1 MQ . Each channel’s output was checked using 
program SETUP71. Any problem thermocouples were either fixed 
or excluded from the calculations in program DRPGB. Before 
energizing the heater and starting the tests, the evaporator 
pool was slowly cooled to a saturation pressure of 2.22 °C by 
allowing a low coolant flow through the condenser. Once 
completed, the heating surface flux was set to approximately 
90 kW/m* to degas the refrigerant; any non-condensible gases 
were collected in the condenser and vented off using the 
vacuum pump for approximately one minute. The apparatus was 
then secured for a period of approximately 8 hours prior to 
Operation to allow the surface to cool and become fully wetted 


in the refrigerant. 


B. OPERATIONAL PROCEDURE 
The following procedure was used to operate the facility 
(this follows closely the procedure used by Sugiyama [Ref.2]): 
1. The R-502 and R-12 refrigeration units were utilized 
until the temperature of the water-ethylene glycol sump 
reached a minimum of -8 °C before starting tests (with 


Comeimued Cooling to -17 °C). 


2. The data acquisition unit, computer and variac control 
panel were switched on. 
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The computer program SETUP was loaded and run. All=@are 
acquisition channels were rechecked. The refrigerant 
pool temperature was then cooled down slowly to 2.2 °c 
by circulating coolant through the eoendemeamw 


The data acquisition/reduction program DRPGB was loaded 
and run. 


The desired heat flux setting for the test tube was 
input to the program DRPGB and the variac was adjusted 
to attain that heat flux. Starting heat flisgaae 
between 500 and 600 W/m, 


Control valve VC was adjusted to control coolantwimae. 
through the condenser to maintain a constant saturation 
temperature Of 2.22 7G cn nee 


At each heat flux setting, the thermocouple values were 
allowed to stabilize for about 5 minutes prior to 
taking data. 


The following data were then taken and used in the 
acquisition/reduction program DRPGB: pool temperature; 
test tube wall temperatures; sump temperature; current 
and voltage sensor readings. 


Data were taken at each desired heat flux setting and 
always at the same saturation temperature (the coolant 
flow in the condenser was adjusted to ensure this). At 
each heat flux setting, program DRPGB generated a 
permanent printout listing: wall temperatures of the 
test tubes; refrigerant pool temperatures; sump 
temperature; wall superheat; heat-transfer coefficient 
and the corrected heat flux (lie. to account for heat- 
transfer lost from each end). 


For each heat flux setting, steps 6 through 9 were 
repeated. The superheat (Twall-Tsat) and the respective 
heat flux were plotted on a log-log scale. 


Table 3 provides a summarized listing of all data runs. 
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TABLE 3. LISTING OF DATA RUNS . 


Data File Tube — Purpose 
} DO821SM0 Pepe aca ane v4 
| DO824SM0 AS OlSVerieeVOrl Ik yssy 
| DioissMo 


/pi217sMo. R - Se Nercocn | eee bet 
D1016SM4.44"~ R-124 Daccesaes 44 6 
pi114smo.sd R- = a ‘eee ON Ho ey 


| DO111ASM10 R-124 Smooth Data 
po1isM10__—| R-124 


25 


Gon. tinudt on. of. Tabwvewsm 

Data File 
po220SMi0 pata | 
D0109ASMiO Data 


| Do109smi0 Data 
Data 
p11242610 Data 


eee R-114 TURBO-B _—si|: Data 
D1006TBO R-114 TURBO-B Repeatability 
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Continuation of Table 3. 


Purpose 
DL002GYX0 Data 
Peer OTEO R-124 Data | 
Data 

Data 
R124 Data 


Peet 2 SH 6 Ro 24 HIGH FLUX Daiwa. 





oy 


Continuation Of ‘Tab le....3.. 


pata Pile 

130HF10 
HIGH 
| D012 7CHF20 
POLZ6GEHE PO Rawle 2 AiG Hae ax Dave a 


The tabulated data was filed using the following file name 
system: 












Example (DO126AHF10) : 
1. First letter D siamply refers terdacs 


2. The following 4 characters are always the date; in tags 
case Jan 26. 


3. The next letter represents order of the day test was 
completed on the respective date; 
no letter sirore cud 
A- 2nd Yun 
B- 3rd run 
Cee 


4. The following 2 letters represented the tube used; 
SM= smooth 
19= 19 fpi GEWA-K 
26= 26 fpi GEWA-K 
HF= HIGH FLUX 
TB= TUReBO-B 


5. The last 1 or 2 numbers represent the percentagemas 
Oil; in this: case Ge 


* this one run was pure R-124 tested at 4.4 °C vice 2.2 °C 
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V. RESULTS AND DISCUSSION 


Pe REPRODUCTION OF EXPERIMENTAL DATA 

The first tests conducted were to verify the 
reproducibility of the single tube data taken by Sugiyama 
[Ref. 2] with R-114. A GEWA-K 26 fin per inch (fpi) tube 
comparison with R-114 is given in Figure 5.1 and shows very 
Similar results with previous data. Uncertainty bands for both 
wall superheat and heat flux are given on the figure. 
Consequently, at low heat flux (< 1000 W/m“) , there 1S more 
scatter due to the lower wall superheats and hence higher 
uncertainty. Overall, the agreement between the two data sets 
1S satisfactory. 

In Figure 5.2, an increasing heat flux data run using 
a smooth tube is compared with Sugiyama’s [Ref. 2] data. It 
Can be seen that in the low heat flux, natural convection 
region and the high heat flux, nucleate boiling region the 
data agreement is good. In the mid-range of heat flux, (3 
kW/m*< Ges kW/m*), the data varies somewhat, depending upon 
how the tube nucleates. Curve #1 nucleates at a relatively low 
Mel superheat (12.5 °C) while curve #2 (both taken by 
Sugiyama [Ref.1]) nucleates at a high wall superheat of about 


35 °C (point C). The present data (curve #3) exhibits partial 
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Figure 5.1 
26 fpi tube (decreasing flux) 
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nucleation (see point A) where 3 of the six thermocouples had 
a low wall superheat (< 14 °C) and the remaining 3 had a wall 
superheat (> 18 °C). Observing the tube surface, it could be 
Clearly seen that half the tube was boiling while the other 
half was in the natural convection mode. At point B, alaieues 
thermocouples were within 0.9 °C of each other, indicating 
complete nucleation over the entire surface. At this point 
(and above), agreement with the previous data is excellent. 
The discrepancies in this mid heat flux region demonstrate the 
random behavior of nucleate boiling incipience; this is 
further investigated later in this chapter. After all sites 
nucleated (Point B), the standard deviation was 0.9 °C with an 
average thermocouple temperature of 12.9°C. The point of 
incipience and its random occurrence at different levels of 
Superheat will be examined later in this chapter. 

Sugiyama [Ref.2] did not conduct tests on a 19 fpi GaWee 
K tube. This is one of the most common tubes used in 
commercial evaporators. Therefore, base-line data for the 19 
fol GEWA-K tube with R-114 and R-114/oil mixtures were taken 
during this thesis before evacuating and cleaning the 
apparatus, and filling with R-124. The R-114 data obtained was 
incorporated in the R-124/R-114 comparison of the 19 fpi tube 
(2s akle auiparcts Gy Ibe 

Figure 5.3 shows the repeatability of decreasing heat 
flux data for the newly constructed 19 fpi GEWA-K tube using 


R-124. The good agreement is an indication of the consistency 
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Figure 5.3 Repeatability Comparison of R-124 Data with 0% Oil 
for 19 fpi GEWA-K Tube (decreasing flux) 
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of the data reduction program DRPGB, the apparatus and the 
operator. Decreasing heat flux runs provide better comparison 


as temperature overshoot and hysteresis effects are avoided. 


B. POOL BOILING HEAT-TRANSFER COEFFICIENT CORRELATIONS 

Among the correlations available to predict pool boiling 
heat transfer coefficients for a smooth tube, most can be 
expressed as heat flux as a function of wall superheat. The 
pure R-124 smooth tube data is fitted very closely Gysaame 
equation: 

G=1 7 -0en ne ‘ae 

A correlation in terms of fluid physical properties and 
tube geometry is given by Chongrungreong and Sauer [Ref. 7]. 
For R-124 at 177 kPa, this reduces to: 

q-45) onan. (2) 

Another simpler correlation based only on system pressure 

is also given by Chongrungreong and Sauer [Ref. 7]: 
G=76 9 nome (3) 

Stephan and Abdelsalem [Ref. 6] present a statistical 
correlation for four different groups of ftluads? jae. 
refrigerants, their correlation 1s presented in the form: 

q= (ClaT)?:?4 (4) 
where Cl is given as a function of refrigerant type and system 
pressure. Unfortunately, as R-124 is a relatively new 


refrigerant, their graphical representation of Cl does not 
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ia@elucesR-124. However, Ehe present data could be fitted to an 
equation of the form given in equation (4): 

oe le cyclayat)) 2 = (5) 
memee, at a pressure of 1 atm~e Cl> 1.34. 

Figure 5.4 shows the present R-124 smooth tube data for 
decreasing heat flux along with the above correlations. At 
lower heat fluxes the prediction of the Chongrungreong and 
Sauer (Ref. 7] correlation (equation (2) above) exhibits 
considerably closer prediction than the other correlations. In 
the more complex nucleate boiling region, the Stephan and 
Abdelsalem [Ref. 6] correlation and the Chongrungreong and 
Sauer [Ref. 7] correlation (equation (3)) display a fair 


prediction of the slope. 


C. SMOOTH TUBE DATA IN REFRIGERANT/OIL MIXTURES 

Although tests were conducted with o11 concentrations of 
0,1,2,3,6, & 10%, only data with 0,3 and 10% have been plotted 
mmemost Of the figures for clarity. Figures 5.5 and 5.6 show 
the heat transfer performance of pure R-124, R-124/3% oil, and 
R-124/10% oil for a smooth tube during increasing and 
decreasing heat flux respectively. It can be seen that at the 
highest heat fluxes, as oil concentration is increased to 3%, 
the heat transfer is enhanced by a maximum of 17% over pure R- 
[meee ror ELurther increases in Oil concentrations to 103%, 
performance drops off. This drop-off in thermal performance 


has also been observed by Jensen and Jackman [Ref. 13] who 
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Figure 5.4 Comparison of Smooth Tube Performance for 
Pure R-124 with Prediction 
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Figure 5.5 Performance Comparison for Boiling R-124/ 0%, 3%, 


& 10% Mixtures for Smooth Tube (increasing flux) 
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Figure 5.6 Performance Comparison for Boiling R-124/ 0%, 3%, 


& 10% Mixtures for Smooth Tube (decreasing flux) 
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Pepe uputeame the “decline to @ steeper oil concentration 
@eadeenc, Inhibiting bubble growth. 

When o11 is mixed with R-124, two competing mechanisms 
@eeur: (i!) O11! accumulates in the vicinity of the vapor bubble 
surface as the more volatile refrigerant 1S evaporated; 
Meee aditfustonm Of this “accumulated oil back into the bulk 
mixture establishes an equilibrium concentration gradient 
near the departing bubble. Mechanisms (1) and (2) should act 
to inhibit bubble growth and lead to an overall decrease 
in the heat transfer with oil addition compared with pure 
refrigerant. It therefore appears that there 1s an additional 
mechanism that has a significant impact on the thermal 
performance of pool boiling of refrigerants with small 
Quantities of oil. For R-124/oil mixtures of 3% and greater, 
Sonificant foaming in the evaporator was observed at mid to 
Peon heat fluxes (= 10000 kW/m?). It is possible that a third 
mechanism which could be affecting heat transfer is this 
foaming that occurred within the bulk R-124/oil mixture. 
Udomboresuwan and Mesler (Ref.14] reported enhancement in pool 
MOlling in the presence of foam. They attributed this 
enhancement to the foam causing the liquid/vapor interface to 
be closer to the tube wall. This would create a thin liguid 
film between a vapor bubble and the heated tube surface, 
thereby providing large heat-transfer coefficients. Higher 
heat fluxes and greater oil concentrations tend to increase 


the foaming, which could explain why this enhancement is seen 
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more predominately at high heat fluxes. However, as oil 
concentration increases, eventually the oil concentration 
gradient becomes so steep that mechanisms (1) and (2) tend to 
dominate and heat transfer performance drops off (i126 .7§@me 
foaming enhancement effect is offset). 

At low heat fluxes, the scatter in the data is attributed 
to uncertainty. For the increasing heat flux data (Fig 5 95)% 
1t appears that the addition of oil tends to delay the wjegmam 
at which boiling commences. Due to the statistical nature of 
incipience, however, one should be wary of drawing much 
conclusion from this figure. The effect of o11 on inéi#emeaee 
1s discussed in greater detail later in this chapter. 

Figure 5.7 compares smooth tube decreasing heat flux data 
from this thesis for R-124 with Sugiyama [Ref. 2] R-114 data 
at oil concentrations of 0%, 3%, and 10%. At all levelow@en 
heat flux, the R-124 data consistently demonstrate improvement 
in heat transfer over R-114. This ranges from 50% at low heat 


flux to 15% at high heat flux. The heat flux from aeoeniaes 


surface is proportional to bubble diameter (d,) ) ssaiaeee 
departure frequency (f,), number of active nucleation jaa 
per unit area, vapor density (po), and latent heat of 
VapOrilZae lon (Neg) . Barthau (Ref. 15] concluded, aemen 


determination of d,, f,, and active nucleation site dena 
for a smooth copper tube in R-114 using optical techniques, 


that as saturation pressure was increased, d, decreased, while 
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f,, nucleation site density and heat transfer inereascaduaaeems 
conducting the tests at 2.2 °C the corresponding saturation 
pressure of 11 psig for R-124 is greater than the pressure of 
O psig for R-114. This suggests that the nucleation. 
density was greater for R-124 than ‘for R-114 at the same 
Saturation temperature. Added to this is the higher latent 
heat of vaporization, vapor density, and thermal conduceimaae 
for R-124 at this saturation temperature. The R-114/R-124 
property comparison at 2.2 °C saturation temperature is 
Summarized in Table 1. It is therefore reasonable to expect a 
higher heat transfer coefficient for R-124 at the tested 


Saturation temperature of 2.2 °C. 


D. BOILING PERFORMANCE OF 19 AND 26 FPI GEWA-K TUBES 
Figures 5.8 ard 5.9 show the heat transfer performance of 
a 19 fpi GEWA-K tube in pure R-124, R-124/3% oil, and R- 
124/10% oil during increasing and decreasing heat flux 
respectively. There is a consistent increase in thermal 
performance for oil concentrations of 3% and 10% over pure R- 
124 at high heat fluxes. This increase 1s substantially 
greater than that obtained with the GEWA-K 26 fpi tube for the 
same oil mixtures, as shown in Figures 5.10 and 5.11 50seu 
increasing and decreasing flux respectively. Surprisingly for 
both tubes, there seems to be a significant influence of oil 


in the natural convection region (Figs 5.8 & 5.10) ime 
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Figure 5.8 Performance Comparison for Boiling R-124/ 0%,3%,& 
10% Mixtures for 19 fpi GEWA-K Tube (increasing flux) 
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Figure 5.9 Performance Comparison for Boiling R-124/0%,3%,10% 
Mixtures for 19 fpi GEWA-K Tube (decreasing flux) 
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Figure 5.10 Performance Comparison for Boiling R-124/ 0%,3%, 
& 10% Oil Mixtures for 26 fpi GEWA-K Tube (increasing flux) 
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Figure 5.11 Performance Comparison for Boiling R-124/ 0%,3%, 
& 10% Oil Mixtures for 26 fpi GEWA-K Tube (decreasing flux) 
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tends to be a decrease in performance for increasing oil 
concentration. This was not so obvious with Sugiyama [Ref. 2] 
R-114 data. Once again, the point of incipience seems to be 
slightly delayed with increasing oil percentage. The thermal 
performance of the finned tube is better than the smooth tube, 
evidenced by thermal enhancements of 2.1 to 2.9 with the 19 
men tube and 1.3 to 2.1 with the 26 fpi1 tube. However, it 
should be noted that the heat flux calculation does not 
account for the increased surface area of the finned tube, but 
rather uses a root diameter. Some of the enhancement is 
therefore due to the increased surface area of the fins (Figs. 
5§.24-5.26 compare heat transfer coefficients with smooth 
tube), with each fin behaving essentially as the copper 
surface of a smooth tube in terms of thermal performance. 
Figures 5.12 and 5.13 compare R-124 and R-114 decreasing 
heat flux data for the 19 and 26 fpi GEWA-K tubes respectively 
for 0%, 3%, and 10% oil concentrations. The 19 fpi tube shows 
considerable improvement over the entire heat flux range while 
the 26 fpi tube shows improvement only in the low heat flux 
region. The improvement of the lower fin density tube can be 
speculated to be due in large part to better circulation and 
mixing between the fins. For both tubes at low oil 
concentrations and high heat flux, the foaming effect provides 
improvement in heat transfer similar to that seen with the 


smooth tube. 
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E. BOILING PERFORMANCE OF HIGH FLUX AND TURBO-B TUBES 
Figures 5.14 and 5.15 show the heat transfer performance 
of a HIGH FLUX tube in pure R-124, R-124/3% 011, and R-1247eeee 
O1l during increasing and decreasing heat flux respeceavemas 
Uncertainty analysis (Appendix D) was performed on the HIGH 
FLUX tube and typical uncertainty bands are also given, The 
uncertainty is largest at low wall superheat and low heat flux 
Since the uncertainty 1s of the order of the measured values. 
It can be seen that the heat transfer performance decreases 
for all heat fluxes as oil concentration increases. Bais 
decrease in performance is due most likely to the oil becoming 
trapped within the porous coating. As the refrigerant is 
evaporated, a very steep concentration gradient is developed 
within the pores; oil diffusion is also retarded because of 
the presence of the interconnected pores; this is particularly 
severe at high heat fluxes with 10% oil mixture, where a very 
large drop-off in performance is seen. In the natural 
convection region, the effect of oil is not as pronounced 
at high flux, exhibiting only minor reduction in heat transfer 
compared to the pure refrigerant. Any advantageous effect of 
foaming is negated due to the pores remaining impenetrable to 
the foam. Figure 5.16 compares R-124 and R-114 decreasing heat 
flux data for the HIGH FLUX tube for 0%,3%, and 10<caiea 
concentrations. The trend of the data for both fluids is very 


Similar. With such a high density of nucleation sites alreag, 
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Figure 5.14 Performance Comparison for Boiling R-124/0%, 3% 
& 10% O1l Mixtures for HIGH FLUX Tube (increasing flux) 
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Figure 5.15 Performance Comparison for Boiling R-124/0%,3% 
& 10% Oil Mixtures for HIGH FLUX Tube (decreasing flux) 
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present on the surface, there does not appear to be a 
Significant thermal enhancement at higher operating pressures 
as seen with the finned and smooth tubes. 

Figures 5.17 and 5.18 show the heat transfer periGiiemee 
of a TURBO-B tube in pure R-124, R-124/3% o11, and R-124/a02 
oil mixtures during increasing and decreasing heat flux 
respectively. The results of oil addition are similar to the 
HIGH FLUX tube. The notable exception is that at hiG@hiiieums 
concentrations and high heat fluxes, the heat transfer 
performance of the TURBO-B tube does not drop off as rapidly 
as the HIGH FLUX tube. Again, only a small effect of Giawiae 
the natural convection region was observed. 


Figure 5.19 compares R-124 and R-114 decreasing heat flux 


o\? 


data for the TURBO-B tube for 0%, 3%, and 10% (jaa 
GQOnGent rabionse ec this comparison 1S similar to the 
results already discussed with the HIGH FLUX tube. 
F. EFFECT OF R-124/OIL MIXTURES ON INCIPIENCE OF NUCLEATE 
BOILING 
Separate tests were conducted to analyze the incipience of 
nucleate boiling or the onset of nucleate boiling (ONB) in 
pure R-124 and how this was affected by the presence of o1l. 
These tests were taken in accordance with the general 
procedures outlined in the preceeding chapter. However, 
smaller heat flux increments were applied to the test tube 


during the natural convection region in order to obtain the 


most accurate wall superheat at the ONB. Also, during these 
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Figure 5.17 Performance Comparison for Boiling R-124/0%, 3%, 
& 10% Oil Mixtures for TURBO-B Tube (increasing flux) 
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Figure 5.18 Performance Comparison for Boiling R-124/0%, 3%, 
& 10% O11 Mixtures for TURBO-B Tube (decreasing flux) 
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incipience tests the refrigerant pool was cooled down ae 
approximately half the rate of the prior increasingmess 
decreasing heat flux runs. Each incipience test was conducted 
with a minimum of 8 hours between between re-testing to allow 
adequate time for the R-124 pool to settle and for the 
nucleation sites to de-activate. The maximum time allowed for 
the aging of the tube surface was 48 hours. Similar tests were 
repeated seven times on the smooth and HIGH FLUX tubes. Figure 
5.20 summarizes the data on a probability chart. Due to the 
statistical nature of ONB, You et al [Ref. 16] proposeaiigaae 
chart as the best means of presenting incipient data. The 
seven incipient wall superheats for each condition are given 
for the smooth and HIGH FLUX tube at 0%, 3%, and 10% oil 
concentrations as the percent probability of nucleation 
occurring. For example, the smooth tube in pure refrigerant 
never nucleated at a wall superheat less than 12°C. Also the 
data show that for a smooth tube in pure R-124, nucleation 
occurred four times at a wall superheat greater than 20 “°C, 
and three times at a lower wall superheat; this corresponds to 
a probability of nucleation occurring at a superheat of 2gmae 
between 3/7 (43)% and 4/7 (57)%. As the superheat further 
increases, the probability increases also until at a superheat 
of 28°C there is a 100% probability of nucleation. Bar-Cohen 
[Ref. 3] also reported widely varying R-113 superheats when 


‘departure from convection’ or nucleation occurred =e 
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Figure 5.20 Probability of Nucleation Occurring with Smooth 
and HIGH FLUX Tube for 0%,3%, & 10% Oil Mixture 
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probability chart of Figure 5.20 can be made more accurapome 
plotting a greater number of incipient superheats for each 
Edoe. 

For the smooth tube, there is a distinct shift CoG ioe 
Superheats with increases in oil concentration. This 
represents a delay in incipience caused by the addition of 
O1l. Memory and Marto [Ref. 10] conducted similar tests with 
R-114 and R-114/oil mixtures and reported the oil 
concentration to have no systematic effect on the point of 
incipient boiling for the smooth tube. Wanniarachcha ete 
(Ref. 9] reported an increase in the wall superheat needed to 
initiate boiling for the smooth tube in R-114 with 103 gee 
noting that 3% oil caused no increase. Memory and Marto [Ref. 
10] precicted the effect of oil on superheat needed to 
nucleate from a vapor-filled cavity by incorporating the 
physical properties of the refrigerant/oil mixture into a 
nucleation parameter (N) given by Marto and Lepere [Ref. 17]: 

N=0°T Ne 


sat/ vapor Bae 


If the parameter were to increase, a delay in the incipient 
point with increasing oil concentrations would appear likely. 
However, it is not clear how the oil effects this parameter; 
with the smooth tube in R-124/oil mixtures it resulted ina 
delay in incipience, while the HIGH FLUX tube exhibited no 


distinct delay in incipience. Also the HIGH FLUX tube has a 
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whole range of sites that can be easily activated while the 
smooth tube only tends to have very small sites which may move 
outside the range of sites that can be activated when a 
refrigerant/oil mixture is used. 

In the case of pure refrigerant, it does appear that with 
the relatively few incipient tests conducted, the wall 
superheats at inciplience for the smooth and HIGH FLUX tubes 
are larger on average for R-114 than for R-124. This initial 
result would support the approximately 50% greater nucleation 
parameter of R-114 over R-124. This subject requires further 
Steudy . 

Figure 5.21 shows an expanded view of the smooth tube’s 
heater sleeve rolled out to show the exact locations of each 
thermocouple (longitudinal and circumferential) imbedded on 
the surface of the sleeve. Several data sets were taken for 
the smooth tube and then repeated after reversing the position 
of the tube within the pool. This procedure was conducted to 
see if there were any noticeable trends in the thermocouple 
Beagamgs which may be an indication of local hot or cold 
regions within the pool. However, this reversal of the tube 
had no consistent effect on the thermocouple readings or the 
point at which nucleation occurred. This appears to underline 
the uniformity of thermocouple output within the pool. A more 


detailed survey of this data needs to be done. 
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Figure 5.21 Expanded View of Thermocouple Locations within 
the Heater Sleeve of Smooth Tube (#1 at 12 O’clock) 
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G. SUMMARY OF THE R-124 DATA 

Eases ome aaa Se2o summarize the effect of all oil 
concentrations (0,1,2,3,6,&10%) in R-124 for the five tubes at 
heat fluxes of 25 and 90 KW/m? respectively; 25 kW/m? 
represents a typical operating heat flux. For each tube, the 
ratio of the heat-transfer coefficient at a given oil 
concentration divided by that for the same tube in pure R-124 
is given. At low heat flux, the smooth tube and the 19 38 
tube give an increase in performance as oil is added, up to 
concentrations of 6%. The smooth tube does the best in this 


respect, giving an increase in heat transfer of about 23% 


oe 


(compared to a smooth tube in pure R-124) at 6 ean. 
concentration. The HIGH FLUX and TURBO-B tubes show a decrease 
in performance as oil is added, while the 26 fpi tube shows no 
change. At the highest oil concentration (10%), performance 
for all tubes is worse than that at 6%, although the smooth 
tube still indicates better heat transfer than when no oil was 
added, which is rather surprising. Similar trends were 
obtained for R-114, although at 10% oil, all tubes had poorer 
heat transfer than with no oil. At higher heat flux, general 
effects are similar, although now the 19 fpi tube shows the 
best performance of about 35% at 6% oil concentration. 
Interestingly, between 3 and 6% oil concentrations, all tubes 
except the HIGH FLUX tube show improved heat transfer over the 


Same tube in pure refrigerant. Although below 3% oil there is 


an unexplained ‘dip’ in the heat transfer coefficient ratio 


q3 


with the TURBO-B tube. The HIGH FLUX surface has very soem 
performance at 10% oil concentration, being only 35% Gree 
in the pure fluid. 

Figures 5.24, 5.25 and 5.26 compare the five tubes 
directly for o1l1 concentrations of 0, 3, and 10% respectively. 
From these figures, enhancements (ie. improvement in heat 
transfer compared to the smooth tube) can be calculategianaas 
desired heat flux. These are provided in Table 4 at heat 
fluxes 25 and 90 kW/m*. The advantages of enhanced surfaces 
are clearly apparent. The smooth tube always has the poorest 
performance (as expected) and is used as a base-line with 
which the other tubes are compared. Between the finned tubes, 
the 19 fpi tube consistently outperforms the 26 fpi tube at 
all heat fluxes and oil concentrations: enhancements range 
between 2.1-2.9 for the 19 fpi tube and 1.3-2.1 for the Zommee 
tube. Indeed, at the highest heat fluxes at 10% oil 
concentration, the 19 fpi tube is the best tube tested 
(enhancement of 2.6). The GEWA-K 19 fpi tube to smooth tube 
area ratio is 2.4, while the GEWA-K 26 fpi tube to smooth tube 
area ratio is 3.2. The HIGH FLUX and TURBO-B tubes are the 
best tubes over most conditions of heat flux and oil 
concentration, with both having fairly equal performances The 
exception is the HIGH FLUX tube at high heat flux and oil 
concentration, where performance drops off rapidly, as already 
mentioned. Over practical heat flux ranges (15-30 kW/m) 


therefore, the porous coated and modified fin tubes are the 
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bese for aioli lveencentrations. Although the TURBO-B and HIGH 
Pingx cubes Dertormance iS degraded at 10% o11 concentration, 
in-service oil separators should maintain a 3% or less 
concentration in an actual evaporator system. At this level, 
the TURBO-B and HIGH FLUX tube are superior at normal chiller 
heat flux between 10 and 30 KW/m-. 

A point of interest worth noting was that the Viton and 
neoprene O-rings used experienced considerable swelling of 
approximately 30% by volume when subjected to R-124. This was 
enough to prevent their re-use when the Teflon inserts were 
removed to change the tube. Also noted was that the O-rings 
only swelled when taken out of the system. When undisturbed, 
they maintained a good seal (this might not be true if the 
tube were sliding past the O-ring). Also significant, after 
removed from use within the apparatus for 60 days, both types 


@eeo-rings did shrink back to thelr original size. 


Table 4. ENHANCEMENT SUMMARY FOR R-124 


Heat Flux 25 kW/m? 90 kW/m? 
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Figure 5.22 Effect of Oil Concentration on Heat Transfer 
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VI. CONCLUSIONS 


For the smooth and finned tubes, significant improvements 
in heat transfer (>50%) were obtained in pure R-124 
compared with pure R-114 at the same saturation 
temperature. This was attributed to the higher saturation 
pressure of R-124 activating more nucleation sites and the 
higher latent heat of vaporization, vapor density, and 
thermal conductivity of R-124 over R-114. 


For the HIGH FLUX and TURBO-B surfaces, performance in the 
two refrigerants was Similar, since the number of active 
nucleation sites is already high. 


For pure R-124, the HIGH FLUX tube gave the largest 
enhancement of any tube over the whole range of heat flux. 


For the smooth and finned tubes, significant improvements 
in the heat transfer (>50%) were also obtained in R- 
124/oil mixtures when compared with R-114/oil mixtures at 
the same saturation temperature. This is also assumed to 
be due to the higher saturation pressures associated with 
R-124. 


Increasing the oil concentration in R-124 by up to 6% for 
the smooth and finned tubes improved heat transfer 
performance compared to the pure refrigerant. For example 
at a heat flux of at 90 kW/m’, performance for the smooth 
and finned tubes increased by 18% and 35% respectively. At 
higher oil concentrations, performance deteriorated. 


Any increase in o11 concentration in R-124 for the HIGH 
FLUX or TURBO-B tubes led to a decrease in performance 
Mehwehateat high heat Eluxes and high oil concentrations, 
the HIGH FLUX tube gave only 35% of its performance in 
pure R-124 and the 19 fpi tube gave the best overall 
performance. 


At chiller design heat fluxes (10-30 kW/m*) and oil 
concentrations of 3% maximum, the HIGH FLUX and TURBO-B 
surfaces still provide the best heat transfer which is 
about the same as with R-114. Therefore, similar or 
slightly better performance is expected with R-124 when 
used as a ‘drop-in’ replacement for R-114 ina refrigerant 
evaporator operated at the same saturation temperature. 


BagiltioOn OL Oil delays the point of incipience (occurs at 
a higher wall superheat) for the smooth tube. For the HIGH 
FLUX tube, oil addition has no measurable effect on the 
wall superheat at incipience. 
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VII. RECOMMENDATIONS 


The remaining enhanced tubes should be tested in R-124 at 
%,3%, and 10% o11 concentrations and compared to the five 
tubes above and their prior performance in R-114 and R- 
114/011 “Mixtures. 


Incipience tests performed with finned tubes and TURBO-B 
tubes should be conducted to compare the influence of oil 
on wall superheat. 


The influence of the variation in time between test runs 
on boiling nucleation (ie. onset of nucleate boiling) Gages 
be investigated more thoroughly. 


The single tube apparatus should be modified to test an 
array of two horizontal boiling tubes and compare these 
results with the R-114 tube array data previously 
recorded. 


A prototype chiller presently using R-114 should be tested 
with R-124 to monitor the performance. 
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APPENDIX A. THERMOPHYSICAL PROPERTIES OF R-124 
The following thermophysical properties of saturated R-124 
are plotted versus temperature in degrees Celsius. They were 


generated from REFPROP [Ref. 12]. 
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Saturation Pressure for R-124 
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APPENDIX B. 


Date 19 Nov 1992 
NOTE: Program name : ORPSB 
Oisk number «= OI 
Date 19 Nov 1992 
NOTE: Program name : ORPGB 
Disk number = 82 
New file name: 01119193 
TC are defective at locations 
Tube Number: ES) 
Data Set Number = 1 Bulk O11 % = 
TIME: 15:40:90] 
TC No: ] 2 5 4 
Temo ye, 3x38 S50 3.20 
Twa Tliad Tliad2 Tvapr Psat 
a eos 2.14 2.14 2.18 18,84 
Thetab Htube Qdo 
1.063 6.475E+92 6.879E+d2 
Cata Set Number =# 2 Bulk 311 % 2 
TIME: 1S:42°22 
TC No: 1 2 5 4 
Temo 3.58 3.839 3, 7e 3.6) 
Twa Tliad Tliqd2 Tvapr Psat 
353 face rear g | 2.25 10.9) 
Thetad Htube Qdo 
1.377 6.422E*02 8.84S5E+02 
Data Set Number = 3 Bulk O11 7 ® 
TIME: 15:45:28 
TO NG: 1 2 3 4 
Tempo 3.85 4.35 4.36 4.19 
Twa Tliaqd Tliad2 Tvaor Psat 
4.09 2.09 aee | 2.17 190.83 
Thetabd Htube Qdo 
1.921 6.S23E*02 1.2S3E+03 
Oata Set Number * 4 Bulk O11 2% «= 
TIME: VS 2472 
TC No: 1 2 3 4 
Temo 4.10 5.03 S.02 4.87 
Twa Tized Tliqgd2 Tvaor Psat 
4.62 2.15 2.8 i Famer ag <1) °:/80 SP: 
Thetab Htube Qdo 
2.350 6.986E+02 1.642E+93 
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ou 
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Tsumo 
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Big eis 
Tsumo 
Sa 
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Tsumpo 


-13.2 


Sece 


3.64 
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REPRESENTATIVE DATA SET 


i 8 
+94, 90050 
< 8 
=99..95 -99593 
“i 8 
= 99 Oo oe 
? 8 
29999 0099. o0 


Data Set Number = Seecurk Oil. 

Glee 15:49:03 

TC No: | 2 s 4 

Temp : 4.32 Ses 5.76 5.44 
Twa Tliad Tived2 Ivaee. Psat 
5.87 2.@S 220 2.c3 10.88 
Tretab Htube Qdp 


2.845 7.835cEtOe 2.2c8E*d3 


Cata Set Number = 6 sulk Ori <= 

Le iy Se 1§:S1:30 

TC No: | ce 3 4 

Temp : 4,42 6262 6.62 yA 
Twa Tliad Tliqd2 Tvaor Psat 
S.6) Zaz 2.16 Cee 10.95 
Trhetab Htube Qdo 


Be Scleeoetose+O2 §5,.045E+03 


Data Set Number — 7 Bulk Qil % = 
TIME: 1525470) 


TC No: ] 2 5 4 
Tempo : 4.75 Tes ee 6.75 
Twa Titad Tliqd2 Tvapr Psat 
6.4 2.15 Ns) 2.28 18.93 
Thetab Htube Qdp 
3.83! 1.@31E+03 3.9SI1E+@3 
Data Set Number = 8 8ulk O11 % = 
TIME: 16:56:43 
TENG? | 2 3 4 
Temp : S.14 8.28 8.15 7.38 
Twa Tliee Tliqd2 Tvapr Psat 
6.66 2.c4 ener 2.35 11.02 
Tnetab HKtube Qdo 


4.302 1.2S36€*03 S$.410E+03 


Oata Set Number = 9 osulk Onl =x. = 

TIME: 15-59-50 

TC No: | 2 3 a 

Tempo : 4,78 6.89 Sli 6.85 
Twa Tired Tliqgd2 Tvapr Psat 
5.97 Zoa4 2uc4 Zes2 0.97 
Thetab Htube Qdp 


SeocO 2.tSl6+95 7.649E+03 


Data Set Number = 10 Bulk Oil % = 
TIME: re: 02:35 
TE No: | 2 3 4 
Temp : 4,83 S.18 ei 4.84 
Twa Tliad Tliqgd2 Tvapr Psat 
4.86 2.50 C65 face A0.35 
Thetadb Htube Qdo 


2-565 842 )150E+@3S 1.072E+04 


Data Set Number ~ 11] Bulk O11 % » 

TIME: 16:17:38 

TC No: | 2 5 4 

Temp : S.14 5.66 s.s93 Sone 
Twa Tliad Tliagd2 Tvaor Psat 
5.29 ciao Eevee Prior: 10.82 
Thetab Htube Qdo 


5.044 5§.0986€+03 1,.S526+04 


oa 
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5 IS) 
4,58 4.56 
Tsump 
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S.12 / aa ae 
Tsumo 
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35.0 


S20 S.22 
Tsumo 
=(2,7 
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$.SS 35S) 
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4.88 4,46 
Tsump 
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Tsump 
“12.4 
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=99..99 


=39 359 


=J9).39 


=39:; 99 


=a, 99 


= 95799 
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Sshs |i tshs 


= sls Its Be 


=99.959 


hs eee be 


=o 9299 


Data Set Nuesnder = 12 Bulk O11 2 = 
TIME: 16:20:31 
fC Noe | Z 3 + 


Temp 5273 6.34 6.28 S.928 
Tw Tliaa Tliad2é Tvege “san 
5.83 2226 2.19 2120 18235 
Tretad Htube Qdo 
2.630 6.328€+83 2.297E+04 

Oata Set fhuscer = 13 “SGlk Ope 

PERE: WG223°15 

TC No: | 2 5 4 

Temp : 6.61 Ta18 Gales 5.64 
Twa Tliad Tliae2 Tvase Psat 
6.62 Zon 2uce 2.28 10.95 
Thetab Htube Qdp 


4.316 7.794€+03 3.364E8+04 


Data Set Number = 14 S8ulk O11 2% = 
TIME: 16:24:49 

FC NG: | 2 3 4 
Temp : 7.42 7.84 7285 Ee i> 


Twa Tliad Tliqad2 Tvapr Psat 
areas 2.35 PHP) Zune CIOeS7 
Thetab Htube Qde 
4.895 9.726€+83 4.761E+04 


Data Set Number = 15 8ulk O11 % » 
TIME: 16:29:15 
Te awe: | 2 3 4 
Temp : 8.17 Be o2 Sa Soi 
Twa Tliad Tliqd2 Tvaor Psat 
7.78 2:29 ence 2.28 10.93 
Thetabd tube Qdoe 
S.S50!1 1.300E+84 7.1SIE+O4 
Data Set Number = 16 8ulk O11 3% # 
TIME: 16:33:55 
TC No: | 2 3 a 
Temp : 8.79 8.96 9.33 8.87 
Twa Tliad Tliqd2 Tvapr Psat 
<i a 6 2.28 2. ie 2.28 18.86 
Thetab Htube Qdo 


6.169 ' 1.5556+04 9 S77e<04 


Data Set Number = 17 8Sulk O11 2 # 

TIME: 16:36:37 

Leone: \ 2 5 4 

Tesp : 7.44 7.44 Tease 7.89 
Twa Tivad Tliqd2 Tvapr Psat 
lace 2s rare li 2.17 18.83 
Thetad Htube Qdoe 


$.070 I1.l238E+84 6.235E+04 


Oata Set Number = 18 Bulk O11 2% = 

TIME: 16:48:00 

TG Ne: ] 2: 3 4 

Temp : 6.63 7.03 rae Bee 
Twa Tliqd Tlrad2 Tvaor Psat 
6.57 o.ce 2.18 2.24 18.98 
Thetab Htube Qdo 


a.Se 9.6S7E+03 4,180E+@4 


a2 


S30 


S 
6.21 

Tsumo 

-11.9 


oee 


6.85 
Tsump 
-11.8 


3.0 


Tsumo 
-11,8 


8.32 
Tsumo 
-11.7 


9.12 
Tsumo 
-11,4 


6.98 
Tsump 
-11,4 


ot 


5.24 =Jona0 


6 " 
5. 997-33 
6 7 
B.49 35292 
6 qi 
7209 294809 
6 7 
7.99 =J39-393 
6 7 
6.56 <=J9.93 
6 7 
5). 36 s- SSR 


=33.43 


=937 43 


a) =)- =) | 


-99 35 


= 99a 


“339304 


“39-33 


Data Set Number © 19 8ulk O11 % » 

TIME: 16:43:34 

TO Ne: s ce 4 

Tees : 6.99 5,65 6.55 Sy o 
Twa i lied Titqde Tevage Fs4t 
6.'4 emia 2.20 2.24 !28.92 
Thetad Htubde Qdp 


3.899 7.859E+@3 3.067E+04 


Cata Set Number = 2@ Bulk O11 % = 

T Lae: 16:48:99 

TC No: 1 Z z) 4 

Temp : 5.62 6.13 6.13 $.32 
Twa lived Tliqd2 Tvapr Psat 
5.78 ieee) 2,14 2.19 128.85 
Thetab Htubde Cdo 


3.506 6e754E+O3 2. 368E+04 


Data Set Number = 21 Gulk Qil 2% = 
TIME: 16:48:34 
TC No: | 2 3 4 
Tempo : S.17 S.73 S.71 S.36 
Twa Tliad Tliqd2 Tvapr Psat 
5.38 265) once eves. 10.91 
Thetad Htube Qdop 


3,045 -S.55SE+03 I692520¢6 


Data Set Number = 22 Bulk O11 2% = 
TLME: ie Slood 
Geno. 1 2 3 4 
Temp : 4.69 5.29 S.1S 4.385 
Twa Tliad Tliqd2 Tvapr Psat 
4.82 Cael cole 22 VO. a7 
Thetabd HAHtube Qdop 
2.613 4.9256+03 1.1826 +04 
Data Set Number #= 23 Bulk Oil % = 
1 On) ar 16:54:33 
TC No: 1 2 3 4 
Temp : 4.32 4.77 4.63 4.45 
Twa Tliad Tliad2 Tvapr Psat 
8,85 2.30 eel CARTS el B's es 4 
Thetabd Htube Qdp 


27 e6es Gb55b 70557. 985E +05 


Data Set Number = ZS Gulk Gil % @ 
TIME: 16:97:59 
TC No: | 2 3 * 
Temp : 2,94 A055) 423 4.12 
Twa Tliad Tliqgd2 Tvapr Psat 
4.38 2.24 2.'5 el lie oz 
Tretsb  Htube Qde 


3363 (2.8396+03 ~5.289E+03 


Data Set Number = 25 Bulk O11 2% = 
TIME: 17:01:02 
Ge fear 1 Zz ) 4 
Temp : 3.56 5.9) 3.84 3.75 
Twa Tliaqad Titee2 tIvaor Psat 
Sere Zee 2.988 2.14 18.80 
Thetab Htube Qdp 


Sc. Nese «05 a 4192703 


23 


6.4) Ses. - 35.99 


$.92 S.19 -99.99 
Tsumo 


aililicy: S 


3.2 


S.42 4.88 -99.99 
Tsumo 
=Aliding 7 


3.0 


%.92 4.44 -99.99 
Tsumo 
-11.8 


4.55 Se 39). 39 


4-all 3.91 -99.99 
Tsumo 


=l2.0 


3.2 


3.74 3.68 -99.99 
Tsumo 


“le. 


=e\e)os| 5) 


=39.993 


mile |Bgo| 


=99.99 


ole bo Ane ee 


gan 99 


-a239 


Oate Set Number - £6 Suib Jil 3 « 
TIME: 17:@35:18 
TC No: | 2 3 4 
Teao : 3.S6 3.83 Bi AT et! 3.76 
Twa Tliad Thiqdc Tvayur Psat 
3.69 Ree Pas) ceer Tecde 
Thetad Htude Qdo 
1.420 1.8809&€*@3 2.59786+03 
Date Set Number © 27 gult O11 2 = 
TIME: 17:06:85 
TC No: | 4 3 4 
Teap : Sia 3.63 3.56 3.82 
Twe Tliad Tliad2 Tvaer Pset 
3.S@ 2.35 2.28 2.24 180.89 
Theteab Htube Qdo 
1.260 1.S71E*@3 1.981E+03 
Date Set Number © 28 Sulkh O11 2 © 
TIME: 17:08:64 
TC No: | Z 3 4 
Taap : neg | 3.42 3,34 3.34 
Twa Tliaqd Tliqgd2 Tveor Psat 
332 2.26 2.18 2.14 186.98 
Thetad Htube Qdo 
1.185 1.1S9€*@3 1.373&+83 
Date Set Number = 29 Sulkh O11 2 © 
TIME: 17:10:55 
TC No: | 2 3 * 
Teap : S25 re 3.18 3.16 
Twe Tliaqd Tliqgd2 Tvaor Pset 
Sole 2.28 2.25 2.23 10.89 
Thetedb Htubde Qdo 
.937 1.@S6E+@3 9.893E+@2 
Date Set Number = 36 Bulk O11 2 © 
TIME: li7s03 336 
TC No:- 1 re 3 4 
Taap : 3.20 3.40 3.36 Sn5 
Twe Tliad Tliqd2 Tvaeor Pset 
3.30 Pred 2.68 2.10 16.77 
Thetebd Htube Qe 
1.205 1.2466 °@3 1.S@1E+@3 
Oete Set Number © 31 Bulk O11 2 © 
TIME: 17:16:@@ 
TC No: | 2 3 4 
Temo : 3.18 3.28 3.260 3.16 
Twe 1a ey-r-| Tliqgd2 Twapr Poet 
3.18 2.30 avew 2.28 10.93 
Theta’ Htube Qdo 
.9@7 9.9228 °O2 9. 001E +62 
Oate Set Number © 32 Bulk O11 2 © 
TIME: 72 is3s 
TC No: 1 z 3 4 
Teap : or3< 3.08 2.92 2.89 
Twe Tliaqa Tiiqgd2 Tvapr Psat 
2.94 Zo Ge We 2.18 18.81 
Theteabd Htube Qdo 
.796 8.888E*@2 7.0746+82 
Oate Set Number = 33 Sult O11 2% © 
TIME: 17:19:16 
TC No: | Pe 3 4 
Teno : 2-32 2.98 2.88 2.86 
Twa Tliad Tliaag2 Tveor Psat 
esi Cole 2.19 2.16 128.82 
Theteb Htubde Odo 
-753 7.746SE*82 §.829E+82 


94 


S 6 zu 
3.69 3.61 -99.99 
Tsao 


~ * a 
ae Une oC5 


-99.39 
Tsumo 
=\2o5 


-939.99 
Tsuse 
-12,4 


3.8 


Sl2 ere | 
Tsump 


wiles 


=99740 


3.60 


5.25 3.28 -99.39 
Tsuaop 


=\<76 


3.@ 


Brae) a62o mogeae 
T $u"®O 


“12.7 


3.8 


2.91 3.03 -99.99 
Tsumo 


-12.7 


3.8 


$ 6 7 
«.8? 7g -39.59 
Tsumo 

“12.8 


8 
239-39 


=39...99 


ces be BB 


=3'9'. 219 


madi 


= clelnee) 


eee ele 


#39:.93 


APPENDIX C. SAMPLE CALCULATIONS 


Data file number D1221SM0 was used to conduct the sample 
Calculations. The saturation temperature was 2.28 °C with a 
eeeete smooth tube heat flux of 31290 W/m? . 

A. TEST-TUBE DIMENSIONS 
we 015875 m 

ee—--Ol27 m 

D,=.01245 m 

f——2032 m 


Lu=.0762 m 


B. MEASURED PARAMETERS 
V=141.95 volts 

T=2.34 amps 

iJ=16.94 °C 

m2-17.15 °C 

T3=16.42 °C 

M4=16.22 °C 

T5=15.66 °C 

T6=15.48 °C 


jee 2.28 °C 


sat 


k.=344 W/mK 


25 


C. OUTER WALL TEMPERATURE OF THE BOILING TUBE 
p=7°D,=7: .015875 m =.04987 m 


Ac=m (D,?-D,7) /4=m((.015875) 7- (.0127)*) /4 


Ae=gelaxlo 2 me 


Q,=VI=141.95:2.34= 332.16 W 


6 
Tayg=bTn/n=), ‘Tn/6=16.3 Se 
, 
Dor lava en ley ae 2 tee 
: . 
Tyorl6.3-1 332,16"(1n(.015875/.01245))/ (2s 
a oes 


G27 Tee G2 12=0. 78-1494. CC 


wo sat 


D. PROPERTIES OF R-124 AT FILM TEMPERATURE 
The thermophysical properties of R-124 were obtained from 
REFPROP [Ref. 12] and are shown in Appendix A. 


Tea Dot lene) 2a 6242) 2 ee eee 


p=3.4593-(.0426)T-+(3.9495x10*)T,°- (4e193x0 epee 


(20 709x108 


26 


u=3.4593- (.0426) -9.2+(3.9485x10-*) - (9.2) %- (4.193x107®) + (9.2) 3+ 
ro70ox10 °) + (9.2)? 


u=3.098x10°* N-s/m? 


memea4 8- (2.8619) T.- (6.7267x10 7) Te*- (7.2852x10°>)T,? 
emit 4 e- (2.8619) -9 .2- (6.7267x10 7) - (9.2) 2- (7.2852x10°7) - (9.2)2 


0=1407.84 kg/m? 
v=u/p=3.098x10°*/1407.84= Oe alse Oia m*/s 
fee do x10 “-(3.5436x10 “)T,- (1.9545x10° ’)T,*+(3.1835x10°7)T,? 
Pees 9 110 ~- (3.5436x10 *)* (9.2) - (1.9545x10°’)- (9.2)7+ 
i s35e10°7)-(9.2)7 


k=7.192x10°* W/mK 


pee 0542+ (2.1405x10 ~)T,+(1.0709x10 ”) T,*- (6.4721x10°*) T,?- 


4324x100. +(4.136x10°11) T° 


cee 5427 (2. 1405x010) -) (9.2) +(1.0709x10°>) « (9.2)*- (6.4721 


atom eee eo 4 ) 9D) +4 (45136x10 77) - (9. 2)° 
eo 4.74 wy eg K 


a=k/(p+c,) = Tegner) (407 84-1074 .74) 


a=4.753x10°? m@/s 


oe 


B=- (Ap/AT) /p=2.094x10°7 (1/K) 


wT 
Kk 
i! 


v/a=4.63 


E. HEAT-FLUX CALCULATION 

The average natural convection heat-transfer coefficient 
at the non-boiling ends of the test tube is calculated time 
the Churchill and Chu (Ref. 4] correlation 
For simplification let R represent the following in the 


‘elope et wale nal 


R=[(g'@+D,?: 6+ tanh (m+ Lu) ) /(v-a-Lu-m) ]?/8 


then 


h=k/D,[.6+.387(R/(1+(0.559/Pr) 9/18) 8/27] } 2 

where 

m=[(h-p)/(k,'Ac)]1/? 

and an iterative process is used in program DRPGB to compute 
the value of h by starting with the value for h Gime 


(W/m?:K). The resulting values for h and m were: 


h=265.43 (W/m*:K) 
m=23.23 (1/m) 
therefore: 


Qe=(h'p:k.:Ac)?/?-6-tanh(m ‘Lu)=7.44 W 


ee) 


F. HEAT FLUX THROUGH ACTIVE BOILING LENGTH 
Q=0,,-2:Q-=332.16- (2:7.44)=317.28 W 
Ab=m'D,:U=7:0.015875:0.2032=1.0134x10°* m* 
q=Q/Ab=317.28/1.0134x10°7=31,308 W/m? 
h=q/6=31,308/13.84=2262.1 W/m*-K 

The following resultS were produced by the data 
acquisition and reduction program DRPGB: 
g=31290 W/m? 


6=13.84 °C 


h=2262 W/m*:°K 


TABLE 5. DIMENSIONS OF BOILING TUBES TESTED 






TUBE Dino ZaeiiieiDae Mmm DOM mm moe mm» bu mm ee 


W/moK 


203.2 344 
GEWA-K a ON: ee), om ea 1A a 20 Se alae © 344 
noe tpl 
GEWA- K ; 2 Ose aero 
Eom LpD1 ai | . | es 7 


Brees Olli1. 6 4 o.3.8 ee ease eg |) 208.2) 7602 1 3908 


HIGH IL ae, ISVs, eZ heros, ZS 22 .\ 7622 30 ate 
FLUX 


* HIGH FLUX tube porous coating is on a copper/nickel base 










oe) 


APPENDIX D. UNCERTAINTY ANALYSIS 

Using methods of uncertainty analysis outlined by Kline 
and McClintlock [Ref. 18] and performed by Sugiyama [Ref. 2], 
four data points are analyzed. The data points chosen were for 
a Smooth tube and HIGH FLUX tube. For each tube, there was a 
data point at a high heat flux and one at a low heat flux. The 
following 1S a sample uncertainty analysis calculation from 
data file number D1220SM0 with a Smooth tube, 0% Gi eee 
heat flux of 86,000 W/m*. The results of the remaining sample 
points are listed in Table 6. Uncertainties are given as a 
percentage of the calculated parameter, i.e. a relative 


UnGerReEatany . 


A. UNCERTAINTY IN HEAT-TRANSFER RATE 


ea volts Ver 2g32> VOLES 


Estimated uncertainties in the measured quantities I and V 
6I=+ 0.025 amps 


O6V=+ 0.05 volts 


fai, 7162 12-3. 32 amps 
V=25 Ve=2337 3 eVolee Qn=VI 


where: = uncertainty in measurement 


100 


60,/Qn= ( (6V/V,) 24 (@s0/ Te) a) AZ 
6Q,/Qn= ame on 305) -4(.025/199)-)-/2 


ew) O,= 1.37 & 
B. UNCERTAINTY IN SURFACE AREA 


Ab= nm:Do'L 
Be= 15.88 (mm) 6D_= .1 (mm) 
m= 203.2 (mm) 6L= .1 (mm) 
6Ab/Ab= ((6D,/D,)* + (6L/L)2)1/2 

Oem Ab= ({.1/15.88)2 + (.1/203.2)°*)1/2 


6AbD/Ab= .63 % 


C. UNCERTAINTY IN WALL SUPERHEAT 


ao TT. -T = ee. jie =O. oC 


wo T 


sat sat sat 


en) tavg Sn liniD,/D,)/ (2+: Lk.) ] 
Tn= thermocouple readings 

eee Hr /n 

fee 21.63 °C 

iiea-21.41 °C 

m-20.63 °C 

T4=20.94 °C 

mee20.06 °C 

T6=20.40 °C 


Tayg=20.85 °C 


EEA 


standard deviations ((£ (in 


Standard deviarweon—= 0255.6 


The logarithmic term in the equation for Two can be 
neglected for the uncertainty analysis as it 1S very Sma 
comparison to the standard deviation. 

Wien ata= 05S eae 

6T o=standard deviation=55"-¢ 

Ok sia) Ae 

OAT /AT=((6T 2 /AT) “a omeee, Ale ee 

SAT /AT=( (. 55/7 18,4621) 2 Gs Ody beae ee 


OO AR=2295. 5 


D. UNCERTAINTY IN HEAT FLUX 

q= (0,-2:'Q¢) /Ab 

Q,=VI Q),=890.56 W 6Q,=+12.2 W 

assuming the same proportion in the uncertainty fons 
q=86000 W/m? Bo= (0 sO mos game 

Or 270s S7isei ow 


Qp= 9.38 W 60¢=0.128 W 


6q/ = ( (6Qn/ (Qn-20¢) ) 2+ (25Q¢/ (Qy-2Q¢) ) 7+ (GAD/Ab) 7] 7”? 
6q/q=[(12.2/(890.56-2:9.38))2+(2+(.128)/(890.56-2:9.38))*+ 
(6239%10 2/ 0moten) 24 


OG /.G= 1 SA 


no 2 


E. UNCERTAINTY IN BOILING HEAT-TRANSFER COEFFICIENT 
=o / AT 

Baf/h= [(6q/q)2+ (6AT/AT) 2} 4/2 

Oh/h= [(.0154)2+(.0295) 2]?/2 


pay n= 3.3 % 


TABLE 6. UNCERTAINTY ANALYSIS OF FOUR DATA POINTS 


Parameters Smooth tube |} Smooth tube | High Flux |e Salto ln Met De BI e 


4019 W/m? 86000 W/m? 2244 W/m? 88860 W/m- 


enjh (%) 





All uncertainty analysis calculations were conducted 
Guring decreasing heat flux runs. The wall superheat 
uncertainty (6AT/AT) for the HIGH FLUX tube was significantly 
greater than that of the smooth tube due to a greater 


deviation in the thermocouple temperature readings. 


iO) 3 


APPENDIX E. SETUP PROGRAM 
The following setup program was used in the preparation of 
the apparatus prior to commencing testing. The program was 
written in Hewlett-Packard Basic 5.0 for both the Hewlett- 
Packard 9300 and 9852A series data acquisition/ controls 


The setup program performs the following: 


1. Monitor the Sump temperature 
2. Monitor the evaporator liquid temperature 
3. Measurement and readout of all thermocouple channels. 


4. Measurement and readout of the power supplied to the 
tube cartridge heater. 


5. Measurement and readout of the power supplied to the 
auxiliary heaters if used. 
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PROGRAM: SETUP 

BATE: AUGUST 3,139! 

PROGRAMMER: LT DEAN SUGIYAMA 

ROGUE VECSS  SEANNIEREAKE JAN 2s, 1992 

Pen / Coy sea?) 

CoM OsleOubedl,cotc? .94569.-767345.929S,°8005595.81 -3247486589 §.97E39E+ 


eo lsce tls 


OATA 3.94078E+14 

REAO C(#) 

ON KEY 1,15 60TO 27 

PRINTER IS 1} 

PRINT 

PRINT 

PRINT USING “4X ,°°SELECT OPTION*S** 

PRINT USING “6X ,**@=MONITOR SUMP**° 
PRINT USING °6X,°*1=MONITOR LIQUIO*** 
PRINT USING “6X ,°°2=CHECK THERMOCOUPLES” *~ 
PRINT USING “6X ,*°°3=CHECK MAIN HEATER* °° 
PRINT USING °6X ,*°°4=*CHECK AUX HEATERS* °° 
PRINT USING “6X ,°"S=EXIT PROGRAM" °° 
PRINT USING °4xX ,°°NOTE: KEY 1 = ESCAPE**” 
BEEP 

INPUT Ido 

IF Ido-S THEN Ido2#5 

IF Ido=@ THEN SQ 

IF Ido=1 TREN 155 

IF Ido#2 TREN 173 

IF Ido#3 THEN 195 

IF Idoe4 THEN 195 

IF Ido#S THEN 231 

PRINT 


PRINT 

PRINT “SUMP TEMPERATURE CES C 
PRINT 

OUTPUT 789;°AR AFI! ALI} VRS° 
OUTPUT 789;°AS SA* 

Sum2Q 

FOR J=1 TO 5 

ENTER 709;E 

Sum=Sum+E 

NEXT J 

Eave=Sum/S 

Temp =FNTvsviEave) 

PRINT USING “4X MOO.00";Temp 
BEEP 

PRINT 

WAIT S 


ILC 


rs! GOTO S2@ 

sat 

bes PRINT 

1S6 PRINT “LIQUIO TEMPERATURE DES C° 
'S$3 PRINT 

SEI OUTPUT 739;°AR AF@B ALO9 VRS" 
162 Sum= 

161 FOR 121 TO 2 

eZ OUTPUT 7aSe AS SA 

163 ENTER 729;€ 

164 Sum=SumeE 

155 NEXT | 

166 EavesSum/2 

167 Temp 2FNTvsviEave? 

168 PRINT USING °44 MCO.00°; Temp 
169 BEEP 


178 «WAIT S 

171 60TO 15S 

1721 

173. PRINT 

176 PRINT “CHANNEL TEMPERATURE OE6 C* 


eet OUTPUT 709; °AR AFQQ@ ALI! VRS° 
178 FOR [#! TO 12 

179 OUTPUT 7@09;°AS SA° 

182 Sum=Q 

181 FOR J#!1 TO § 

182 ENTER 7@9;E 

183 Sum=SumetE 

184 NEXT J 

(8s EavesSua/S 

186 Temp 2FNTvsviEave) 

187 PRINT TAB(3);1sTABC IS): Temp 
188 NEXT I 


189 BEEP 

192 WAIT S 
191 60TO 173 
194) 

a5 PRINT 


196 OUTPUT 709;°AR AFLA ALZ22 VRS" 
197 FOR I#!1 TO 3 

198 OUTPUT 7834 AS SA: 

199 Sum=Q 

222 FOR J=1 TO S 

2) ENTER 7@3;E 


cQ2 SumsSumt& 


106 


203 
204 
20S 
289 
210 
Zan 

zalke 
216 
2a7 
2K 


NI ro 
ri - 
Y wo 


mrs —- OV woos Oo wo & Uf — 


FIFI TICTIN TIN MINN PO Po fo 
W tH OF WW fa fi roars N FIM PO PD 


Os 
~ 


NEXT J 

IF T=] TREN Vo!lt=Sum/§ 

IF T22 ANDO Idoe3 THEN 

PRINT “MAKE SURE VOLTAGE BOX IS SET TO MAIN HEATERS" 
Amo=Sum/S 

END IF 

IF {23 AND Ido=4 THEN 

PRINT “MAKE SURE VOLTAGE BOx IS SET TO AUX HEATERS’ 
AmozSum/S 

alia) ie 

NE«T I 

AmpzABS( Amo?! .9182) 

YVolt=ABS(VoltecS) 

Power=VolteAmp 

ResistencesVolt/Amp 

PRINT 

BEEP 

PRINT “VOLTAGE(V) CURRENT(A) RESISTENCE( onms) POQWER(W)* 
PRINT 

PRINT USING °1X,S(MO000.00 4X)" Volt Amp Resistence Power 
WAIT § 

SeTO' 135 

BEEP 

PRINT 

PRINT “THAT'S ALL FOLKS!" 

END 

DEF FNIvsviv) 

Com (Ce GC 7? 

T20(Q) 

FOR I#t TO 7 

T2zT+C(T evur*t 

NEXT | 

TeT+S&. S26897E-cr Tet 3. °81199E-2-TeS  O689CS9E-5) 
RETURN T 

FNEND 


LOW 


APPENDIX F. PROGRAM DRPGB 
The data acquisition and reduction program DRPGB is 
written in Hewlett-Packard Basic 5.0 for the HP 9300 series 


computer and listed on the following pages. 
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ni 
Q 


182 
192 


FILE NAME: ORP6B 

OATE: Octobar 19, 1964 

REVISED: July 11, 1991 BY Oaan Sugiyama 

REVISED: Aug 16, 1992 BY Gaorge Bertsch (including R-124 propertias ) 
COM /Idp/ Idp 

PRINTER IS ! 

CALL Select 

INPUT “WANT TO SELECT ANOTHER OPTION (1#Y ,Q8N)?° ‘sel 

IF Iselel THEN GOTO 72 

BEEP 

BEEP 

PRINTER IS 1} 

PRINT "OATA COLLECTION/REPROCESSING COMPLETEO™ 

ENO 

SUB Main 

Com /Ipeop7 Ift 

COM /Idp/ Idp 

COM /Cev1Ct2) leal 

COna/Wil/ 02-071 .0e Llu .Keu 

OtmeEmt(12).1¢iz2) Olaci3) 02a 13) O14a(13),00a(13) .Lati2) Luatt3) Keuati3), 


Etc i9),Tm$(4)01S) 


200 
212 
222 
230 |! 
242 
252 
262 
272 
288 
292 
320 
312 
322 
332 
342 
352 | 
35) 
362 
372 
382 
392 
422 
412 
422 
432 
43} 
432 
433 
442 
452 
462 
472 
482 
492 
5200 
512 
$22 
$32 


540 | 
552 | 


OATA @.1@086091 ,25727.94369 ,-76734S.8295 ,78025595.861 
OATA -92474B86569 ,6.976BBE+!! ,-2.66192E+13 .3.9407BE+14 
REAO Cie) 

OATA “Smooth*® ,*High Flux” ,"Tnermoaxel-E” ,"Thermoexel -HE” 
OATA Smooth,High Flux ,Turbo-8B,High Flux Mod,Turpo-B Mod 
REAO Tn$(e) 

PRINTER IS 701 

BEEP 

IF Idp#4 THEN PRINTER IS | 

IF Idpe4 THEN GOTO 2662 

INPUT "ENTER MONTH, OATE ANC TIME (MM:00:HH:MM:SS)” Oates 
OUTPUT 709:3°TO";Oatas 

OUTPUT 709;°TD" 

ENTER 709;:0ata$ 


PRINT 

PRINT " Month, date and time :";Oate$ 
PRINT ° Oata :° ,OATES( TIMEOATE ) 

PRINT 

PRINT USING "1X ,""NOTE: Program name : ORPGB*"" 
BEEP 


INPUT "ENTER OISK NUMBER” On 

PRINT USING “16X,"°°Oisk numbar = °° 227°30n 

BEEP 

INPUT “ENTER INPUT MOOE (02#3054A ,)eFILE)* Im 

BEEP 

Tft=2 

INPUT “SELECT FLUID (0#114,10124)° lft 

BEEP 

INPUT “SELECT HEATING MOOE (Q*ELEC; !eWATER)° Ihm 
BEEP 

INPUT “ENTER THERMOCOUPLE TYPE (Q@NEW,1"#0L0)"* ,Ical 
IF ImeQ THEN 

BEEP 

INPUT “GIVE A NAME FOR THE RAW OATA FILE” ,02_filas 
PRINT USING “16X,""Naw file nama: °° ,14A";02_filasé 
Sizale22 

CREATE BOAT O2_file$ Sizel 

ASSIEN @File2 TO O2_filas 


OUMMY FILE UNTIL Nrun KNOWN 


ION 


Soe Ol_filaS="Oummy: 
572 CREATE BOAT O1_filaS ,Sizel 
588 ASSIGN @Filal TO Ol_filas 
590 OUTPUT @Filal;Oates 
62 IF Imhm=@ THEN 
612 BEEP 
622 INPUT “ENTER NUMBER OF OEFECTIVE TCS (Q@*QEFAULT, 4 MAX.)° Idte 
63 IF Idtc*#@ THEN 
642 Ldtct=@ 
652 Ldtc2s@ 
651 Ldtc3<@ 
652 Ldtc4<@ 
668 PRINT USING "16X,°"No defective TCs axist*** 
672 ENO IF 
682 IF Idtce!l THEN 
690 BEEP 
700 INPUT “ENTER DEFECTIVE TC LOCATION (1-8)° ,Latc! 
710 PRINT USING “16X,°°TC 18s defective at location °° OO*sLdtc! 
722 Ldtc2=@ 
re Ldtc3#2 
122 Latc4e@ 
730 ENO IF 
740 IF Idtc#2 THEN 
750 BEEP 
76 INPUT “ENTER OEFECTIVE TC LOCATIONS (1-8)° Ldtcl ,Ldtec2 
772 PRINT USING “16X,°*TC are defactive at locations ** 00,4X OO°;:Late! Ldte2 
Pare Ldtc3<Q 
Tz Ldtc4=Q 
782 ENO IF 
78) IF Idtes3 THEN 
782 BEEP 
783 INPUT "ENTER DEFECTIVE TC LOCATION (1-8)" ,Ldtel ,Ldte2 .Ldte3 
784 PRINT USING "1GX,°"TC 18 cefactive at location °° ,3(A4X 00) iLdtel Ldte2d Ua 
te3 
78S Ldtc4=#@ 
786 ENO IF 
Tae IF Idtce=4 THEN 
788 BEEP 
789 INPUT "ENTER OEPECTIVE TC LOCATION (1-8)° Lote! ,Ldte2 Ldte3 Ldte4 
792 PRINT USING “16X,°°TC is dafective at location °° ,4(4xX O0O)"1tLatel ,Ldtco tea 
tes. Ldtes 
792 ENO IF 
793 IF Idte>4 THEN 
880 BEEP 
B12 PRINTER IS 1 
B22 BEEP 
830 PRINT “INVALIO ENTRY" 
B42 PRINTER IS 781 
85 GOTO 610 
BBQ ENO IF 
870 ENO IF 
882 OUTPUT OFilalibLdtc! ,Ldte2 Ldtc3 ,Latcé 
B90 ' Im=1 option 
322 ELSE 
912 BEEP 
22 INPUT “GIVE THE NAME OF THE EXISTING OATA FILE* O2_files 
932 PRINT USING “16xX,°°Old file mame: °° 14A°:O2_file$ 
942 ASSIGN @Fila2 TO O2_files 
950 ENTER @Fila2iNrun 
362 ENTER @File2:Oolds 
372 PRINT USING "!16X,""This cata sat taxen on: “" 14A" Dold$ 
380 ENTER OFila2iLdtc! ,Ldtc2,Ldtc3 ,.atesd 
392 IF Ldtc1>@ OR Ldtc2>@ OR Ldtec3>@ OR Ldtc4>0 THEN 
1200 PRINT USING “16X,""Thermocouples were cefactiva at locations: ** ,4(30,4x)"5 
Létel bete2 Edtes-bates 
'@12@ ENO IF 


agak le 


1022 
1032 
1042 
150 
1268 
1061 
1062 
1Q07Q! 
1QB0 
1292 
1180 
1112 
1120 
1130 
114 
115 
116 
117 
118 
1190 
1200 
1210 
1220 
1230 
1240 
1252 
1260 
1270 
128 
129@ 
1300 
131 
1320 
133 
1342 
1352 
1362 
1372 
1380 
1390 
1480 
1412 
1420 
143 
1440 
1450 
1460 
1472 
1480 
1490 
1580 
1512 
1520 
1530 
1542 
155 
1560 
1570 
1582 
1539 
1620 
1612 
1620 
1630 


1642 
165 


ENTER @Fiiecilet 

ENO IF 

Idtc#@ 

IF Latct>@ THEN Idtceldtc+! 

IF Ldtc2>0@ THEN Idtceldtc+! 

IF Ldte3>@ THEN Idtceldtc+! 

IF Ldtc4>@ THEN Idtceldtc+! 

IF Ime#Q@ AND Ihme! THEN 1595 

BEEP 

INPUT “WANT TO CREATE A PLOT FILE? (08N,1#Y)" [plot 
IF Iplot#! TREN 

BEEP 

INPUT “GIVE NAME FOR PLOT FILE" P_files 
CREATE BOAT P_fileS,4 

ASSIGN @Plot TO P_file$s 

ENO IF 

IF Inmel THEN 

BEEP 

INPUT “WANT TO CREATE Uo FILE? (@=eN,1¢Y)° ,luf 
IF Iuf@#! THEN 

BEEP 

INPUT "ENTER Uo FILE NAME® Ufiles 

CREATE BOAT Ufile$,4 

ASSIGN ®@Ufile TO Ufiles 

ENO IF 

BEEP 

INPUT “WANT TO CREATE Re FILE? (Q"N ley) Ire 
IF Ire#! THEN 

BEEP 

INPUT “ENTER Re FILE NAME® Refiles 

CREATE BOAT Refile$,10 

ASSIGN @Refile TO Refiles 

END IF 

ENO IF 

PRINTER IS 1 

IF Im#@ THEN 

Been 

PRINT USING “4X ,°°Select tube numder””*” 

IF Inm=@ THEN 

PRINT USING “6X ,°°Q@ Smooth 4 inch Ref°*" 
PRINT USING “6X ,.°°!1 Smooth 4 inch Cu (Press/Slide)*** 
PRINT USING °6X,°°2 Soft Solder 4 inch Cu"*° 
PRINT USING °6X,°°3 Soft Solder 4 inch HI6H FLUX*°* 
PRINT USING °6X,°°4 Wieland Hard B inch**°” 
PRINT USING °6X,°°S HIGH FLUX B ineh*°* 
PRINT USING °6X,°°6 GEWA-K 428 Fins/in*** 
PRINT USING °6X,°°?7 GEWA-K 26 Fins/in’*" 
PRINT USING °6X,°°8B GEWA-T 19 Fins/in’®°” 
PRINT USING °6X,°°39 6EWA~T OR GEWA~-TY 26 Fins/in‘’*” 
PRINT USING °6X,°°1@ THERMOEXCEL~E°°* 

PRINT USING °6X,°°1!1 THERMOEXCEL-HE*** 

PRINT USING °6X,°°12 TURBO-8°** 

PRINT USING °6X,°°13 GEWA-K 19 Fins/in®’* 
ELSE 

PRINT USING °6X,°°@ Smooth tube” °* 

PRINT USING “6%. °°4 High Flux ~** 

PRINT USING °6X,°°2 Turto-B°** 

PRINT USING °6X,°°3 High Flux Mod*”” 

PRINT USING °6X,°°4 Turbo-B Mod*** 

ENO IF 

INPUT Itt 

OUTPUT @FilelsItet 

ENO IF 

PRINTER IS 7@1 

IF Itt<!1@ THEN PRINT USING “16X,°°Tube Number: a Ort tat 
IF itt>9 THEN PRINT USING “16X.°"Tupe Numper: 2 Par ieke 


ne aL 


1600 
1672 
1682 
1692 
1782 
17G1 
1728! 
173! 
174Q) 
1752 
1762 
177 
1782 
179@! 
182Q! 
181@ 
1822 
1832 
184Q! 
18S@! 
1868 
1872 
1882 
1892! 
192Q! 
191@ 
1920 
1932 
194@! 
195Q! 
1960 


me 


Ire imme! Teen PRINi USING IJloa, iuoe type: ,$30 t1fdee cr 
SEEP 

INPUT “ENTER OUTPUT VERSION (@"LONG ,1*SHORT ,2"NONE)® Lov 

JEEP 

INPUT “SELECT (Q"LIQ ,PeVAP ,22(LIQ+VAP )/2)° Ilav 


CIMENSIONS FO TESTEO TUBES 

ELECTRIC HEATED MODE 

OisDiameter at tnermocouple positions 

DATA .Q111125,.83119125,.0111125,.0129542,.@12446 ,.2129542,.81200965 
DATA .@1@2@96S,.01157,.01157,.@1157,.8@1157,.@1157,.@18@@965 

READ Dia(e) 

OleOla(ltt) 


O2"Diametar of test section to the base of fins 

DATA .2@1587S,.@1587S5,.@1587S, .@15824, .@15875, .015824, .81272 
OATA .@127,.@138,.0@138,.0138,.0138,.@138,.8127 

READ O2a(e) 


ODielnsidce diameter of unanhanced ends 

DATA .@127;,8127,.0127,.0132,.0127,.0132, .Olllizs ore 
DATA .8118,.0118,.@118,.0@118,.@118, .@ttiti2s 

READ Diale) 


CoeOutside diameter of unennanced ands 

DATA .8@1587S,.0@1587S,.015875, .015824, .01587S, .@1S824, .e@1278, .01272 
DATA .O1331, 01331) ,.91331,.013351-201S6 02? 

READ Doa(e) 


LeLength of ennanced surface 
DATA .1016,.1016,.1016,.1016,.2032,.2032,.2032,.2032,.2032 22032 ..202 eee 


Sue One wea Ooe 


oo a 
1982! 
199Q! 
2222 


READ La(e) 


LueLangth of umennmancad surface at tha’ ends 
DATA .@254,.0254, .0254 ,.0254 , .O@762,,076¢2,.@76c ,.@762, .O0762,.0762,,0 acne 


b25.0'60,.0762 


2012 
2820! 
2232! 
2242 
2252 
2262 
2272! 
2082 ! 
2290! 
2122 
2118 
cleu 
21 3e 
2142 
2.50 
2162 


JR ref) rar 
RJ ps fl to 


READ Lua(e ) 


KousTmermal Conductivity of tuba 

DATA 398 344 344 45 344 45 344 344 398 398 398 398 ,398 344 
READ Kcua(e) 

IF Inmel THEN 


Data statemants for watar heating moda 


DATA @.21587S ,@.01587S ,@.0169 2.0138 2.0169 2 ,@,2@ ,2@,@ 
READ OD2ale) 

DATA @.@127 2.8127 ,@.814S 2.20127 0.2145 ,2,0,2,8,0 
READ Dia(e) 

DATA @.015875 ,@.@1587S5 @.0169 2.215875 2.0169 8,0,2,2,2 
READ Doa(e) 

DATA @.3248 2.3048 2.3048 0.3248 2.3048 2 ,2,2,2,2 
READ La(e) 

DATA @.0254 2.0254 2.8254 2.0254 2.0254 2,2 ,8,2,0 
READ Lua(e) 

CATA 398 45,398 ,45 ,398 2,.2,2,2 2 

READ Kcua(e) 

ENO IF 

O2=OZat ltt) 

DieDia( ltt? 

DoeDoal( ltt) 

LeLacItt? 

LuslLua( Itt) 


Keuskcua(l[ tt) 
Xne.§ 


dds2 


2308 Fre.3 

2310 IF Itt#@ THEN Cfe#!l.70E+3 
2320 IF Itt>@ THEN Cf#3.7037E+12 
23320 A=PI*(00*°2-01°2)/4 

2342 P#=PI*0o a 

2350 IF Imm#! THEN 


2362 8EEP 
23702 INPUT °TUSE INITIATION MOOE. (1*HOT WATER ,2=STEAM ,3*"COLO WATER)® ,Itim 


2380 IF Itim=! THEN PRINT USING “16X,"°Tube Initiate: Hot Water’”” 
2390 IF Itim=e2 THEN PRINT USING °!16X,°"Tube Initiate: Steam” 

2400 IF Itim=3 THEN PRINT USING “16X,°°Tube Initiate: Cold Water’** 
2412 INPUT “TEMP/VEL MOOE: (@=T-CONST ,V-OEC:1"T-OEC ,V-CONST; 2*T-INC ,V-CONST)°, 
Itv 

2420 IF Itv=@ THEN PRINT USING “16X,"°Temp/Vel Mode: Constant/Oecreasing’”* 
2430 IF Itv#l THEN PRINT USING “16X,°°Temp/Vel Mode: Decreasing/Constant”*” 
2440 IF Itv#2 THEN PRINT USING “16X%,°"°Temp/Vel Mode: Inereasing/Constant”** 
245@ INPUT “WANT TO RUN WILSON PLOT? (18#Y @8N)° [wil 

2462 IF Inme! ANDO Iwil#@ THEN 

2472 IF Itt#®@ THEN C:1#.032 

2480 IF Ittel OR Itt#3 THEN Ci#.259 

2490 IF Itt#2 OR Itt*4 THEN Ci®. 62 

2500 BEEP 

2512 INPUT °ENTER CI (OEF: WHe.032 HF=.0S59 7TB=.062)° Ci 

252@ PRINT USING “16X,°°Sieder-Tate “°° 

2530 PRINT USING °16X,°° Constant OLE a 3 Ct 

2540 ENO IF 

25S@ ENO IF 

2560 IF Inmel ANO Ime! ANO Iwile! THEN 


2570 IF Itt#@ THEN C1#.032 

2580 IF Itt#!l OR Itt#3 THEN Cr2#.2S9 
2590 IF Itt#2 OR Itt#4 THEN Ci*.062 
2622 ASSIGN @File2 TO ® 

2612 CALL Wilson(Cf ,C:i) 

2620 ASSIGN OFile2 TO O2_files 

2630 ENTER O@File2:Nrun ,Oold$ Ldtcl ,Ldtc2 ,Ldtc3 Latc4 Itt 
2648 ENO IF 

2650 Nsub#@ 

2650 IF Idp*4 THEN [hm=} 

2670 IF Inmel THEN Nsub#8 


2682 Ntc#6 

263@ IF Inme@ THEN Ntce#l2 
2700 Ja} 

2710 Sx#Q 

2720 Sys 

2730 Sxs#0 


2740) Sxy=@ 

2752 Repeat: | 

2760 IF Im=@ THEN 

2782 Ido#2 

2792 ON KEY 1,15 RECOVER 2752 

2802 PRINTER IS 1 

2810 PRINT USING °4xX,°°SELECT OPTION®’” 

2820 PRINT USING °&X ,°"°Q*eTAKE DATA’’® 

2830 IF Inme@ THEN PRINT USING “6X ,°"1*SET HEAT FLUX*** 
2840 IF Immel THEN PRINT USING “6X ,°"1@SET WATER FLOW RATE*** 
285@ PRINT USING °&X ,"°°2=sSET Tsat’’* 

2860 PRINT USING °4X,°°NOTE: KEY | = ESCAPE’*” 

2870 S8EEP 

2882 INPUT Ido 

2890 IF Ido>2 THEN Ido#2 

2900 IF Ido#=® THEN 4442 

2310! 

29320! LOOP TO SET HEAT FLUX OR FLOWMETER SETTING 


2930 IF Idoel THEN 
2348 IF Ihm=d@ THEN 


Irs 


2950 
2960 
2o7@ 
2982 
293 
3282 
3818 
3828 
3030 
3048 
3858 
3058 
3878 
3888 
3098 
389] 
3892 
3188 
3118 
5128 
313 
3148 
515e 
3168 
317 
3188 
3198 
32Q@ 
3218 
3228 
3238 
3248 
3258 
3262 
3278 
3282 
S29@ 
3388 
2310 
5526 
5330 
3348 
335 
3368 
3372 
3388 
3330! 
348! 
3418 
3428 
3438 
3448 
3458! 
3468 
3478 
3488 
3498 
3588 
3518 
3528 
3538 
3548 


OUTPUT 7@9);°AR AFZ2@ AL21 URS” 

BEEP 

INPUT “ENTER DESIRED Qdo° ,Oaae 

PRINT USING “4x ,*°DESIRED Qap ACTUAL Qap’’* 
Erre!8Qe 

FOR Ie! TO 2 

OUTPUT 789;:°AS SA” 

Sum=@ 

FOR Ji#! TO S 

ENTER 7@091E 

Sum=#SumerE 

NEXT Ji 

IF Ie! THEN Volt#Sum/S 

IF I1*#2 THEN Amp=Sum/5 

NEXT I 

PRINT °Vs:* Volt 

PRINT “Is:° Amo 

Amo2ABS(Amoe!]. 9182) 

Volt#ABS(Volte2S) 

Aado@VolteAmp/(P1«O2¢L ) 

IF ABS(Aqde-Oqda)>Err THEN 

IF Aqda>Dadp THEN 

BEEP 400@,.2 

BEEP 400@,.2 

BEEP 4008,.2 

ELSE 

BEEP 252,. 
BEEP 2S5Q,. 
BEEP 258,. 
END IF 
PRINT USING °4X .MZ.30E ,2x .MZ.30E° sDade ,Aqde 
WAIT 2 

GOTO 3888 

Svoc = 

BEEP 

PPLINT USING “4% MZ.30E 2x .MZ.30E° ;Oqda wade 
Err #580 

WAIT 2 

60TO 30808 

ENO IF 

ELSE 

BEEP 

INPUT “ENTER FLOWMETER SETTING” .Fms 

60TO 2812 

ENO IF 

END IF 


hI hI PQ 


LOOP TO SET Tsat 


IF Ido#2 THEN 

IF Ikdt#! THEN 3472 

BEEP 

INPUT "ENTER DESIREO Tsat” ,Otla 

PRINT USING °4xX,°" OTsat ATsat Rate Tv Rate’** 
Ikdt=1 

Old!=@ 

Old2=@ 

Nn=1 

NeseNn MOO 15 

Nn#Nnel 

IF Nese#l THEN 

IF Inm*2 THEN PRINT USING °4X °° Tsat Tlal Tld2 Tv 
IF Imme! THEN PRINT USING °4xX,°"° Tsat Tldl raz Tv 


t Loade” lout ~~ 


3550 


3560 
3578 


END IF 
IF Imme@ THEN OUTPUT 789;°AR AFO@S ALT! VRS” 
IF [mmel THEN OUTPUT 789;°AR AF@ ALS VRS" 


alee 


Tsump 


Tsumo 


Tinle 


3582 
Bode 
3622 
3612 
3622 
3632 
36542 
3652 
3662 
3672 
3682 
3692 
3702 
3712 
3720 
3732 
3742 
3750 
3762 
3770 
3782 
3792 
3802 
3812 
3822 
3832 
3842 
3850 
3862 
3872 
3882 
3892 
3900 
3912 
S320 
3932 
3942 
3952 
3962 
3972 
3982 
39392 
4202 
4212 
4022 
4232 


FOR lei TO 6 

IF Ihme® AND I1>4 THEN 3862 
Sume=2 

OUTPUT 729;°AS SA” 

FOR Ji#!l TO 22 

ENTER 709:Eliaq 
Sum=Sum+Eliaq 

NEXT Ji 

Eliqe#Sum/22 
TLld@FNTvsv(Eliaq) 

IF I#] THEN TidleTld 

IF I=2 THEN Tld2=Tld 

IF I#3 THEN TveTld 

IF [#4 THEN TsumpeTla 

IF 1eS THEN Tinlet#=Tld 

IF I1#6 THEN TouteTlad 

NEXT I 

IF Imhm=! THEN 

OUTPUT 729;"AR AFQ@Q ALQ@Q@ VRS* 
OUTPUT 789:°AS SA” 

Sum= 

FOR Kke! TO 22 

ENTER 789;:E 

SumeSum+E 

NEXT Kk 

Emf (7 )=ABS(Sum/2Q) 
TpoilesEmf(7)/3.96E-4 

END IF 
Atide(Tldl+Tld2)e.5 

IF ABS(At]d-Dtld)>.2 THEN 
IF Atid>Dtld THEN 

BEEP 4200,.2 

BEEP 42002,.2 

BEEP 4202,.2 

EGSe 

BEEP 252,.2 

BEEP 252,.2 

BEEP V2sO-.2 

END IF 

ErrieAtid-Old! 

Oldt=Atid 

Erer2eTv-0ld2 

Old2eTv 

IF T1ld1>100. THEN 4262 

IF IhmeQ THEN PRINT USING °4X S(MDDD.D0D ,.2X)°:DOtld,Tld! ,Tld2,Tv,Tsumop 
IF Inmel AND Idp#@ THEN PRINT USING °4X ,7(M00.00 ,2xX)*10tld Tldl ,Tid2,Tv,Ts 


ump ,tTinlet Toile 


4042 


IF Immel AND Idp#4 THEN PRINT USING °4X ,S(MDD.00 2X) ,3(M30.00 ,2x)"1s0tld,Tl 


d! ,Tld2 ,Tv ,Tsump ,Tainlet ,Tpile,Tout 


4252 
4062 
4272 
4282 
4232 
4122 
6112 
4122 
4130 
4142 
4152 
4162 
4172 
4162 
4192 
4222 


WAIT 2 

60TO 3500 

ELSE 

IF ABS( Atl d-Dtid)>.1 THEN 
IF Atld>Dtld THEN 

BEEP 3222,.2 

BEEP 3002,.2 

ELSE 

BEEP 822,.2 

BEEP 822,.2 

END IF 

ErrieAtld-O0ld! 

OldteAtld 

Err2eTv-01ld2 

OldgeTv 

IF Imme2 THEN PRINT USING °4X S(MO00.00 .2x)°:s:Otla,Tldl ,Tld2,Tv,Tsump 


18g Bs 


4210 IF Immel THEN PRINT USING °4@x ,S(¢MO0.00,2x) ,5(450-00 (14 frei ee 
~Tsump ,Tiniat,ipise,iout 
22@ «WAIT 2 
4230 60T0 3522 
4248 ELSE 
4250 BEEP 
426Q@ €ErrieAtld-Old!) 
4270 Oldt#Atia 
4280 Err2eTv-Old2 
42392 Old2sTv 
4300 IF Ihm=@ THEN PRINT USING °4X ,S(M000.00 ,2%)"°:Otld,Tldi ,Tld2 Tv ,Tsumo 
4310 IF ITmme!l THEN PRINT USING °4X ,8(M00.00 ,2X)°s:O0tld,Tld! ,Tld2,Tv ,Tsumo ,Tinlet 
thi Le tout 
4322 WAIT 2 
4332 GOTO 3522 
4340 END IF 
435@ ENO IF 
4362 END IF 
4370! ERROR TRAP FOR Ido OUT OF BOUNDS 
438@ IF Ido>2 THEN 
4390 BEEP 
4482 GOTO 2812 
4412 END IF 
442Q! 
4432! TAKE DATA IF Ime@ LOOP 
4440 IF lkolel THEN 4482 
4452 s8EEP 
4462 INPUT “ENTER BULK DIL %° Boo 
4470 Ikole! 
4482 IF Inm9@ THEN OUTPUT 789;:°AR AFOQ ALI VRS” 
4490 IF Inmmel THEN OUTPUT 729;:°AR AFQ ALS VRS* 
4502 IF Inmme@ THEN Ntc#t2 
4518 FOR Teli TO Ntc 
4520 OUTPUT 78@9;°AS SA” 
4532 Sumed 
4540 FOR Ji©#1 TO 22 
4552 ENTER 7Q9;E 
456@ SumeSum+& 
4570 IF [#(S$-Nsuo) OR [Te(!Q-Nsub) THEN Et(Ji-1 )e#£ 
4580 NEXT Ji 
4590 Kadl#Q 
46@@ IF I[e(S3-Nsub) OR 1e( 1Q@-Nsub) THEN 
4612 EaveesSum/22 
4622 Sum=Q. 
4630 FOR Jee@ T0 19 
4642 IF ABS(Et( Jk -Eave)<S.Q0E-6 THEN 
4652 SumeSumrEt( Jk) 
466@ ELSE 
4670 Kdlekdl+!) 
4682 ENO IF 
46392 NEXT Jk 
47282 IF I#(9-Nsub) OR Ie(1Q0-Nsub) THEN PRINT USING °4X,°°Kdl = ** OO°:Kdl 
4710 IF Kdl>1@ THEN 
472@ BEEP 
4732 BEEP 
4742 PRINT USING °4X °"°Too mucn scattering in data - repeat data sat*”” 
4752 60TO 2802 
476@ ENO IF 
4772 ENO IF 
478Q@ Emf(] )eSum/(28-Kdl ) 
4792 NEXT I 
4822 IF Ihmel THEN 
4812 OUTPUT 709: °AR AFQQ ALOO VRS” 
4822 OUTPUT 729:°AS SA° 
4832 Sums 
4842 FOR Kkel TO 22 


4éS2 
4a5ed 
4872 
4880 
4832 
4922 
4912 
4920 
4932 
4942 
4950 
4960 
4972 
4982 
4992 
5220 
5212 
$222 
5232 
5042 
5250 
5262 
5270! 
5280! 
52392 
518 
5112 
5120 
$132 
5140 
5152 
5160 
5172 
5182 
5192 
5222 
5212 
5220 
5250 
5242 
5252 
Sy 
5260 
5272 
5282 
5292 
5322 
5312 
$322 
S330 
5342 
5352 
5360 
$372 
5382 
5390 
5400 
5410 
5422 
5430 
5442 
$452 
5460 
5472 
5482 
54932 


ENTER 7OS:5 

SumMssumee 

NEXT Kk 

Em? (7 )#ABS(Sum)/22 

ENO IF 

IF ImmeQ@ THEN 

OUTPUT 709;"°AR AF2@ ALZ2! VRS“ 
FOP I#!1 TO 2 

OUTPUT 70S AS SA 

Sum=2 

FOR Jisl TO 5 

ENTER 7@09;E 

Sume@SumeE 

NEXT Ji 

IF I#1 THEN VreSum/§ 

IF [#2 THEN Ire#Sum/5 

NEXT J 

ENO IF 

eect 

IF Ihm=Q@ THEN ENTER @File2:Bop ,Told$ ,Emf(e) Vr Ir 
IF Inmet THEN ENTER @File2:Bop ,Told$ ,Emf(e) Fms 
ENO IF 


CONVERT emf'S TO TEMP ,VOLT ,CURRENT 
Twa2d 
FOR I#i TO Nte 
IF Idgtce>@ THEN 
IF I=tdtc!) OR IeLdtc2 OR IelLdtc3 OR IelLdtc4 THEN 
TCT )#-99.99 
60TO $240 
ENO IF 
ENO IF 
IF Itt<4 ANO Imm=@ THEN 
IF I1>4 ANO I<93 THEN 
TCT 2-99.99 
GOTO $242 
ENO IF 
ENO IF 
TCT @FNTvav(Emf( I) ) 
NEXT I 
IF Itt<4 THEN 
WARNING: 4 BAO THERMOCOUPLES ARE NOT POSSIBLE FOR ITT<4 
FOR I=! TO 4 
IF IeLdtcl OR IeLdtc2 THEN 
Twaelwa 
ELSE 
TwaeTwatT( 1) 
ENO IF 
NEXT I 
TwaeTwa/(4-Idtc) 
ELSE 
IF Inhme=i THEN $452 
FOR [=i TO 8 
IF Ieidtc! OR Ietdtc2 OR Ietdte3 OR IelLdtc4 THEN 
Twaelwa 
cece 
TwaeTwatl( 1) 
ENO IF 
NEXT I 
TweTwa/(8-Idtc) 
ENO IF 
Tla=@T( 9=Nsub ) 
Tld2=T( 1@-Nsub ) 
Tidae(TlatTld2)¢.5 
Tv#T( 11-Nsub) 
IF Itt<3 ANO Ihm=Q@ THEN 


de OF 


S502 Tld2e-99,99 

SSIQ Tve(T(1@)47( 11) 72 
S52@ END IF 

SS3@ Tsumpel( 12-Nsub ) 
S540 IF Inme@ THEN S572 
SSS@ TinlateT(13-Nsub) 
SS6Q@) —-TouteT(14-Nsub ) 
S570 IF Inme#@ THEN 

S580 Amp#A8S( Ire! .9182 ) 
SS$9Q@ VolteA8S(Ur)e25 
S600 QeVolteAmp 

S610 END IF 

562@ IF Itt*@ AND Imme@ THEN 
S630 KcueFNkKcul Tw) 


S64@ ELSE 

S6S@ Kcuekcual( Itt) 
S66@ END IF 

S67Q! 


S68@' FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
S690 IF Imme@ THEN TweTw-QelL06(02/01 )/(2ePleKcueLl ) 
S702 IF Ilave@ THEN TsateTlda 

S71@ IF Ilavel THEN Tsate(TldatTve.5 

S$72@ IF Ilaeve2 THEN TsateTv 

S73@ IF Immo! THEN 

S748 TavgeTinilet 

S75Q@ 6rade37.9853+.104388*Tavg 

S76Q@ Tdrop#A8S(Emf(7))e1 .E+6/(18¢6rad) 

S770 TavgceTinlet-Tdrope.S 

S78Q@ IF A8S(Tavg-Tavgc)>.@! THEN 

S790 Tavge(Tavet+Tavge )*.5 

S802 60TO 5752 

S8!1@ END IF 

582! 

Sa3@' COMPUTE WATER PROPERTIE 

S342 IF Imme! THEN 

S850 KwerFNKw( Tavg) 

S86Q@ MuwaseFNMuwl! Tavg ) 

S87@ CoweFNCow! Tavg) 

S88Q@ PrweFNPrwulTavg) 

S89@  RnoweFNRhow! Tavg?) 

5900 Twielavg 

$919! 

S920! Compute MODT 

S$93@ Mdote3.96S7E-3+Fmse¢( 3,61955E-3-Fmse¢(8.820Q06E-6-Fmae( 1. 23688E-7-Fmse4.31897 
E-10))) 

5940! MdoteMdote(1.@36S-Tinlete(!.96644E-3-TinleteS.252E-6 ))/1.8037 
S$95Q@ Kdt#@ 

S36Q@ QeMdot*CoweTdrop 

S970) LetdeTdron/LOG((Tinlet-Tsat)/(Tinlat-Tdrop-Tsat)) 
S98Q@ VUoeQ/( Pl eDoeLeLmtd) 

S99Q@ RweDoelLD6(Do/D1)/(2. ekcu) 

6@@0 Twelsatt+FreLatad 

B@I@ VweMdot/(RmowePleOi*2/4) 

6@2@ ReweRnoweVweD1/Muwa 

6@30 HieCieku/DieRaw® .8ePrwu (1/3. )°( Muwa/PNMuw( Twi) )*.14 
6040 TwiceTavg-Q/( PI eDoeLeH ) 

6850 IF ABS(Twi-Twic)>.@1 THEN 

6@6Q Twrel(TwitTwicde.§ 

6878 6GOTD 6232 

628@ END IF 

6990 Twre(TwreTwicde.5 

6108 Hoe#l/( 1/Uo-Do/( Dierks “Ry? 

6'11@ END IF 

612@ ENO IF 

6130 IF Immel THEN 

6140 TmetaneQ/(HoePl eDoeL ) 


dies 


6152 
6162 
6172 
6182 
5198 
6200 
6212 
6222 
6232 
6240 
6250! 
6260! 
6272 
6282 
6292 
6322 
6312 
6322 
6332 
6342 
6352 
6360 
6372 
6380! 
638) 
6382 
6383 
6384 
6385 
6386 
6387 
6388 
6390! 
6420! 
6412 
6422 
6432 
6442 
6452 
6462 
6472 
6482 
6492 
65202 
6512 
6511 
6S2Q! 
6530! 
6542 
6552 
6562 
657! 
6580! 
6592 
662 
6612 
6620! 
663! 
6642 
6652 
6662 
6672 
6682! 
6690! 
6722 
6712 


Twwisateinetao 

ELSE 

Thetabelw-Tsat 

ENO IF 

IF Thetab<@ THEN 

BEEP 

INPUT “TWALL<TSAT (Q@=CONTINUE, 1*ENO)* , lev 
IF ITeve® THEN GOTO 2760 

IF lev=! THEN 7222 

ENO IF 


COMPUTE VARIOUS PROPERTIES 
Thilme(TweTsatde.S 
Rno#FNRno(Tftitlm) 
MueFNMul( Thilm) 
KeFNK(TFilm) 
CpoeFNCpiTfilm) 
BetasFNBeta(Tfilm) 
Hf QeFNHfg(Tsat) 
Ni#Mu/Rno 
Alpnaek/(RnoFfCp ) 
PereNi/Aloha 
PsateFNPsat(Tsat ) 


PRINT "Rno:” ,Rno 

PRINT “Mu: * Mu 

PRINT “°K: °K 

PRINT °Cp:° Co 

PRINT “Tfilm:* ,Tfilm 

PRINT ‘Hfg:* Hfg 

PRINT °“Psat:* Psat 

PRINT “Tsat:° ,Tsat 

COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
FOR UNENHANCEO ENO(S) 

Hoare! 3 

Fe=(HpareP/(KcueA) )*.SeLu 
Tann@#FNTann( Fe) 

Tneta=TnetabdelTann/Fe 

Xxe(3.8) eBetaeThetans0o*3eTann/(FeeNi*Alpna) )*.166667 
Yye( )4(,559/Pe )*(9/16))°(8/27) 
Hoarcek/Oo*( .6+,.387¢Xx/Yy)*2 

IF A8S((Hbar-Hbare )/Hbarc)>.00) THEN 
Hoare(HoartHbarc )¢.5 

GOTO 6422 

END IF 

PRINT “NAT CONV h:* |Hbar 


COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENO(S) 
Qle(HbarePeKcueA)* .SeThetabeTann 

Qe#Q-24Q1 

AsePleQ2eL 


COMPUTE ACTUAL HEAT FLUX ANO BOILING COEFFICIENT 

Qdp=Qc/As 

Htupe*Qdp/Thetab 
Caf=(CpoeTnetab/Hfg)/(Qdp/( MueHf de! .014/(9.81 oRho )*.5)7(1/3. Pre *1.7) 


RECORD TIME OF OATA TAKING 
IF Im=@ THEN 

OUTPUT 709;°TD° 

ENTER 783; Told$ 

ENO IF 


OUTPUT OATA TO PRINTER 


PRINTER IS 721 
IF Iov#@ THEN 


LAE, 


6720 PRINT 

6730' PRINT USING °10x ,°°Data Set Numper = ~~ ,000 2X ,°“Sulk O17) 2 se mop ee 

4A°:sJ 80p ,Told$ 

6731 PRINT USING °10X,°°Oata Set Numper © °° 000,2X,°°Bulk O11 % = “* 00.0°sJ,8 

op 

6732 PRINT ° TIME: ° ,TIMES( TIMEDATE ) 

6740 IF Inm=Q@ THEN 

6750 PRINT USING *18X,°°TC No: | 2 3 4 S 6 7 

gree 

6760 PRINT USING °10X,°"°Temp :°°* 8¢1x% MOD.00)°:T(1),702) 103) T4010 5) eee 

7),1¢8) 

6778 PRINT USING °18X,°° Twa Tiiad Tliad2 Tvapr Psat Tsump” °° 

6780 PRINT USING °10X ,2(¢M00.00,1X),1x ,M00.00,1x ,2¢( 1X ,MOO.00) 2% AOD VD); cee 

ld2,Tv ,Psat ,Tsump 

6790 PRINT USING °10X,°° Thetab Htube Gda >” 

6800 PRINT USING °10X .MO0.30,1X ,MZ.30E ,1X .MZ.30E°;: Thetab ,Htube Qdp 

6810 ELSE 

65828 PRINT USING "10X,°° Fes Vid Tsat Tinl Terop Tnretab Q Uo 
Hott 

6830 PRINT USING °10X ,4(20.00,1X),2.30,1X ,00.00,1X% ,3(MZ.30E ,1X)°sFms Vw Tsat fi 

nlet ,Tarop ,Tnetaod ,Qap ,Uo ,Ho 

6840 ENO IF 

6850 ENO IF 

6860 IF lTovel THEN 

6870 IF J=#!l THEN 

6880 PRINT 

6890 IF Ihm=@ THEN 

6300 PRINT USING "18X,°" RUN No Oils Tsat Htube Qdp Thetab’”” 

63910 ELSE 

692@ PRINT USING °10X,°° FMS 8). Se eal HTUBE QOP THETAB”*°" 

6332 ENO IF 

6940 ENO IF 

6352 IF Ihm=*@ THEN 

63960 PRINT USING "12X ,30,4X ,00,2X ,MO00.00 ,3( 1X MZ.30E)°sJ ,80p ,Tsat ,Htube ,Qdp,The 

tad 

69370 ELSE 

6980 PRINT USING °12X,30,4X ,00,2xX ,M00.00 ,3( 1X ,M2.30E)°:Fms 800 ,Tsat tude ,Qan,T 

netad 

6390 ENO IF 

7080 ENO IF 

7010 IF Im=#@ THEN 

7020 S8EEP 

7030 INPUT “OK TO STORE THIS DATA SET (18#Y ,@9N)?° Dk 

7040 ENO IF 

7050 IF Okel OR Ime! THEN JeJ+l 

7060 IF Okel ANO ImeQ THEN 

7070 IF Inhm=@ THEN DUTPUT @FilelsBoo ,Told$ ,Emf(e) Vr Ir 

7080 IF Inmet THEN OUTPUT @Filels8on ,Told$ ,Emfie) Fms 

7038 ENO IF 

7100 IF Iufel THEN OUTPUT ®@UfilesVu Uo 

7110 IF Ires#l THEN OUTPUT @RefilesFms Rew 

7120 IF (Im=1 OR Ok#!) ANDO Iplot#! THEN OUTPUT @PlotsQdp ,Thetab 

7130 IF Ime® THEN 

7148 SEEP 

71S8@ INPUT “WILL THERE 8E ANOTHER RUN (1#Y ,08N)?° ,Go_on 

7160 NruneJ 

7178 IF Go_one® THEN 7220 

7188 IF Go_on<>® THEN Repeat 

7190 ELSE 

7200 IF J<Nrunel THEN Repeat 

721@ ENO IF 

7220 IF Im#@ THEN 

1250 GEEP 

7240 PRINT USING “10X,°°NOTE: °° .22Z,°° cata runs were stored in file ~° 108A iJ- 

1 O2_f ied 

7250 ASSIGN OFilel TO 


26 


7260 OUTPUT O@File2iNrun-! 

7278 ASSIGN OFilel TO Ol_files 

728@ ENTER OFile!:DateS Ldtcl ,.dtc2 [tt 
7290 DUTPUT OFile2:0ate$ ,Latc! ,.Ldtc2 Itt 
7300 FOR [#1 TO Nrun-l 

7310 IF Inm*@ THEN 

7320 ENTER @Filel:B0p ,Told$ ,Emf(e*) Ur Ir 
7330 OUTPUT @File2iBop ,Told$ ,Emfle) Ve Ir 
734@ ELSE 

7350 ENTER O0File!;:Bop ,Told$S Emfie) Fms 
736Q@ OUTPUT @File2:Bop ,Tola$S ,Emfie) ,Fms 
Po70) ENG IF 

7380 NEXT I 

7390 ASSIGN @Filel TO « 

7400 PURGE “OUMMY" 

7410 END IF 

742@ BEEP 

7430 PRINT 

7442 IF Iplot#! THEN PRINT USING °10X ,""NOTE: “° ,22,°° X-Y pairs were stored in 
plot date file ~" ,10A°;J-I ,P_files& 

7450 ASSIGN @File2 TO * 

7462 ASSIGN OPlot TO « 

7472 IF Iuf#l THEN ASSIGN ®Ufile TO * 
7480 IF Ire#1 THEN ASSIGN @Refile TO 
7490 CALL Stats 

7500 BEEP 

7510 INPUT “LIKE TO PLOT DATA (18Y ,@9N)7° Ok 
752@ IF Ok#!l THEN CALL Plot 

7530 SUBEND 

754! 

7550! CURVE FITS DF PROPERTY FUNCTIONS 

756@ DEF FNKcu(T) 

7570! OFHC COPPER 258 TD 300 K 

758Q@ TkeT+273.15 Geo” K 

7590 Kea3a-.112eTk . 

76@@ RETURN K 

761@ FNENO 

762@ DEF FNMU(T ) 

Tec? ‘COM /larop/ Ift 

7622 IF Ift*#@ THEN 

7630! 178 TO 368 K CURVE FIT OF VISCOUSITY 


7640 TkeTt+273.15 lee to kK 
765Q MueEXP(-4.46364(1011.47/TK ) 01. QE-3 
fool SEND IF 


7652 IF Iftel THEN 

gees!) 225 30 3573 K 

7657 Muses. 4993-4. 26E-2¢T+3.948SE-40T* 2-4, 193E “Ge T° 342. 0709E-8eT 4 
7658 MusMuel.@E-4 

7659 ENO IF 

7660 RETURN Mu 

76870 FNEND 

768@ DEF FNCo(T) 

7681 COM /Iprop/ Ift 

7682 IF Ift#@ THEN 

7690! 188 TO 400 K CURVE FIT OF Co 

7700 4 TkeT+273.1S te TO % 

7710 Cp=#.4018841 .GSQQ7E-ZeTk+! SI494E-SeTk “2-6.678S3E-10¢eTk “3 
7720 Co=Co*1222 

7721 ENO IF 

7722 #IF Ift#=1 THEN 

7723! 223 10 373 kK CURVE FIT Co(Rt2@) 

7725) «=Cpel.054242.1405E-ZeT+! .O709E-SeT* 0-6 . 4721E-B8eT*3-1.4324E-JeT*4+4,.136E-11 8 
T*S 

7726 Co#=Cp #1202 

G72? ~END IF 

7730 RETURN Cop 


Ze 


7742 FNEND 

775@ OEF FNRno(T) 

7751 COM /iprcep/ Ift 

77S2 #IF Ift#@ THEN 

776@ Tk=sT+273.15 i¢ tC 

7770 Xwl-(1.8eTK/753.9S) 'K TOR 

778Q@ Ro=36.32+61.1464140X°(1/3)016.41801SeX+17. 476838 exX* Se] .119828ex*2 
779@ Ro=Ro/ .062428 

7791 ENO IF 

7792 «IF Ift=#=) THEN 

7793! 223 TO 373 K (R-124) 

7796 Rowl434.8-2.8619¢T-6.7267E-3eT°2-7.2852E-SeT *3 
7797 ENO IF 

7800 RETURN Ro 

7818 FNENO 

7820 OEF FNPr(T) 

783Q@ Pr=FNCo(T )eFNMu( T)/FNK(T) 

764@ RETURN Pr 

78SQ@ FNEND 

7860 O&F FNK(T) 

7861! COM /Iproo/ Ift 

7862 IF Ift=@ THEN 

7870! T<36@ K WITH T IN C 

788@ K=.@71-.@0@261 eT 

7881 ENO IF 

78682 IF Ift=! THEN 

(985! 223 10 373 K CR=128) 

7B8B8S Kae7.SI91E-2-3.S436E-4eT-!) .994S5E-7eT°243.183S5SE-JeT*3 
7886 END IF 

789@ RETURN K 

79@@ FNENO 

791@ OEF FNTann(X) 

7920 P=sEXP(X) 

793Q@ Q=1/P 

7942 Tanhn=(P-Q)/(P+4Q) 

795@ RETURN Tanh 

796@ FNENO 

797@ OEF FNTvsviV) 

7980, COM /Ce/ C17) heal 

799@ T=C(@) 

8@@@ FOR I=! TO 7 

B8@1@ TsTe+C( 1 )ev*l 

8@2@ NEXT I 

8030 IF Ical=! THEN 

8040 TxT-6§.7422934E-2¢T¢(9.0277043E-3-Te( -9.325991 7E-S)) 
8@S@ ELSE 

8@60@ T=T*8.626897E-2+Te¢(3.76199E-3-TeS.Q6892S9E-S ) 
8@7@ ENO IF 

8@8@ RETURN T 

8@9@ FNENO 

81@@ OEF FNBeta(T) 

811@ Rop=sFNRno( T+. 1! ) 

812@ Rom=FNRno(T-. 1) 

8132 Setaz-2/(RoptRom)e(Rop-Rom)/.2 

8140 RETURN Seta 

81S@ FNENO 

816@ OEF FNHFQ9(T) 

8161 COM /Iprop/ Ift 

8162 IF Ift=@ THEN 

8170 Hf gm] 3741 344E4S-Te( 3. 3094361 E+2e7¢! 2165143) 
817) ENO IF 

8172 IF Ift=! THEN 

Bl73! 223 TO 373 K CURVE FIT (RI2Z4) 

8176 Hfg=1S9.93-.4609eT-1.4S58E-3¢T* 2-1 .371SE-SeT*3 
8177 Hf g=Hfge122@ 

8178 ENO IF 


22 


8182 
8198 
8280 
8221 
8202 
8219! 
8222 
8230 
8242 
824) 
8242 
8243! 
8245 
8247 
8248 
8250 
8260 
8272 
8282 
8290 
8332 
8310 
8320 
8330 
8342 
8352 
8360 
8370 
8380 
8332 
8402 
8419 
8420 
84390 
8449 
8458 
8460 
84790 
8482 
8492 
8500 
8512 
8520 
8530 
8542 
8550 
8560 
8572 
8582 
8530 
8620 
8610 
8620 
8632 
8642 
8652 
8660 
8672 
8682 
8690 
8720 
8712 
8722 
8732 
8742 
8750 


RETURN Hfg 

PF NEND 

DEF FNPsat (Tc) 

EOM lores, lft 

IF Ift#O THEN 

@ TO 8@ deg F CURVE FIT OF Psat 

Tfe!l .8eTce+32 

PaeS .94552S+TFe( .153552082+T£e( 1. 4840963E-3+TF eS. 61S0671E=6) ) 
PgePa-14.7 

END IF 

IF Ift#!l THEN 

223 10 373 Kk CR=124) 
Pae(162.68+5.946eTc+9.2081E-29Tc* 246. 9359E-4eTc*3 0. 14504 
Po=Pa-14.7 

ENO IF 

IF Pg>@ THEN | +=ePS16,-e1n Ho 

PsatePg 

Ezoc 

PsatePoe#29.92/14.7 

ENO IF 

RETURN Psat 

FNEND 

DEF FNHsmooth(X ,80p ,I sat) 

DIM ACS) ,8¢(5),C°5),0¢5) 

DATA .20526 ,.29322,.319048 ,.§5322 ,.79909,1.00258 
DATA . #4515 ,.72992,.73189,.71228, 68472 ,.64197 
DATA .41092,.17726,.25142,.54806 ,.81916 ,1.8845 
DATA .71403,.729!13,.72565,.69669!1 ,.66S5867 ,.61889 
READ Ale) BC) C(*) Die) 

IF Bop<6& THEN I*Bop 

IF 8op*6 THEN [#4 

IF 8o0pe!10 THEN [#S 

IF Isatel THEN 

HS@EXP( ACI )+8( 1 )*LO06(X)) 

ELSE 

HS@EXP(C(1)#D( 1 )*#LOG(X)) 

END IF 

RETURN Hs 

FNEND 

DEF FNPoly( x) 

Gome-Gely, ACO, 10) ,CC10) 8CS) Nop Jeent Ope .llog,ifn,ljorn Nyon 
Xl =x 

Poly#B8(@) 

FOR Ie! TO Nop 

IF Iloge!l THEN X1l*#LO6(X) 

PolysPoly+B( I )eX1*!I 

NEXT I 

IF llogel THEN Poly#EXxP(Poly ) 

RETURN Poly 

FNEND 

SUB Poly 

DIM RC1Q) $018) Sy(12),Sx612),Xx( 180), ¥y( 100) 
EORMeyCaly/ ‘acl .18) CC ie) 8(5) N,lornt Ope ,llog,Ifn,ljoin,Njoin 
CDM /Xxyy/ Xp(25) ,Yp(25) 

FOR I#@ TD 4 

8( 1 )#Q 

NEXT I 

BEEP 

INPUT “SELECT (Q*F ILE ,1eKEYBSDARD ,2=#PROGRAM)” Im 
ImelIme) 

BEEP 

INPUT “ENTER NUMBER OF X-Y PAIRS" Np 

IF Ime@)] THEN 

SEEP 

INPUT “ENTER DATA FILE NAME“ ,O_ files 

BEEP 


eae 


8760 INPUT “LIKE TO EXCLUCE DATA PAIRS (lev Q@8N)?”" Jeg 
8772 IF led=l THEN 

8780 S8EEP 

879@ INPUT “ENTER NUMBER OF PAIRS TO BE EXCLUDED” ,Ipex 
8802 ENO IF 

8810 ASSION OFile TO O_files 

8820 ELSE 

883@ S8EEP 

8840 INPUT “WANT TO CREATE A DATA FILE (187 ,@8N)?° Yes 
8850 IF Yese#! THEN 

886@ S8EEP 

8870 INPUT “GIVE A NAME FOR OATA FILE® ,0_files 
8880 CREATE 8DAT O_file$,5 

8892 ASSIGN @File TO O_files 

8922 ENO IF 

891@ ENO IF 

8922 S8EEP 

8930 INPUT “ENTER THE ORDER OF POLYNOMIAL” ,N 
8940 FOR [#@ TO Ne2 

8958 Sy(I)#@ 

8960 Sx(I)#@ 

8972 NEXT I 

8982 IF lede!i ANO Ime! THEN 

8990 FOR I#!1 TO Ipex 

9000 ENTER @Filesx,yY 

9010 NEXT I 

9822 ENO IF 

9032 FOR I#! TO No 

9040 IF Im=l THEN 

9050 IF Dpo#2 THEN ENTER ®File:X,Y 
9060 IF Qpo<2 THEN ENTER OFilesyY,Xx 
9070 IF Dpo#l THEN Yer/x 

9080 IF llogel THEN 

9290 IF ODpo#e2 THEN Xtex/Y 

9122 XetQ6(Xx) i 

9110 IF Opo#2 THEN YeLOG( xt) 

9120 IF Opo<2 THEN YeLD6( Y) 

3130. END IF 

91428 ENO IF 

915Q@ IF Im=2 THEN 

9160 8EEP 

9170 INPUT “ENTER NEXT X-Y PAIR® X.Y 
9162 IF Yesel THEN OUTPUT @File:x,Y 
$19@ ENO IF 

9200 IF Im<3 THEN 

9210 Xx(I)ax 

9220 yy I dey 

9230 ELSE 

9240 xXeXxp(I-1) 

9250 YeYp(I-1) 

926@ ENO IF 

9270 R(Q)#Y 

9280 Sy(Q@reSy(Q)+yY 

3290 54 18x 

93QQ Sxl eSx( 1 +x 

9312 FOR Jel TON 

9320 RJ eR(J-1 dex 

9330 Sy( J deSy( J d+R( J) 

9340 NEXT J 

9350 FOR J=2 TD Ne2 

9360 St J eS( Jal dex 

A370) Sat LyesxGl est) 

9380 NEXT J 

9390 NEXT I 

9420 IF Yes#!l AND Ime#2 THEN 

9412 S8EEP 


124 


9420 PRINT USING “124,UL, ~ X~Y Dairs were stored in ".ie (10M GND U_tiles 


$430 ENO IF 

9440 Sx(Q@)=@Np 

945@ FOR I*#@ TON 

946@ Ctl eSyi1) 

9470 FOR Je@ TON 

94BQ ACI, J )eSx(IeJ) 

9490 NEXT J 

95@@ NEXT I 

951@ FOR I[*@ TO N-1 

9520 CALL Divide(l) 

953@ CALL Subtract( I+!) 

9540 NEXT I 

95S@ B(N)#@C(N)/ACN,N) 

9560 FOR I=#@ TO N-I 

9570 B(N-I-I)#C(N-1-1) 

958@ FOR J=@ TO I 

959@ B(N-I-I )@B(N-I-I)-ACN-1}-1 Nod *BCN=J) 
9600 NEXT J 

9610 BCN-I-I @B(N=1-1)/ACNH1<I N-I AT) 

9620 NEXT I 

9630 ‘PRINTER IS 701 

964@ !PRINT Bie) 

S65@ ' PRINTER IS 70S 

S66@ IF Iprnt*@ THEN 

9670 PRINT USING °12X,° "EXPONENT COEFFICIENT*"° 
96B8@ FOR I*#@ TO N 

9690 PRINT USING “15x ,00,SX MO. 70E*:1 ,B( 1) 
9700 NEXT I 

971@ PRINT ° “ 

9720 PRINT USING “12x ,°°OATA POINT x 7 Y(CALCULATEO) OISCR 
EPANCY°"° 

973@ FOR I=! TO No 

9740 Ye=xB(Q@) 

9750 FOR J=! TON 

9760 YeseYer+B( J exx(1)*J 

9770 NEXT J 

97B@ OsYy(l)-Ye 

9790 PRINT USING °1SxX ,30 4X ,4(M0.SDE ,1K)"31 ,Xe(1) Yy(1),Ye ,0 
SBQ@Q@ NEXT I 

981@ END IF 

9B20 ASSIGN @File TO « 

983@ SUBENO 

9840 SUB Oivide(M) 

desO) COMe7Coly/ ACI@,16) CC 10) BCS) N,larnt ,Opo,llog.lfn,ljoin,Njoin 
9B6@ FOR I=M TO N 

9870 Ao#ACT ,M) 

968B@ FOR J=m TO N 

9B89@ ACI J ACI ,J)/Ao 

9900 NEXT J 

9910 Cl )#C(l)/Ao 

99320 NEXT I 

9930 SUBENO 

9940 SUB Subtract(K) 

995@ COM /Caly/ A(1@,10) ,C( 10) ,B(S).N Iprnt ,Opo Illog,Ifn,Ijyoin,Njyoin 
9960 FOR I=k TON 

937@ FOR Jek-1 TON 

9980 ACI J )@A(K=-1,J)-ACI J) 

93890 NEXT J 

1@Q@@@ Ci 1 d=C(K=-1 )-C(T) 

10010 NEXT I 

1002@ SUBENO 

10030 SUB Plin 

10040 COM /Caly/ AC1@,10),C(10) BCS) N,Iprnt ,Opo Ilog,Ifn,Ijoin,Njoin 
10@S@ COM /Xxyy/ Xx(25) ,Yy(25) 

1@@6@ PRINTER IS 785 


1a sae 


19070 BEEP 

1@@8@ INPUT “WANT TO PLOT Uo vs Vw? (17 28N)” [uo 
10092 IF Iuos@ THEN 

10100 BEEP 

1Q@1t@ INPUT “SELECT (Q=#n/nQ@z same tubde,!en( HF )/n(sm)* ,Irt 
10120 BEEP 

1013@ INPUT “SELECT h/n RATIO (I*FILE ,O@=sCOMPUTED)® ,Inrat 
10142 IF Inrat#@ THEN 

101S®@ BEEP 

10160 INPUT “WHICH Tsat (196.7 ,@5-2.2)° ,Isat 
10178 ENO IF 

12180 Xmined 

10192 Xmax=#1Q 

10220 Xsteps2 

10210 IF Irt#®@ THEN 

10228 YminsQ 

10232 Ymaxel.4 

10242 Ysteose.2 

1@25@ ELSE 

10262 Ymin=@ 

10272 Ymaxel!lS 

10282 YstepsS 

189298 ENO IF 

10@3@0 ELSE 

10318 Opo#2 

10328 Ymin=Q 

1@330 Ymaxe12 

109342 Ysteps3 

10352 Xmin=@ 

10362 Xmaxe4 

1Q@372 Xstepe!| 

1@38@ ENO IF 

10390 IF Ihrat#! THEN 

12422 Ymin=Q 

10410 Ymax#18 - 

10428 Ystepe3 

12432 xXmins@ 

1Q@442 Xmax#9 

12452 Xstep=2 

1046@ ENO IF 

10472 BEEP 

12482 PRINT “IN: SPIsIP 2302 ,2200 ,8302 68200; - 
10492 PRINT °SC @,100,@,1@@;TL 2,@:° 
10500 Sfx=#1Q00/(Xmax-xXmin) 

1QS1@ SfystQQ/( Ymax-Ymin) 

10520 PRINT “PU @,@ PO° 

1@S53@ FOR XaeXmin TO Xmax STEP Xstep 
10540 Xe(Xa-Xmin)eSfx 

10550 PRINT “PA°:X," Q@4 XTi° 

19562 NEXT Xa 

1Q@S7Q@ PRINT “PA 102@,0:PU;s° 

10@58@ PRINT “PU PA @,2@ PO* 

10592 FOR YawYmin TO Ymax STEP Yasatep 
19602 Ye( Ya-YmindeSfy 

1Q@61@ PRINT “PA @,°:Y,°YT" 

10620 NEXT Ya 

1Q@63@ PRINT “PA @,10@ TL @ 2° 

19640 FOR Xaexmin TO Xmax STEP Xstep 
1@6SQ@ Xs(Xa-Xmin)eSfx 

10660 PRINT °“PA°;X "1003 XT° 

1Q@673 NEXT Xa 

1@668@ PRINT “PA 120,100 PU PA 180,2 PO° 
109699 FOR YaeYmin TO Ymax STEP Ystep 
10700 Yel Ya-YmindeSfy 

1071@ PRINT “PO PA 10@,°,Y,°YT° 
1072@ NEXT Ya 


ZG 


1Q732 
10740 
12752 
10762 
10772 
10782 
107932 
12802 
1281 
12822 
10830 
12842 
1@BS@ 
1086 
10870 
12882 
18892 
10922 
10912 
12922 
12932 
12942 
12352 
12962 
1097 
12982 
10992 
1122 
11812 


11@2@ 
11032 
11@4@ 
11852 
11 @6@ 
1107 
11082 
11892 
111@@ 
1111 
1112 
1113Q) 
11142 
11152 
1116@! 
11172 
1118 
1119 
112080 
11218 
Ilkez @ 
1123 
1124Q 
11252 
1126 
11272 
112B@ 
11292 
1132@ 
11312 
11320 
11330 
11342 
11352 
11362 
11372 


PRINT "PA 100,102 PU 

PRINT “PA @,-c 35h 1.5,2° 

FOR XaeXmin TO Xmax STEP Xsteo 
Xe(Xa-Xmin)eSfx 

PRINT ‘PA°GX,” ,Q@4° 

IF [uo#@ THEN PRINT “CP -2,-l11LB°iXa;s"° 
IF Iuo=!l THEN PRINT “CP -1.5,-!1LB°:Xai~° 
NEXT Xa 

PRINT “PU PA @,@° 

FOR YaeYmin TO Ymax STEP Ysteo 

IF ABS(Ya)<1.E-S THEN Yas 
Ye(Ya-Ymin)*#Sfy 

PRINT “PA @,°3Y,°° 

IF Iuo#@ THEN PRINT “CP -4,-.2S5:LB"1Ya: °° 
IF Iuoe#l THEN PRINT “CP -3,-.2S:LB":sYas°” 
NEXT Ya 

Xlabal$=#"O1l Percent’ 

IF Tuo#@ THEN 

IF Irt*#@ THEN 


Ylabel9="h/hQz° 

ELSE 

YlabelS="h/hsmooth’ 

ENO IF 

PRINT °SR 1.5,2:PU PA SQ@,-10 CP" :-LEN( Xlabal$)/2;:°Q@:LB° s;Xlabals;°° 
PRINT “PA -11,50 CP @,°;-LEN(Ylabel$)/2¢S/6;°OI O,1;LB°s;YlebelS:"* 
PRINT “CP @,0° 

ELSE 


PRINT °SPQ:SP2° 
PRINT “SR 1.2,2.4:PU PA -8,35;01 @,1sLBUsPR 1 ,8.S:LBorPR -1 ,0.5:LB (kW/m 


PRINT "PR -! ,@.5:5R 1,1.5:L82;SR 1.5,2:PR .5,.5:L8.:PR .S,0@:L8K)° 
PRINT “PA 42,-10:01 1 ,O@;LBViPR .4,-13;LBwsPR 1,.5:LB( m/s)" 
PRINT °SPQ;5P1° 
ENO IF 

lon=2 
BEEP 

INPUT “WANT TO PLOT DATA FROM A FILE (1#Y ,Q"N)?° Okp 
Ien=@ 

IF Okp*#! THEN 
BEEP 

INPUT “ENTER THE NAME OF THE DATA FILE” 0 files 

IF Tuoe®@ THEN 
BEEP 

INPUT “SELECT (Q*LINEAR, !#LO6(X ,Y)* ,Ilog 

ENO IF 
ASSIGN O@File TO O_file$ 
BEEP 

INPUT “ENTER THE BEGINNING RUN NUMBER® Md 
BEEP 

INPUT “ENTER THE NUMBER OF X-Y PAIRS STORED’ Npairs 
IF ITuo#@ AND Imrate@ THEN 


BEEP 

cNP UTS ENTER OESITRED HEAT FLUX" (Q 
ENO IF 

BEEF 

PRINTER IS | 


PRINT USING °4X ,"°Select a symbol:°°° 

PRINT USING “4X ,°"! Star 2 Plus sitgn’’” 
PRINT USING “4X ,°°3 Circle 4 Sauare**° 

PRINT USING *4X,°°S Romous”** 

PRINT USING *4X,°°6& Rignht-side-up triangle’*” 
PRINT USING °4X,°°7 Up-side-down triangla’** 
INPUT Sym 

PRINTER IS 7@5 

PRINT “PU DOI’ 

IF Sym=t THEN PRINT “SMe° 


Lay 


11380 IF Syme2 THEN PRINT “SM+" 

11390 IF Sym=*3 TREN PRIN: 3Mo 

11400 Nn=#4 

1141@ IF Ilog*! THEN Nn=! 

11420 IF Md>1 THEN 

1143@ FOR I=#!1 TO (Ma-1) 

1144@ ENTER OFileixa,Ya 

1145@ NEXT I 

11462 ENO IF 

11472 IF Inmrat#@ THEN 

11482 Q1=Q 

11492 IF Ilog=#! THEN Q=L06(Q) 

11$@@ END IF 

11518 FOR I#!1 TO Npairs 

1192@ IF Tuo#@ AND Inrat#@ THEN 

11530 ENTER @FilesXa Ble) 

11542 Ya=B( 2) 

1155@ FOR Kel TO Nn 

11562 YasYar+Bik eQ*K 

11578 NEXT K 

11§8@ ENO IF 

11$9@ IF Iuo#l OR Imrat#!1 THEN 

1160@ ENTER ®File:Xa,Ya 

11618 IF Iuo#l THEN Ya=#Ya/128202 

1162@ END IF 

11638 IF Iuo#@ AND Imrat#@ THEN 

11648 IF Illog=! THEN Ya=EXP( Ya) 

11652 IF Tlog#=@ THEN Ya=Q1/Ya 

11662 IF Irt=*®@ THEN 

11672 IF Xea#@ THEN 

11688 YorYa 

11692 Ya=! 

11700 ELSE 

11718 Yaesyra/Yo 

1172@ ENO IF 

1173@ ELSE 

11742 Ham=FNHsmootn(Q Xa Isat) 

11758 YasYa/Hsm 

11760 END IF 

11770 END IF 

11782 Xx(I-1 exe 

11799 Yy(I-1)#Ya 

11800 X=(Ka-XmindeSfx 

11810 Ye( Ya-Ymin)eSfy 

11820 IF Sym>3 THEN PRINT °SM° 

11830 IF Sym<4 THEN PRINT °SR 1.4,2.4° 

11842 PRINT °PA* X,Y,"° 

1185@ IF Sym>3 THEN PRINT °SR 1.2,1.6° 

11862 IF Sym=4 THEN PRINT °UC2 ,4,99,8,-8,-4,0,0,8,4,20,1° 
11870 IF SymeS THEN PRINT °UC3 ,0,99,-3.7°6 .-3,6,3,6.3,763~- 
11880 IF Sym=6& THEN PRINT °UC@,5.3,99,3,-8,-6 .2,3,81° 
1189@ IF Sym=7 THEN PRINT °UC@,-S.3,99,-3,8,6 .0@,-3,-8:1" 
1198@ NEXT I 

1191@ BEEP 

11920 ASSIGN ®File TO « 

11930 END IF 

11948 PRINT “PU SM° 

1195@ BEEP 

11968 INPUT °WANT TO PLOT A POLYNOMIAL (1t#Y¥ @92N)?° Okp 
11970 IF Okp=1 THEN 

11982 BEEP 

11999 PRINTER IS 1 

12@@Q PRINT USING °4X ,.""Select lina type:°°° 
12018 PRINT USING °6X,°°2 Se lide line: 
12220 PRINT USING “6X,°"1 Dasnec’** 
12030 PRINT USING °6X,°°2,,,5 Longar lina - dasn°*°* 


TZ 


12042 
l<e02a0 
122862 
1287 
12882 
12092 
1218@ 
12110 
12122 
12132 
12148 
12158 
12162 
12170 
12182 
12192 
1228 
12212 
e220 
ee 30 
12242 
12252 
1226 
lzz70 
12280 
1229 
12380 
12312 
12320 
2430 
12348 
1235 
12362 
12370 
12382 
Vzag0 
12422 
12412 
12422 
12430 
12448 
12452 
12460 
12472 
12482 
12492 
12582 
12518 
12520 
12530 
12548 
12552 
12562 
1257 
12582 
12890 
12622 
12612 
12622 
12630 
12640 
12652 
1266@ 
12672 
12682 
1269 


INPUT Ipn 

rFrINTem [La 705 

BEEP 

INPUT “SELECT (QeLINEAR,1*LO06(x% ,Y))* ,Ilog 
Iprnte!} 

CALL Poly 

IF Iuo#l THEN 

BEEP 

INPUT “DESIRE TO SET X Lower and Upper Limit (10Y ,Q0N)?° Ixlim 
IF Ixlime@ THEN 

X11#2 

Xul=7 

ENO IF 

IF Ixlimel THEN 

BEEP 

INPUT “ENTER X Lower Limit’ X11 
BEE 

INPUT “ENTER X Upper Limit” ,Xul 
ENO IF 

ENO IF 

FOR XeeXll TO Xul STEP Xstep/2S 
Ienelcnel 

YaeFNPoly(Xe) 

IF Iuoel THEN YeeYa/ 1200 
Ye(YerYmindeSfy 

X=a(Xe-xmin)eSfx 

IF Y<@ THEN Y#2 

IF Y>10@ THEN GOTO 12412 

PueQ 

IF Ipmel THEN IdfelIcn MOO 2 

IF Ipne2 THEN Idfelcen MOO 4 

IF Ipne3 THEN Idfelcn MOO 8 

IF Ipne4 THEN Idfelcnm MOD 15 

IF Ion=S THEN Idfelen MOO 32 

IF Idf@l THEN uel 

IF Pue® THEN PRINT “PAS X,Y ,°PO° 
IF Puel THEN PRINT “PA* X Y,°PU" 
NEXT Xe 

PRINT “PU” 

60TO 110708 

ENO IF 

SEEP 

INPUT "WANT TO QUIT (1#Y Q°N)?° Iquit 
IF Iquitel THEN 12492 


60TO 11872 

PRINT “PU SPQ" 

SUBENO 

SUB Stets 

PRINTER IS 721 

J=Q 

K=aQ 

BEEP 

IF Iplote#l THEN ASSIGN @File TO P_file$s 
SEEP 


INPUT “LAST RUN No?(@*QUIT)* Nn 
IF Nn#@ THEN 12952 
NneNn-J 

S$x-2 

Sy*@ 

$z22 

Sxs2Q 

Sy s2Q 

$zse2 

FOR Iel TO Nn 
JeS+| 

ENTER @F116:0,T 


LEAS, 


12702 
12712 
12728 
1273 
12742 
12785 
12762 
lev 72 
12788 
1279 
12802 
12810 
12820 
12830 
12848 
12852 
1286 
12872 
12882 
1289 
1292 
123912 


12922 
1293 


RheQ/7 

Sx@Sx+*Q 

Sxs#Sxs+Q°2 

Sy#SyrT 

Sys*Sys*T°2 

$z2*52+H 

SzseSzs+H*2 

NEXT I 

QavesSx/Nn 

Taves#Sy/Nn 

HavesSz/Nn 

SdevaeSQR( ABS( (NneSxs—-Sx72 )/(NneCNnA1)))) 
Sdevt#SQR( ABS! (NneSys-Sy~2)/(Nne(Nnal )))) 
Sdevn=SQR( ABS! (Nn®$28-S2°2 )/(Nme@(Nnml )))) 
Sr#!1Q0*Sdevn/Have 

Sa*!Q0*Sdeva/Qave 

St#10Q*Sdevt/Tave 

IF Kel THEN 12932 

PRINT 

PRINT USING “11X,"“"°OATA FILE: °° ,14A° Files 
PRINT 

PRINT USING “°11X,° "RUN Htube SdeveH Qdp SdevQ Tmetab SdevT”” 


Ke | 
PRINT USING °11X 00 ,2(2x .0.30E ,1X ,30.20) 2X 00.30 1X (30.20 iJ (have saeeee 


,9Q,) ave ,St 


12942 
1295 
12962 
123972 
12982 
12992 
13222 
13212 
13820 
13032 
13242 
13252 
13260 
13072 
13282 
132390 
13182 
13118 
1312e 
13132 
13142 
13152 
13162 
13172 
13182 
13192 
13202 
13212 
Iie ey 
13230 
13242 
13288 
13262 
13278 
13282 
13259 
13322 
13312 
13322 
13332 


6O0TO 12572 

ASSIGN @Filel TO » 

PRINTER IS 1 

SUBENO 

SUB Coef 

COM /Caly/ AC1@,10) C(10),8(S) .N Iornt Ose Ilog lfm ljormenjorm 
BEEP 

INPUT “GIVE A NAME FOR CROSS-PLOT FILE® ,Cofs 
CREATE 8DAT Cof$,2 

ASSIE6N @File TO Cofs 

BEEP 

INPUT “SELECT (Q@=*LINEAR ,}#LO06(X ,Y))° Llog 
BEEP 

INPUT “ENTER OIL PERCENT (-1*S7T0P)° 80p 

IF Bop<®@ THEN 13!22 

CALL Poly 

OUTPUT OFilesBop ,8e) 

60TO 13062 

ASSI6N @File TO ° 

SUBENO 

SUB Wilson(Cf Ci > 

COM /Wii7 ©2701 06.76 Lu .Keu 

OIM Emf( 12) 

| WLISON PLOT SUBROUTINE OETERMINE CF AND Cl 
BEEP 

INPUT “ENTER DATA FILE NAME” Files 

SEEP 

PRINTER IS 1 

PRINT USING “4X ,““Select ontion:*”” 

PRINT USING °4xX,.°° Q@ Vary Cf ana C:i°°° 
PRINT USING °4xX,°° 1 Fix Cf Vary Ci°°° 
PRINT USING. °4x%.°° 2 Vasey Cl Fix Ci! 
INPUT “ENTER OPTION” ,Icfiax 

PRINTER IS 701 

IF Iefix#@ THEN 13328 

IF Icfix>@ THEN BEEP 

IF Icfix#!l THEN INPUT “ENTER Cf° Csf 

IF lefix#2 THEN INPUT “ENTER CI° Ci 

PRINTER IS 1 

INPUT “Want To Vary Coeff?(l#yY ,@=N)° Iccoef 


P30 


13340 IF Iccoef#l "HEN INPUT “ENTER COEFF’ |R 
13358 PRINTER IS 721 

(3360 IF Icfixe@ OR Icfixe2 THEN Cfa=.004 
13370 IF Icfixel THEN CfaeCsf 

13382 CiaeCi 

13392 Xne#.8 

13400 Fre.3 

13410 J 20 

13420 Rre3. 

13432 IF Iccoafe!l THEN RreR 

1344@ PRINTER IS 1 

13450 PRINT Oo ,O01 ,Kcu 

13460 ASSIGN @File TO Files 

13470 ENTER @FilesNrun Dates ,_dtc! ,Ldtc2 ,Itt 
13480 RweOo0eL06( 00/01 )/(2*Kcu) 

13492 Sx#@ 

13500 Sy=#2 

13510 Sxy= 

13520 Sx2=@ 

13530 Sy2s@ 

13548 FOR Tel TO Nrun 

13550 ENTER @Filei8oe ,Told$ ,Emf(*) Fms 
13562! CONVERT EMF'S TO TEMPERATURE 
13870 FOR Jel TOS 

1358@ T(J @FNTvsviEmf(J)) 

13590 NEXT J 

136@Q@ Tsate(T(1)4#T(2))¢.§ 

1361@ TavgeT(S) 

13620 Grade37.9853+.104388eTavg 

13630 ToropeEmf(7)e!1 .E+6/(10.¢Grad) 
13648 Tavgceel(S)-Tdrope.S 

13650 IF A8S{ Tavg-Tavge )>.@!1 THEN 
13662 Tavge(Tavgt+tTavge )e.5 

13670 GOTO 13622 

1368@ ENO IF . 

13690! 

13700! Compute properties of water 
13712 KweFNKu( Tavg) 

13720 MuwaeFNMuwl Tavg) 

13732 Cow#FNCpow(\ Tavg) 

13740 PrweFNPrw(Tavg) 

13752 Rnow*FNRhow(! Tavg ) 

13760! 

13770! Comoute properties of Freon-114 
13780 LmtaeTdrop/LO6( (TiS )-Taat)/(T(S)-Tdrop-Tsat )) 
1379@ IF J )2@ THEN 

13800 TweTsat+FreLata 

13818 ThetabeTw-Tsat 

13820 Jj!) 

13832 ENO IF 

13840 Tfe(TweTsatde.S 

13850 Rno#FNRho( Tf) 

13860 Mue#FNMu( Tf ) 

13870 KeFNK( TF) 

1388@ Co=FNCe(TF ) 

1389Q@ BetasFNBeta( Tf) 

13980 Hf g*FNHfg( Tsat ) 

13918 Nt#Mu/Rho 

13920 Alpna#K/(RhoeCe ) 

13932 PreNi/Alpha 

1394Q! 

13950! Analysis 

1396@1 COMPUTE MOOT 

13970 AsP1e(00°%2-01°2)/4 

13980 PePIeDo 

1399@ Mdot#3.96S7E-3+Fmse(3.619SSE-3-Fmse(8.82QQ06E-6-Fmsae(1.23688E-7-Fmse4.31897 


IEehee 


E=119))) 

14002 Q=Mdot*Cow*lorop 

14210) COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
14@2@' FOR UNENHANCED END(S) 

14030 Hpare=1920 

14040 Fee(HbareP/(KecueA) )* .Selu 

140950 Tann=FNTanniFe) 

14262 Theta#TnetaneTann/Fe 

14070 Xx#(93.8leBetaeThetan*eDo 3ZeTann/(FeeNieAlona))*. 166867 
14080 Yye(1e(.SS9/Pr )*(9/16))°( 8/27) 

14092 HbarcekK/Doe!( .6+.387¢Xx/Yy)*2 

14100 IF ASS((Hbar-Hbarc )/Hbar )>.2281 THEN 

14112 Hbare(HoareHbarc de. 

14128 GOTO 14040 

1413@ ENO IF 

1414Q! 

14150! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENDS 
14168 Qle(HbarePeKkcueA)* .SeTnetaoeTann 

14170 Qe#Q-2eQt 

1418@ As=PTeO2eL 

14198! COMPUTE ACTUAL HEAT FLUX 

14200 QdoeQc/As 

14210 IF Iefixe@ OR Icfix>! THEN Csfel/Cf*(1./Rr) 
14220 Thetad=Csf/Coehfge( Qdo/( Muth )e( .814/(9.81eRno))*.S)*(1/Rr ePr 71.7 
14238 HosQdo/Thetab 

14240 OmegashHo/Cf 

14258 UosQ/(PIeDoeLeLatad) 

14260 VweMdot/(RnoweP1eD: °Z/4) 

14278 ReweRnoweVweD1/Muwae 

14288 TwieTwrQeRw/(PleDoel ) 

14290 Gamaeku/DieRew* .8ePrw (1/3. 04 Muwa/FNMuw( Twi) )*.14 
14300! PRINTER IS | 

14310 Ywe 1. /Uo-Rw )*Omega 

14320 XweDmega*Do/( GamaeD: ) 

14330 Sx#Sx+Xw 

143428 SysSy*Yu 

14352 SxyeSxytYweXw 

14360 Sx2eSx24XweXy 

14370 Sy2#Sy2*YueYw 

14382 NEXT I 

14392 ASSIEN @File TO 

14400 Me(SxeSy-NeuneSxy )/( Sx e¢Sx-NruneSx2 ) 

14418 Ce(Sy-SxeM)/Neun 

144220 IF Icfixe@ OR Icfixe3 OR Iefixed THEN 

14432 Cicel/M 

14440 Cfcesl/C 

14450 ENO IF 

14460 IF lIefixel THEN 

14470 Cricel/m™ 

14482 CfceCf 

14490 ENO IF 

14502 IF Icfixe2 THEN 

14510 CiceCi 

re520 Cfest7C 

1453@ ENO IF 

14548 IF ABS((Ci-Cic)/Ci¢e)>.001 DR ABS((CF-Cfe)/Cfe)>.@01 THEN 
14550 Cr=(Ci+Cicde.§ 

14560 Cfe(Cf+Cfere.S 

14572 PRINTER IS 1 

14580 PRINT USING "2X.°° Csf = “° M2.3DE,2%.°° Ci = °° (AZ. 2De wesc. 
14592 PRINTER IS 7 

14682 60TO 13462 


1461@ END IF 

14620 PRINT 

14630 PRINTER IS 701 

1464@ PRINT USING "23X,°° Cf Cras 


iS 2 


14650 
14662 
14672 
14682 
14692 
14722 
14712 
14720 
14730 
14742 
14752 
147628 
14772 
14788 
14792 
14820 
14812 
14820 
14830 
14840 
14850 
1486 
14878 
14880 
14892 
14980 
14912 
14928 
14932 
14940 
14952 
14962 
14972 
14982 
14990 
15082 
15812 
15822 
15032 
15042 
1585 
1506 
1587 
15282 
158392 
15180 
15112 
19120 
15132 
1514 
1515 
1516 
1517 
15182 
1519 
15222 
1521 
15220 
19230 
1524 
15252 
15268 
15272 
15280 
15232 
15302 


PRINT USING “8X ,~ "ASSUMED PME eSDe Ieee oe ota. 14 


PRINT USING °8x ,“"CALCULATED PML 2 ODESSA Mew ouUe @os't 253 
PRINT 

Sum2eSy2-2eMeSxy-2eCeSyomM* 2eSx2+Z2emeCeSx+Nerunel*2 
PRINT USING *18X,°°Sum of Squeres © °° ,2.30E°;Sum2 
PRINT USING “1OX ,°*Coeffictent = °° .0.000";Rr 
SUBEND 

DEF FNMuw( T) 

A@247.8/(7T+133.15) 

Mue2.4E-Se1Q*A 

RETURN Mu 

FNEND 

DEF FNCow(T ) 
Cowed.211208S8-Te( 2, 26826E-3-Te( 4. 42361E-S*2.71428E-7eT)) 
RETURN Cowe!Q20 

FNEND 

DEF FNRhow( T ) 

Ro#999.52946+Te( .01269-Te(S$.482513E-3-Te1l,. 234147E-S)) 
RETURN Ro 

FNEND 

DEF FNPrw(T ) 

PrweFNCow( T)ePNMuw( T)/FNKw(T ) 

RETURN Prw 

FNEND 

DEF FNKuw(T ) 

Xe(7+273.159/273.15 
Kwe-.92247+xe(2,8395-xXe(1.8087-Xe( .$2577-.87344ex ))) 
RETURN Kw 

FNEND 

SUB Plot 

Com /Coly7 AC 10 ,19)26010),8(S) Nop .lornt ,Oo0 ,llog,ifn,ijoin,Njoin 
OIM Bsi3) 

INTEGER 11 

PRINTER IS 1 

Idv=Q 

8EEP 

INPUT “LIKE DEFAULT VALUES FOR PLOT (1#Y ,Q@eN)7° Idv 
Opo=0 

8EEP 

PRINT USING °4X ,""Select Option:°*”*’ 

PRINT USING "6X,°°@ gq versus celta-T**" 

PRINT USING “6X ,"°t Wh versus delta-T*"*’ 

PRINT USING °“6X,°°2 bh versus q**” 

INPUT Opo 

8EEP 

INPUT “SELECT UNITS (@#SI ,1*ENGLISH)* ,Iun 

PRINTER IS 725 

IF Idv<>! THEN 

BEEP 

INPUT “ENTER NUMBER OF CYCLES FOR X-AXIS" Cx 

SEEP 

INPUT “ENTER NUMBER OF CYCLES FOR Y-AXIS’ Cy 

Beer 

INPUT “ENTER MIN X*VALUE (MULTIPLE OF 10)" .Xmin 

SEEP 

INPUT “ENTER MIN Y-VALUE (MULTIPLE OF 18)" ,Ymin 

EESe 

IF Opo#@ THEN 

Cys3 

Cxo3 

Xmine.] 

Ymin=1QQ 

END IF 

IF Opoe! THEN 

Cy#2 

Cxe2 


i3e 


1$3'2 
15320 
15330 
153428 
15252 
15362 
15372 
15380 
153928 
15488 
15412 
15422 
15430 
15442 
15450 
15462 
15478 
15480 
15492 
15500 
15510 
1552 
15530 
1554 
15550 
1556 
15570 
15580 
1559 
15682 
15612 
15622 
15632 
15648 
15658 
1S662 
15872 
15680 
15690 
15702 
15712 
15720 
15730 
15740 
15752 
15762 
15770 
1S782 
15792 
15822 
1$812 
15820 
15830 
15840 
15850 
15862 
15870 
15880 
15890 
15982 
1$912 
15922 
15932 
15940 
15950 
15962 


xmine.l 

raine1@e2 

ENO IF 

IF Opoe2 THEN 

IF Iune@ THEN 

Cys3 

Cx#2 

Xmin#1882 

Ymin#=122 

ELSE 

Cys3 

Cx=3 

Xmin#=102 

Yminel®@ 

ENO IF 

ENO IF 

ENO IF 

BEEP 

PRINT “IN:SPI; IP 2308 ,2200 ,8300 68001" 
PRINT °SC @,100,0,100:TL 2,0;° 
Sfxe1Q@0/Cx 

Sfy#l@Q/Cy 

BEEP 

INPUT “WANT TO BY-PASS CAGE? (18Y @9N)” Ibyp 
IF Ibyp2#! THEN 16792 

PRINT °PU @,@ PO* 

Nne#Q 

FOR Iel TO Cx! 
KatexminelO@*(I-1 ) 

IF ITeCx+!1 THEN Nnel 

FOR Jel TO Nn 

IF Jel THEN PRINT °TL 2 0° 
IF Je2 THEN PRINT °TL 1 O° 


XKaexate!s : 
X@eLET(Xa/Xmin)eSfx 
PRINT °“PA°IX,* Os XTi°® 
NEXT J 

NEXT I 

PRINT “PA 180 ,@1PUi1° 


PRINT “PU PA @,@ PO0° 

NneJ 

FOR I=! TO Cy! 
YateYminel@*(I-1 ) 

IF IeCy+! THEN Nnel 

FOR Je! TO Nn 

IF Jel THEN PRINT “TL 2 Q° 
IF J=2 THEN PRINT °TL I @° 


YaeYateJ 
YeL6T(Ya/Ymin)deSfy 
PRINT "PA Q@,°1Y,°YT° 
NEXT J 

NEXT I 

PRINT °PA 0,100 TL @ 2° 
Nne9 

FOR Iel TO Cx+! 


Katexminel@*( 1-1) 

IF IeCx+! THEN Nn#=! 

FOR Jel TO Nn 

IF Jel THEN PRINT “TL @ 2° 
IF J>1 THEN PRINT “TL @ 1° 
KaexateJ 


K@LBT( Xa/XmindeSfx 

PRINT “PA":X.° 1003 XT° 

NEXT J 

NEXT I 

PRINT °PA 188,108 PU PA 108,08 PO” 


1§9720 Nne3 

1§98@ FOR ie! TO Cyrl 

16990 Yat#YainelO*(I-1) 

16@@@ IF I#Cy+!1 THEN Nn®#! 

16@1@ FOR J#1 TO Nn 

16020 IF J#l TREN PRINT “TL @ 2° 

16@3@ IF J>1 TREN PRINT “TL @ 1° 

16048 YaeYateJ 

16050 YeLGT( Ya/YmindeSfy 

16@6@ PRINT “PO PA 1080,°,Y,°YT° 

16070 NEXT J 

16@8@ NEXT I 

16092 PRINT “PA 108,188 PU° 

165100 PRINT “PA @,-2 SR 1.S,2° 

16410 LislGT(xmin) 

16128 FOR I#1 TO Cx+1 

16130 XawXminelO*(I-1 ) 

16140 X=@L6T(Xa/Xmin)eSfx 

16150 PRINT “PA‘:X .* 83° 

16160 IF Ii1>#@ THEN PRINT °CP -2,-2:LBIQ@:;PR -2 ,2:L8°s125°° 
16178 IF 11<8@ THEN PRINT “CP -2,-2;:LBIQ@sPR O@,2:LB°s1i3°* 
16180 TislIi+l 

16190 NEXT I 

16200 PRINT “PU PA @,2° 

16210 TistL GTC Yman) 

16220 Y1Q0#12 

1623@ FOR I=! TO Cy+! 

16248 YaaYminel@*(I-1) 

16250 YeLGT( Ya/YmindeSfy 

16260 PRINT “PA @,°:3Y,°° 

16278 PRINT “CP -4,-.2S5:LB10:PR -2,2:LB°sI23°° 

16288 Tieli+l 

16298 NEXT I 

16300 BEEP 

16310 INPUT “WANT USE DEFAULT LABELS (IY ,@2N)?° Id] 
16320 IF Idl<>1 THEN 

16338 BEEP 

16340 INPUT “ENTER X-LABEL” xXlabel$ 

16358 BEEP 

16362 INPUT "ENTER Y-LABEL” Ylavel$ 

16370 ENO IF 

16388 IF Opo<2 THEN 

16390 PRINT “SR 1,2:PU PA 40 ,-14;° 

16400 PRINT “LO(TsPR -1.6,3 PO PR 1.2,0 PUIPR .S,-4:LBwosPR .S,1;° 
16410 PRINT “LO-TsPR .S,-1lslLBsatsPR .S.13° 

16420 IF Iun#@ THEN 

16430 PRINT *LB) (K)° 

164428 ELSE 

1645@ PRINT °LB) (F)* 

16468 END IF 

16478 ENO IF 

16488 IF Opo#2 THEN 

16490 IF Tun#=@ THEN 

15500 PRINT “SR 1.5 ,2:PU PA 40 ,-141L6q (W/m;SR 1,1.5:PR @.5,1:L62:SR 1.8,2:PR 
@.S,-1;:L8)° 

16510 ELSE 

MescOer einige sk elasecgieu rw. o4 —leé,lGq (6tu/mcrPR .S,..S+tL6.3PR .S,—.S9" 
TeSS@ereRins LettteR 1.5.1 95R 1.1.5sL62:S5R 1.5,2:PR .S.-1:5L6)1- 
16540 END IF 

16SS@ ENO IF 

16560 IF Dpo#@ TREN 

16870 IF ITun#@ TREN 


16580 PRINT °SR 1.5,2:PU PA -12 40:01 @,\)1LBq (W/msPR -1 0.95:S5R 1,1.5:L62:5R 1 
.§,21PR 1,.5:3L8)" 
16590 ELSE 


[soon sk lso,cteu PA —12.32401 @,13LBa (6tu/mhriPR -.S,.StL6.9PR .5..51° 


ers 


166 '@ 
16622 
16630 
1664 
166S@ 


ui 
Lea) 
ay 


PRINT “L&ft;:SR 1.1.5:PR =o .5Sncee-en 
END IF 

ENO IF 

IF Op0>@ THEN 

IF Iun=@ THEN 


+. oS ees 


16660! PRINT °*SR 1.5,2:PU PA -12 38:01 O,1:L Br (W/msPR -1,.5:5R 1,1.5:L82;55R | 
<5 ,cIPR .S..5t" 


16670 
16682 
1669@ 
187280 
1671@ 
16722 
1673 
16742 
16752 
1876@ 
16772 
16782 
1679 
16822 
1681 
16822 
1683 
16840 
1685@ 
16862 
16872 
16882 
16852 
16902 
16391 
16922 
1633@ 
16942 
163952 
16982! 
163972! 
16382 
16992 
17000 
17212 
1702@ 
17Q@3@ 
17Q@4@Q 
17852 
17082 
17872 
17Q82@ 
17092 
17182 
1711 
17128 
17132 
1714 
17152 
171682 
1717 
17182 
17192 
1722@ 
17212 
17220 
17230 
1724 


PRINT °SR 1.2,2.48PU PA -12,37:01 O,1:LBAsPR 1 ,@.5:LBo:PR -1 0.5118 (W/m 


PRINT "PR -1,.5:SR 1,1.5:L82:5R 1.5,2:PR .5,.5:3:LB8.:PR .5,Q@:L8K )° 
ELSE 

PRINT “SR 1.5,2:PU PA -12,28:01 @,1:LBn (8tu/mhriPR -.5, Si lB01PRe Seon 
PRINT “LBft:PR -1,.5:5SR ¢,1.5st82:1SR 1'.5,2:PR .5,.5:L8.:PR .5,.5:L5F )" 
ENO IF 

ENO IF 

IF Idl*s@ THEN 

PRINT "SR 1.5,2:PU PA S0,-16 CP°:-LEN(Xlabel$ )/2;°O@;L8° sXlabel$:i°” 
PRINT “PA -14,5@ CP @,°s-LEN(Ylapel$)/2¢S/65;°DI O@,1:L8°:Ylanel$: °° 
PRINT “CP @,@ DI° 

END IF 

Ipn=#2@ 

X}l#1.E+6 

Xule-1.6*§ 

IensQ@ 

Lfne#Q 

Tjyoinel 

BEEP 

INPUT “WANT TO PLOT DATA FROM A FILE (18° ,@8N)?° Dk 

IF Dkel THEN 

BEEP 

INPUT “ENTER THE NAME OF THE DATA FILE* ,DO_files 

ASSIGN ®@File TO O_file$s 


BEEP 
BEEP 
INPUT “ENTER THE BEGINNING RUN NUMBER® Ma 
BEEP ' 
INPUT “ENTER THE NUMBER DF X-Y PAIRS STDORED® ,Noairs 
BEEP 
INPUT "CONNECT DATA WITH LINE (18Y ,Q@0N)?7"° Ic] 
BEEP 
PRINTER IS 1 
PRINT USING *4X ,°°Select a symool:°*°*°* 
PRINT USING “8X ,°"! Star 2 Plus sign""* 
PRINT USING °68X,°°3 Circle 4 Square"** 


PRINT USING °“6X,°"°S Romous"*" 
PRINT USING °&X,°°& Right-side-up triangie’®* 
PRINT USING °S&X,°"? Up-side-down triangle’ °° 
INPUT Sym 

PRINTER IS 795 

PRINT °PU DI’ 

IF Symel THEN PRINT “SMe° 

IF Syme2 THEN PRINT “SM+° 

IF Sym*3 THEN PRINT °SMo° 

IF Md>! THEN 

FOR Ie! TO (Md-1) 

ENTER OFilesYa,Xa 

NEXT I 

END IF 

FOR Ie! TO Npairs 

ENTER OFilerYa,Xa 

IF I@!l THEN Qte#Ya 

IF I*Npairs THEN Q2*Ya 

IF Dopo#!l THEN Ya#Ya/Xa 

IF Dpo#2 THEN 

Q=Ya 

YaeYa/Xa 


ISHS 


172S5Q@ tae@ 

17260 €NO IF 

172708 IF Xa<Xll THEN Xl1#Xa 

17282 IF Xa>Xul THEN Xul#Xa 

17292 IF Iun#! THEN 

17302 IF Opo<2 THEN XaeXael.8 

17310 IF Qpo>@ THEN Yae#Yae.!76} 

17320 IF Opo#®@ THEN Yas#Yae.317 

17330 IF Opo#2 THEN Xa=Xae.317 

17348 ENO IF 

173SQ XeL6T( Xa/XmindeSfx 

17362 YeL6T( Ya/Ymin)eSfy 

17378 Kj#0 

1738@ CALL Symb(X ,Y ,Sym,Icl ,Ky) 

17392 60TO 17522 

17400 IF Sym>3 THEN PRINT °SM° 

17410 IF Sym<4 THEN PRINT “SR 1.4,2,4" 

17420 IF Icl#@ THEN 

1743@ PRINT “PA* X,\Y,°° 

17440 ELSE 

174S5Q@ PRINT “PA* XxX ,Y,°PQ" 

1746@ ENO IF 

17470 IF Sym>3 THEN PRINT “SR 1.2,1.6° 

17480 IF Sym=4 THEN PRINT °UC2,4,99,0,-8 -4,0,20,8,4 ,Q:° 
17490 IF Sym*S THEN PRINT °UC3,2@,99,-3,7-6,-3,6,3,6,3,-6:° 
17S5@@ IF Sym=6 THEN PRINT °UCO,S.3,99,3,-8,-6 ,@,3,8:° 
17510 IF Sym=7 THEN PRINT °UCO@ .-S.3,99,-3,8,6 ,@,-3 ,-8:’ 
17520 NEXT I 

17530 PRINT °PU* 

1754@ 8EEP 

175S@ INPUT “WANT TO LABEL? (18Y ,Q@8N)° ,Ilab 

1756@ IF Tlab=! THEN 

17572 PRINT °SPQ:SP2° 

1758@ 8EEP 

17530 IF Klab#@ THEN 

176Q@ XlabeS 

17612 Ylab#85S ; 

1762@ INPUT “ENTER INITIAL X,Y LOCATIONS” ,Xlab ,Ylab 
1763@ Xtt=Xlab-S 

17640 YtraVYlabes8 


17550 PRINT “SR 1,1.5° 

I7GBS PRINT “SMIPA™ ,.Xtt ,Ytt ,~LB r Heat File® 
l7o7O.Ytteartt=3 

17SBO (PRINT “PA® Xtt ,Y¥tt, LE O11 Flux Name” 


17690 IF Syme! THEN PRINT °SMe° 

17700 IF Sym*2 THEN PRINT °SM+° 

17718 IF Sym=3 THEN PRINT °SMo° 

17720 Klab@=! 

1773@ ENO IF 

17748 Kyl 

17750 CALL Symb(Xlab ,Ylab Sym,Icl ,Kj) 
17760 PRINT °SR 1,1.54:SM° 

17778 IF Sym<& THEN PRINT °PR 2,2° 

17782 BEEP 

17792 INPUT *ENTER BOP° Bop 

17800 IF Bop<!1@ THEN PRINT “PR 2,0:LB°1Boo:’* 
17818 IF Bop>9 THEN PRINT "PR .S ,Q@:LB°s8op:"" 
17828 Inf#Q 

17830 IF Q1!>Q2 THEN I[hf#} 

17842 IF Inf*O@ THEN PRINT “PR 4 Q:LBInc’ 
17850 IF Inf#! THEN PRINT “PR 4 0;:LBQec* 
17862 PRINT °PR 2,Q@:LB°:O0_file$i°° 

17878 PRINT °SP@:;SP1;5R 1.9 ,2° 

17882 YlabeYlab-S 

1789@ ENO IF 

179@@ BEEP 


esi 


1791 
17328 
1793 


1794! 


17982 
17962 
17972 
17982 
17992 
18822 
18212 
1802Q 
18232 
18042 
18052 
18262 
18070 
1828 
18292 
18182 
1811 
18122 
1813 
18142 
18152 
18162 
18172 
1818 
18192 
18222 
18212 
18220 
18232 
18240 
18252 
18262 
18272 
18282 
1829 
18322 
18310 
18320 
18332 
18340 
18352 
18360 
18372 
18382 
18392 
1842@ 
1841Q 
1842Q 
18432 
18442 
18452 
18462 
18472 
18482 
18492 
18500 
1851 
18528 
18530 
18542 
18552 
18562 


ASSION @File TO 
Aiee eel. 2 
Xulexulel.2 

60TO 8042 

ENGI 

PRINT “PU SM* 
SEEP 


INPUT “WANT TO PLOT A POLYNOMIAL (1#Y ,Q@2N)7° G0_on 
IF 60_on#*! THEN 

8EEP 

PRINTER IS | 

PRINT USING °4X,°"Select line type:*°* 

PRINT USING °6X,°°Q Solid line*** 

PRINT USING °6X,°°1 Dasned*** 

PRINT USING “6X ,°°2,,,S Longer line - dash°** 
INPUT Ipn 

PRINTER IS 72S 

SEEP 


INPUT “SELECT (Q@*LIN,1eL06(xX,Y))° ,Ilog 
lorntel 

CALL Poly 

IF Ifne@ THEN 

8EEP 

INPUT “ENTER NUMBER OF FILES TO JOIN?* ,Njoin 
ENO IF 

TjoinsQ 

IF Ifn<Njyoin THEN Ijoine!l 

IF Ifne@ OR Tjyoine! THEN 

FOR Ijy#@ TO 3 

B8s( ly reBst( ly +B ly) 

NEXT Ij 

Ifnelfnel 

ENO IF 

IF Njyoinelfn THEN 

FOR [)}#@ TO 3 

BC Ty @Bs(ly)/Nyoin 

Bs( Ij )=@ 

NEXT [yj 

lfn=Q 

ELSE 

60TO 16852 

ENO IF 

BEEP 

INPUT “ENTER Y LOWER ANDO UPPER LIMITS* ,Y11,Yul 
FOR Xxe@ TO Cx STEP Cx/220 


XaekKmine |Q°Xx 

IF Xa<X1l1 OR Xa>Xul THEN 18642 
lenelenst 

Pu=Q 


IF Ione! THEN Idfelcen MOO 2 

IF Ipn=2 THEN Idfelen MOD 4 

IF Ion=3 THEN Idfelen MOO 8 

IF Ipne4 THEN Idfe#elen MOO 16 

IF Ipn#S THEN Idfelen MOD 28 

IF Idf=! THEN Puel 

IF Doo#@ THEN Ya#FNPoly( Xa) 

IF Dpoe2 ANDO Iloge@ THEN Ya=#Xa/FNPoly( Xa) 
IF Opo#2 AND Ilog#! THEN YasFNPoly( Xa) 
IF Opo#l THEN YasFNPoly( Xa) 

IF Ya<Ymin THEN 18640 

IF Ya<Yll OR Ya>Yul THEN 18642 

IF Tune! THEN 

IF Opo<2 THEN XaeXae!l.8 

IF Do0o>@ THEN YawYae.176! 

IF Dpo#@ THEN Ya#Yae,.317 

IF Opo#2 THEN XaeXae.317 


3.6 


las7Q 
185388 
18598 
18602 
18618 
18620 
18632 
18642 
18652 
18662 
1867 
18682 
18698 
18782 
18712 
18722 
18732 
18748 
18752 
18768 
18772 
18782 
18792 
18880 
18812 
18828 
18830 
18842 
18858 
188628 
18878 
1888 
18892 
1890 
1891 
18922 
18930 
18942 
1895 
18962 
18972 
1898 
18992 
19222 
198128 
19820 
1983 
19842 
19852 
19862 
19878 
19288 
19892 
19102 
1911 
19122 
19132 
19142 
19158 
1916 
19178 
19182 
19198 
19282 
1921 
19228 


ENO IF 

Yecol(Ya/iminvesty 
KX@LGT(Xa/Xmin)eSfx 

IF Y¥<d THEN Y#® 

IF Y>10@ THEN GOTO 18642 

IF Pue®@ THEN PRINT “PA* X,Y ,°PO* 
IF Puel THEN PRINT “PA’ X,Y,"PU" 
NEXT Xx 

PRINT °PU° 

60TO 16858 

ENO IF 

SEEP 

INPUT “WANT TO PLOT REILLY'S OATA? (1#Y ,QeN)° |Irly 
IF Opo#@ OR Opo#! THEN 

Xlle3 

Xul#2Q 

ENO IF 

IF Opos2 THEN 

X11=10002@ 

Xul=100008 

ENO IF 

IF Irlye#l THEN 

Y¥11=22 

Yul=7@ 

SEEP 


INPUT “ENTER LOWER ANO UPPER Y-LIMITS FOR PLOTTING" ,Y11 ,Yul 


FOR Xx#@ TO Cx STEP Cx/20 
Xaexmine1O*Xx 

IF Xa<Xll OR Xa>Xul THEN 18982 
Xl*LO6( Xa) 


IF Opo#@ THEN Yle-2.5402837E-14+X19( 4.972015 1-xXle2.S134787E-1 ) 
IF Opoe) THEN Yle-2.9402837E-1+x1°(3.97201S1-xX1¢2.5134787E-1 ) 
IF Opoe2 THEN Yle-3.70738Q1E-1+xX1°(8.7259190E-1-x19°6.8826842E-3 ) 


YaeExP(yYl) 

YL6T( Ya/YmindeSfy 

KeLOT( Xa/Xmin)eSfx 

Ipus? 

IF Y¥<Y1l1 THEN Ipuel 

IF Y>Yul THEN 60TO 18982 

IF Ipue@ THEN PRINT “PA* XY ,"PO° 
IF Ipuel THEN PRINT “PA” xX Y,°PU" 
NEXT Xx 

PRINT °PU* 

ENO IF 

SEEP 


INPUT “WANT TO PLOT ROHSENOW CORRELATION? (1#Y ,Q=N)° 


IF Iromse#l THEN 

Ylle1S 

Yul#82 

SEEP 

INPUT “ENTER Taat (O4ag CC)” ,Tsat 
Csf=.00420 

BEEP 

INPUT “ENTER Csf (OEF#@.004)° Csf 
TfeTsat+2 

FOR Xx#Q@ TO Cx STEP Cx/200 
XaexminelO*xx 

IF xXa<Xll OR Xa>Xul THEN 19422 
XLeLO6( Xa) 

IF Opo<2 THEN Tf#Tsat+Xa/2 
Rho#*FNRno( Tf ) 

KeFPNK( Tf ) 

MusFNMu( TF ) 

Cp*FNCp( TF) 

Hf geFNHfg( Tsat ) 

Ni*Mu/Rho 


abe ie, 


Irons 


3232 
19242 
19252 
19268 
19272 
19282 
19292 
1932 
19312 
19320 
3550 
19342 
19352 
1936 
19372 
19382 
19332 
19422 
19412 
19422 
1943 
19442 
19452 
19462 
19472 
19482 
19492 
13520 
19512 
19920 
1953 
19542 
19552 
19856 
19572 
1958 
19532 
19682 
19612 
19622 
1363 
19642 
19652 
19662 
1967 
19682 
19692 
139722 
1971 
19722 
19732 
19742 
19752 
19762 
19772 
139782 
19792 
19880 
19B12 
19822 
19832 
19842 
19852 
19862 
139B72 
19882 


Pe=stoe aah 

OmegescstehfosCosit .@)4703.81 "kno. os staedee - & geet 
IF Opoed THEN Yae(xa/Omega?)*3 

IF Qpo#! THEN Yae( Xa/Omega)*3/Xa 
IF Opo#2 THEN Yae#xa*(2./3)/Omaga 
IF Opoe2 THEN 

TKcoeTsateXa/Yae,S 

IF ABS(TF-Tfe)>.Q@1 THEN 

TRa( TKeTFe)e.§ 

60TO 19172 

ENO IF 

ENO IF 

YeLGT(Ya/YmindeSfy 

X@L6T( XKa/XmindeSfx 

Ipued 

IF Y<Y11 THEN Ipue! 

IF YOYul THEN 19422 


IF Tpue@ THEN PRINT “PAS x Y ,*PO° 
IF Ipuel THEN PRINT “PAS x Y¥,°PU° 
NEXT Xx 

PRINT °PU° 

ENO IF 

BEEP 

INPUT “WANT TO QUIT (ley QeN)° ,Iqt 


IF Iqte!l THEN 19492 

S0TO 1688 

PRINT “PU PA @,2 SPQ° 

SUBEND 

SUB Symp( X,Y ,Sym,Icl ,Kj) 

IF Sym>3 THEN PRINT °SM“ 

IF Sym<4 THEN PRINT “SR 1.4,2.4° 

YadeQ 

IF Kyo! THEN Yads.8 

IF TeleQ@ THEN 

PRINT “PA* X .Y+Yad,*” 

ELSE 

PRINT “PAS .X YeYad,°PO° 

ENO IF 

IF Sym>3 THEN PRINT “SR 1.2,1.6° 

IF Syme4 THEN PRINT “UC2 4,99 .8,-8 ,-4,0,.0,8,4,283° 
IF Sym@S THEN PRINT °UC3 0,55 -3,-5-3.6. 57072 7con 
IF Syme& THEN PRINT °UC@,S.3,95,3,-38,-8,8,3,81" 
IF Syme? THEN PRINT °UC@.-S.3,99,-3,8.6,.8,-3,-8:° 
IF Kyel THEN PRINT “SMiPR @,-.8° 

SUBEND 

SUB Purg 

BEEP 

INPUT “ENTER FILE NAME TO BE OELETED" Files 

PURGE Files 

GOTO 196392 

SUBEND 

SUB Tden 

COM cect (7) ee 

OIM Emf(}) 

DATA @.10086091 ,29727.943639 ,-767345.829S5 ,78025S95.81 
OATA -9247486S89 6.97688E+11 .-2.66192EF+13 ,3.94078E+14 


READ Cie) 

BEEP 

INPUT “GIVE A NAME FOR FILE TO 8E CREATED" ,Files 
BEEP 

INPUT “SELECT TUBE (QeWH ,leHF 2eWT)° Itt 

BEEP 

INPUT “SELECT THERMOCOUPLE TYPE (@=NEWU 1 =020) sical 


IF Itt<2 THEN 018.0127 
CREATE BDAT Filed ,4 
ASSIEN @File TO Files 
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9692 
19988 
1991 
19928 
19938 
19948 
19952 
19968 
19978 
19988 
19990 
22000 
20012 
20220 
20830 
20042 
28050 
220062 
2202072 
28088 
228390 
22100 
20112 
228120 
22132 
22148 
20150 
22160 
20172 
22182 
2019 
28280 
28210 
28220 
20230 
28248 
20250 
20262 
28272 
22282 
20292 
20302 
220318 
22320 
28330 
20342 
220350 
20362 
20378 
20382 


GUTRU Mer ier .tt 

J=2 

SEEP 

INPUT “ENTER MONTH, OATE ANO TIME (MM: O0:HH:MM:SS)° Oates 
OUTPUT 709;:°TO° Dates 

OUTPUT 709: °*TO° 

ENTER 709;:DateS 


PRINTER IS | 

PRINT 

PRINT ° Month, date and time: “:OateS$ 

PRINT 

PRINT USING °1@X,°° Fms akon Tev Vw" 2 Térop ~~ 

IF K#@ THEN 

PRINTER IS 701 

PRINT 

PRINT ° Month, date and time: “:DateS$ 

IF Itt#Q@ THEN PRINT USING “"10X,°°Tube Type: Wieland Smooth*** 
IF Itt#! THEN PRINT USING °1QX ,°°Tube Type: High Flux*°° 
IF Itt=#=2 THEN PRINT USING “10X,"°°Tube Type: Turbo-B°** 
PRINT 

PRINT USING °10X,°° Fas Tin Tev Vw*2 rdecp: -- 
PRINTER IS 1 

Ke] 

ENO IF 

Seer 

INPUT “ENTER FLOWMETER REAOING’ ,Fms 


OUTPUT 729;°AR AF® AL4 VRI°* 
FOR L=@ TO 4 
CUTEUT. 7091 As SA” 


IF L>@ ANO L<4 THEN 228262 
5=2 

FOR I=@ TO 9 

ENTER 729;:E 

S=S+E 

NEXT I 

IF L=@ THEN Emf(@)#A8BS( 5/12) 


IF Leda THEN Emf( 1 )#A8S(5/12) 
NEXT L 

OUTPUT 729;°AR AFQQ@ AL@S@ VRI“ 
OUTPUT 789;:°AS SA° 

Etp=2 

FOR I#@ TO 9 

ENTER 709;:Et 

EtpeEtp+Et 

NEXT I 

Etp=Etp/1Q 

TineFNTvsv(Emf(l )) 
TaveFNTvsviEaf(Q) ) 
Grad=37.9853+.104388eTin 
Mdot#=3.96S7E-3+Fmse(3.619SSE-3-Fmee(8.82Q006E-G-Fmse(|.23688E-7-Fmse4.31897 


E-10)))) 


20390 
204280 
20412 
22420 
20430 
20440 
20450 
20460 
20472 
20482 
22490 
22580 
20512 
20822 
20530 


VwemMdot/( 1@0@ePIe01°*2 )e4 

TdropeEtpe! .E+6/(1@e6rad ) 

PRINT USING °1@X ,3(00.00 ,4X%) 1X ,2.00,4X .MZ.40°: Fema Tin,Tev Vw'2 ,Tdroo 
BEEP 

INPUT “WANT TO ACCEPT THIS OATA SET? (18Y ,Q@=N)° Ok 


JeJj+| 

IF Ok#@ THEN 

Jej~} 

60TO 28130 

ELSE 

OUTPUT OFiler:Fms ,Emf(e) Etp 


PRINTER IS 721 

PRINT USING °1@X ,3(00.00 ,4X),1X,2.00,4X MZ. 40°F ms ,Tin,Tev Vw°2 ,Tdrop 
PRINTER IS | 

SEE 


a 


72540 
23550 
20560 
28572 
22582 
20599 
2Q@692 
tiles 

20612 
20620 
20632 
20640 
20652 
2@662 
20679 
22682 
2969 


2” 680O_gnal -4EN 70 
ENO IF 
ASSIGN @File TO ¢ 


PRINTER IS 791 
PRINT 
PRINT USING °10X,°°NOTE: °° ,22,°° Gate sets era stored in file °° ISA°: JF 


PRINTER IS | 
SU8END 


$U8 


Uoart 


PRINTER IS 1 
BEEP 

INPUT “Enter Uo File Name” Files 
BEEP 


INPUT 


"Number of Date Runs” ,Nrun 


INPUT “Oo You Want a Plot File7(1#Y ,O@eN)" plot 


22702 8EEP 

29710 IF Iplotel THEN 

20720 INPUT “Give Plot File Name” .P_files 
20730 CREATE BOAT P_file$,4 

20740 ASSIEN @Plot TO P_files 

20752 ENO IF 

29760 PRINTER IS 721 

20772 PRINT 

28782 PRINT 

20790 PRINT USING “1OxX,°° Water Vel Ue 
22800 ASSIGN @File TO Files 

20810 IF Iplotel THEN ASSIGN @Filel TO P_file$s 
229828 FOR [el TO Nrun 

20830 ENTER OFile:Vu Uo 

20840 IF Iplotel THEN DUTPUT @FileliVu Uo 
20850 PRINT USING *°15xX ,0.00,6X .MZ2.30E" (Vw Uo 
20862 NEXT I 

20870 ASSIGN @File TO e 

2Q@B8B@ ASSIGN @Filel TO « . 
208390 PRINT USING °18X,°°NOTE: °° ,22,°° data sets ere stored in file “~* ,1SA™INeu 
n,Files 

29988 IF Iplotel THEN 

2@91@ PRINT USING °10X,° "NOTE: °** 22,.°° X-Y¥ Pairs ere stored in file “° 1SA7=Nru 
nm Poti les 

28928 ENO IF 

20930 PRINTER IS 1 

20942 SUBEND 

20958 SUB Select 

22969 COM /Ido/ Ido 

22978 SEEP 

28982 PRINTER IS 1 

22990 PRINT USING “4X ,°"Select option:*** 


21202 
21012 
21028 
21030 
21049 
21258 


21260 
21287 
21082 
210990 
2110@ 
21118 
ee a 
21a 
21142 
21159 
21162 


PRINT USING “6X,°° 
PRINT USING °6X,°° 
PRINT USING °6X,°° 
PRINT USING “6X,°° 
PRINT USING “6X,°° 
PRINT USING “"6X,°" 
PRINT USING “6X,"" 
PRINT USING °6X,°° 
PRINT USING °6X,°° 
PRINT USING °6X,°" 
PRINT USING “6X,"" 
INPUT Idp 


is 
IF 
Le 
IF 
is 


Idp#@ THEN CALL 
Idpe#! THEN CALL 
Idp*2 THEN CALL 
Idpe3 THEN CALL 
Idp#4 THEN CALL 


@ Teking date or re-processing previous date’”" 
1 Plotting deta on Log-Log 
2 Plotting deta on Linear 
3 Mene cross-plot coefft file**” 
4 Re-circulete weter 

S Purge 

6 T-Orop correction’®** 
5 

8 

Z| 

Q 


Print Uo File’°** 
Modify X-Y file"”"” 


Move 
Comb/Fixup~ ~~ 


Mein 
Elo 
ets 
Coaef 
Main 


ot 16M 
2119 
21200 
21210 
21220 
21230 
21242 
21258 
21268 
21270 
21280 
21290 
21302 
21318 
21328 
21330 
21340 
21350 
21360 
21378 
2138 
21390 
21488 
21418 
21420 
21430 
21448 
21450 
21462 
21470 
21488 
21498 
21520 
21518 
21520 
21530 
21542 


3X3 f3 


To tepes. (MEN CAbe eurS 

IF lap=6 TREN CALL Tacn 

IF Ido#7 THEN CALL Uoprt 

IF Idp#8 THEN CALL Xymod 

IF Idp*9 THEN CALL Move 

IF Idp2#!8 THEN CALL Comb 

SUBEND 

SUB Xymod 

PRINTER IS 1 

BEEP 

INPUT “ENTER FILE NAME” ,File$ 
ASSIGN @Filel TO Files 

BEEP 

INPUT “ENTER NUMBER OF X-Y PAIRS” No 
BEEP 

INPUT “ENTER NEW FILE NAME* File2$s 
CREATE BOAT File2$,5 

ASSIEN @File2 TO File2$ 

BEEP 

INPUT “ENTER NUMBER OF X-Y PAIRS TO BE DELETED” ,Ndel 
IF Ndel#@ THEN 21410 

FOR Ie#!1 TO Ndel 

BEEP 

INPUT "ENTER DATA SET NUMBER TO BE DELETED” ,Ndi I) 
NEXT I 

FOR Je=1i TO No 

ENTER OFileliX,Y 

FOR Ie! TO Ndel 

IF Nd(I)@#J THEN 21498 

NEXT I 

OUTPUT OF ile2:x,Y 

PRINT J,xX,Y 

NEXT J 

PRINTER IS 701. 

ASSIGN @Filel TO 

ASSIGN @File2 TO 

SUBEND 

SUB Move 


21550! FILE NAME: MOVE 


21560! 


21578 


DIM Bop(66) ,A(66) ,8(66) ,C( 66) ,0(66) ,E(66) ,F( 66) ,6(66) ,H( 66) ,J(66) K( 66) Lé 


66) ,M(66) 


21588 
21598 
21602 
21618 
21628 
21630 
21648 
21658 
(I) 
21668 
21678 
21688 
21698 
21780 
21710 
21720 
21738 
21748 
21758 
21760 
21778 
21782 
NCI) 
21798 


DIM Told$(66 (14) ,.N( 66) ,Vr(66) ,Ir(66) 
BEEP 

INPUT *OLO FILE TO MOVE* ,02_files 
ASSIEN @File2 TO DO2_files 

ENTER OF ile2:Nrun ,DateS Ldtc! Ldtc2 Itt 
FOR I*!1 TO Nrun 

ENTER @File2:Boo( 1) ,Told8(I) 


ENTER @File2:A(I) 8¢1),C¢1) OCT) ECT) FCI) 6¢1) ACT) JCD) KOT) UCT) MCI) N 
ENTER @File2:Vr( 1) ,Iril) 

NEXT I 

ASSIGN @File2 TO e 

BEEP 

INPUT “SHIFT OISK AND HIT CONTINUE® ,Ok 

BEEP 


INPUT “INPUT 8DAT SIZE" Size 

CREATE BDAT D2_fileS Size 

ASSIEN @Filel TO O2_files 

OUTPUT @FileliNrun ,Date$ .Ldtc!l Ldtc2 Itt 

FOR I! TO Nrun 

OUTPUT @Filel;:Boo( I) ,Told8(I) 

OUTPUTCORT ai sAtl > (1) CCl) OCL) ECL) FCl) G¢2f) MCL) .J¢1T) KOT) LCL) MCD), 


OUTPUT @FiletsVr( I) ,Ir(l) 
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2'8@@ NEXT ! 

21812 ASSIGN BFilel TO e 

21820! RENAME °TEST® TO O2_filas 

21830 BEEP 2000,.2 

21840 BEEP 4000,.2 

21850 8EEP 4808,.2 

2185@ PRINT “OATA FILE MOVEQ”’ 

2187@ SUBEND 

21888 SUB Comb 

21892! FILE NAME: COM8 

2190Q! 

21918 OIM Emf( 12) 

2192@ SEEP 

2193@ INPUT “OLD FILE TO FIXUP” O2_files 
21948 ASSIGN @File2 TO O2_files 

2195@ Oi_frileS=°TEST" 

2196@ CREATE 8DAT O!_file$,22 

21970 ASSIGN O6Filal TO Oi_filas 

21398@ ENTER @File2itNrun ,OateS i.dtcl ,.Ldte2 ,Itt 
21990 Nrunme@22 

22000 IF K#@ THEN QUTPUT @FilelsNrunm DataS ,L_dtc! Ldatc2 Itt 
22012 FOR Te! TQ Nrun 

22020 ENTER O@FileZ2:80p ,Told$ ,Emfle) Ve Ir 
22038 OUTPUT OFilai:8o0p ,ToldS Emf(e) Vr Ir 
22040 NEXT I 

22050 ASSIGN OFile2 TO « 

22062! RENAME “TEST’ TO O2_files 

22070 BEEP 2000,.2 

22Q82 BEEP 4020,.2 

22092 BSEEP 4000,.2 

2210@ 8SEEP 

22118 INPUT “WANT TO AQO ANOTHER FILE (1@Y Q9N)?" Oka 
22120 IF Okael THEN 

c213@ Ke} 

22148 8EEP : 

2215@ INPUT “6IVE NEW FILE NAME* Nfilas 
22168 ASSIGN OFile2 TO Nfilas 

22178 GOTO 21980 

2218@ ENO IF 

22198 ASSIGN OFile2 TO « 

22200 SUBENO 


Oata : 19 Jan 1993 


NOTE: Program nama : ORP6B 
Or1sk numbar = QQ 
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